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TORQUE EQUILIBRIUM ATTITUDE CONTROL 
FOR SKYLAB REENTRY 

INTRODUCTION 


With the decision in December 1978 to discontinue efforts to keep 
the Skylab Space Station (Fig. 1) in orbit came the decision to terminate 
the highly successful end-on-velocity- vector (EOVV) mode of operation 
11] and to reestablish the solar inertial (SI) mode as soon as practical. 
This would put the spacecraft in an ideal power attitude and greatly 
reduce ground management of systems. It would also increase the drag 
on the spacecraft and cause it to reenter sooner. Unfortunately the SI 
attitude was not going to be maintainable much below 280 km (150 n.mi.) 
since the growing density of the atmosphere was going to cause the aero- 
dynamics to grow to the point where the storage capacity of the control 
moment gyroscopes (CMG’s) for angular momentum was inadequate. Since 
control of attitude to 150 km or below would be required to be able to 
influence reentry and the amount of thruster gas was far too low to 
consider control with thrusters, another new attitude control scheme had 
to be developed. Aerodynamic torques are proportional to density and 
become nearly overwhelming at 150 km, and therefore any attitude control 
scheme which would work at 150 km would have to accurately take into 
account the aerodynamic disturbances. Thus it became clear that we 
must look at the torque equilibrium attitudes (TEA's) if any existed and 
plan our control schemes about these. Thus the important early ques- 
tions were: 

1) Are there equilibrium attitudes? 

2) If so, is adequate solar power available? 

An investigation with the mathematical aerodynamic model of the 
Skylab vehicle indicated there were no aerodynamic trim or equilibrium 
attitudes but there may be attitudes where aerodynamic, gravity gradient, 
and gyroscopic torques balance. Indeed, 12 such attitudes were found. 
Most were not useable as control attitudes since there would be insuffi- 
cient solar energy available to power the spacecraft and battery power 
would be completely inadequate for the several weeks required. Only 3 
of the 12 TEA’s appeared viable if a control scheme could be developed 
for them. So work began on a TEA control scheme. A candidate scheme 
which was promising was developed, but there were large uncertainties 





because of lack of confidence in the aerodynamic coefficients. This lack 
of confidence was due to the lack of test data confirming the aero moment 
model. As a result we did not have high confidence that the new TEA 
scheme would work. Later flight performance would show that our models 
were much better than we expected. The various mathematical models 
and tools needed for TEA control are described here. The operational 
aspects (time lines, command history, etc.) of TEA control as well as 
more complete background information and descriptions of the other 
essential Skylab systems (power, telemetry, etc.) are contained in 
Reference 2. 


COORDINATE SYSTEMS 


The coordinate systems which are pertinent to Skylab TEA control 
are defined in this section. Each system has some special geometrical or 
physicul feature which simplifies the solution of a particular problem. 

The following coordinate systems are described in this section: 
Orbital, Vehicle, Attitude Reference, Solar Inertial, and Z-Local Verticul. 
Each coordinate system consists of a set of mutually orthogonul uxes 
exhibiting right-handedness. 

An inertial (with respect to rotation only) coordinate system is a 
system which retains its orientation with respect to the celestial sphere, 
although the origin may be moving along any general curvilinear path in 
space. Similarly, a vehicle fixed system retains its orientation with 
respect to the vehicle. 

Orbital Coordinate System (0) 

The Orbitul Coordinate System (X Q , Y o < Z q ) is a precessing coor- 
dinate system with its origin at the Earth center of mass. The rate of 
precession about the Earth's north pole is approximately -5 <' .ees/day. 

The Z axis lies in the orbital plane, positive through the ascending 
o 

node of the orbit. The X q uxis also lies in the orbital plane 90 degrees 

ahead of the Z q axis. Since the Skylab orbit was in the X q Z o plane at 

all times, the Y axis was parallel to the orbital ungular momentum vector 
o 

completing the right-handed system (Fig. 2). 


% 




SCUPTIC 
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Figure 2. Coordinate systems. 


Vehicle Coordinate System (V) 

The Vehicle Coordinate System (X y , Y y , Z y ) is a vehicle-fixed 

system with its origin at the center of mass. The X y axis lies along the 

long axis of Skylab and is positive in the direction of the Multiple Docking 
Adapter (MDA) . The Z y axis is positive toward the Apollo Telescope 

Mount (ATM) and the Y y axis completes the right-handed system. 

Attitude Reference Coordinate System (A) 

The Attitude Reference Coordinate System (X A , Y., Z A ) is a 

AAA 

movable system with its origin coincident with the Vehicle Coordinate 
System origin. The axes of this system represent the instantaneous 
desired orientation of the Vehicle Coordinate System axes. 
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Solar Inertial Coordinate System (I) 

The Solar Inertial Coordinate System (Xj, Yj , Zj) is only a pseudo- 

inertial system since it makes one revolution per year. !t was used 
during the Skylab mission to point the instrument? in the desired direc- 
tion. The origin is coincident with the origin of the Vehicle Coordinate 
System origin. The 7,. axis is positive toward the center of the Sun. 

The Xj axis lies at an angle v z from the orbital plane (this angle is cal- 
culated on-board such that the principal X axis is in the orbital plane 
to minimize the build-up of angular momentum) and is positive toward 
the sunset terminator. 


Z-local Vertical Coordinate System (L) 

The Z- Local Vertical Coordinate System (X. , Y. , Z. ) is a rotating 

L \j 

system with its origin at the center of mass of Skylab (the rate of rota- 
tion is one revolution per orbit). The X axis is positive in the dircc- 

I i 

tion of flight and lies in the orbital plunc. The L, uxis is parallel to 

the local vertical direction and is positive outward, away from the Earth. 
The Y. uxis is purallcl to the orbit normal and is positive towurd orbital 
North. 


SKYLAB ATTITUDE AND POINTING CONTROL 
SYSTEM (A PCS) 


The actual control of the Skylab attitude to the attitude reference 
was done exactly as in the original mission [31. However, only the 
pointing control system (PCS) of the APCS was used; the experiment 
pointing control system (EPCS) was disabled. The major parts of the 
APCS were the rate gyros, the sun sensors, the star tracker (it had 
failed during the original Skylab mission), the Apollo Telescope Mount 
Digital Computer (ATM DC), the Workshop Computer Interlace Unit 
(WCIU). double-gim baled CMC's, and cold-gas (compressed nitrogen) 
Thruster Attitude Control System (TACS). 

Six control modes were addressable: (1) STANDHY. (2). SOLAK 

INERTIAL (SI), (3) EXPERIMENT POINTING, (4) ATTITUDE HOLD /CMC , 
(5) ATTITUDE HOLD /TACS, and (6) ZLV. TEA control was programmed 
to be a substate of the ZLV mode. The basic ZLV' (for Z axis along the 
local vertical) attitude was with the Z axis along the local vertical, point- 
ing up, and the X axis in the orbital plune. pointing in the direction of 
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the velocity vector. Any angular offset from the basic ZLV attitude 
(offset identified by the quaternion Q AL ) could be commanded via a set 

of three Euler angles (x) with a 2-3-1 rotation sequence. Changes in 
the TEA attitude reference were achieved by commanding changes to Q^> 

tAL1 ■ lx x l x CX.1, 'x y i 2 - 1 


Since none of the original APCS capabilities were eliminated by the 
addition of the TEA control, the torque equilibrium attitude control could 
also be held by TACS only (this option was contemplated for low-drag 
attitude at low altitudes in case low-drag was needed for reentry control, 
but active TEA control had proven impossible). 

Tile TACS consisted of software plus six cold-gas thrusters, two 
uncoupled ones for control about the Y axis (one would be fired for +Y, 
the other for -Y control); four other thrusters for coupled X and Z 
control [4j. 


CMG CONTROL SYSTEM 


The CMG control system was composed of three orthogonally 
mounted, double-gimbaled CMG’s with angular momentum magnitude H of 
3050 Nms (2280 ft-lb-sec) as shown in Figure 3. The CMG control law 
utilized three normalized torque commands and thr> CMG momentum status 
to generate proper CMG gimbal rate commands [5], The CMG control law 
consisted of three parts; CMG steering law, rotation law, and gimbal 
stop avoidance logic. There also were some other routines for specialized 
situations like caging the CMG's to a desired momentum state [4]. 


1. [ ]-quantities are usually 3x3 rotation matrices, with the further 

definitions (where Xj represents any angle; s = sin, c = cos); 
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Figure 3. CMG mounting arrangement. 


The CMG control law had the ability to operate with either three 
or two CMG’s for redundancy. Since CMG Mo. 1 had failed during the 
original Skylab mission, the CMG control law was always in the two-CMG 
option. 


BANK, ATTACK, AND ROLL ANGLES 


The aerodynamic forces and torques are dependent on the aero- 
dynamic angle of attack and roll angles. To complete the specification of 
an arbitrary attitude, a third rotation angle is required. This angle is 
the so-called bank angle. The sequence is 
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[VL] 


t Pl x [a] 2 t ei x 


( 1 ) 


These angles are obtained from the following relations; 
a = arc cos (VL^) 

6 = arc tan (VL 12 /-VL 13 )| 


four quadrant 


( 2 ) 

(3) 


p = arc tan (VL 21 /VL 31 > 


(4) 


The Skylab onboard computer works in quaternions so that the matrix 
[VL1 must be computed from the given quaternions Qyj* Qy A » ^AL : 





(5) 


The matrix [VL] is derived from Q„. by taking the upper left 3><3 
matrix from 


[VL] = (upper left 3 x 3 of) 


Q - ^ O 

S'L ^VL 


( 6 ) 


The physical bank, attack, and roll (BAR) angles will differ 
slightly from those above, because the orbit of Skylab is not truly cir- 
cular and the air tends to move with the Earth surface so that the 
velocity of Skylab relative to the air is not exactly just the orbital 
velocity. Thus when equations (A2) and (A3) are used, values slightly 
different from those in equations (2) and (3) must be used. These 
values were only used for evaluation purposes and proved to be quite 
useful, so the differences were insignificant. 


2. See Appendix E for nomenclature. 
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TORQUE EQUILIBRIUM ATTITUDES 


As mentioned previously the aerodynamic coefficients are only a 
function of the roll angle and the angle-of-attack. Figure 4 shows the 
zero-moment curves for the three components with respect to these 
angles (see Table 1 for mapping of the vehicle axes). As can be seen, 
there is no set of angles, where all three curves intersect. However, 
they come close in several areas. To get a true three-axis equilibrium, 
other external torques are required. Gravity gradient (GG) torques and 
gyroscopic torques were found to be sufficiently large for all altitudes 
of concern to create 12 TEA’S altogether. Gyroscopic torques come into 
play, since the aerodynamic and GG torques are constant with respect 
to the rotating local vertical coordinate system, and the total angular 
momentum of the Skylab was selectable within certain limits (imposed by 
the finite storage capacity of the CMG system). Since the aerodynamic 
torques do not change (relative to body-fixed axes) when the Skylab is 
rotated about the relative wind velocity vector and GG torques do not 
change in body-fixed axes for a rotation of 180 degrees about any axis 
perpendicular to the local vertical, there is always a pair of TEA'S with 
the same angle-of-attack /roll angle combination, but with bank angles 
differing by 180 degrees, i.e. , there were actually only six basically 
different TEAs with respect to the aerodynamic torques (they are 
indicated in Fig. 4 by asterisks). The TEA’S are shown for an altitude 
of 200 km (108 n.mi.) and zero total angular momentum (Table 2). 

From solar panel power considerations a sun-pointing, inertially 
fixed attitude, as given in the SI mode, is the best and TEA is rather 
bad. Nine of the TEAs were completely hopeless when it was established 
that 28 percent of full Sun illumination was required, on the average, to 
supply the needed power ( 100 percent is the power received when the 
Sun is perpendicular to the solar panels and the vehicle is in an all- 
daylight orbit, as is the case for high Sun elevation angles with respect 
to the orbital plane). Only the remaining three TEA'S were usable, and 
these only when the optimum angular momentum was used (each of the 
TEA’S exists in a volume of the BAR angle space when the total angular 
momentum is varied within its available volume). Even then, some of the 
TEAs did not have enough power for certain Sun elevation angles. The 
three TEA'S were named T121P, T121G, and T275, and they are shown 
in Figures 5, 6, and 7 where the point of view is slightly south of the 
orbital plane and the vehicle is moving from the front lower left to the 
back upper right. T275 was a low-drag attitude (the ballistic coefficient 
was approximately 275) with the MDA trailing; T121G and T121P were 
high-drag attitudes (the ballistic coefficients were about 121). In T121G 
the Skylab was approximately in a GG equilibrium attitude with the MDA 
pointing upward (the same attitude which Skylab had been left in when 
it was deactivated in early 1974); in T121P the MDA was almost perpen- 
dicular to the orbital plane and pointing South. For T121G and T121P 
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Figure 4. Combined zero torque contours. 
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TABLE 1. MAPPING OF VEHICLE AXES 


Axis 

Angle of Attack (deg) 

Roll Angle (deg) 

+x 

0 

any 

-x 

180 

any 

+y 

90 

90 

-y 

90 

-90 

+z 

90 

0 

-z 

90 

±180 


TABLE 2. TEA’S FOR 6 = 3.11E-10 kg/m 3 
(-200 km or 108 n.mi.) 




Bank 

Angle of 

Roll 





Angle 

Attack 

Angle 



No. 

ID 

(deg) 

(deg) 

(deg) 

CD 

BC 


T121P/SU 

ND 

-78.6 

+101.4 

106.0 

+124.9 

6.66 

132.5 


T121G/UN 

DS 

-173.2 

+6.8 

93.8 

+124.0 

6.94 

127.3 

5 

6 

T275/BU 

BD 

+70.2 

-109.8 

147.8 

1-89.2 

3.52 

250.5 

7 

8 

BN 

BS 

-51.4 

+128.6 

168.1 

-117.6 

2.72 

324.9 

9 

10 

DN 

US 

-11.9 

+168.1 

110.6 

-60.7 

6.37 

138.9 

11 

12 

NU 

SD 

+86.4 

-93.6 

71.8 

-51.2 

7.11 

124.2 


Noti: These TFA's are for zero angular momentum and a CP offset 

T 

of [+0.06, -0.07, -0.11] m. See Appendix G for 
nomenclature. 


11 























UP 



Figure 5. (Concluded). 
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Figure ?. T275 (no angular momentum). 


the solar panels were trailing and they were statically stable with respect 
to aerodynamic torques. T121G was also statically stable with respect to 
the GG torques (therefore stable in all axes) whereas T121P was in an 
unstable GG equilibrium. 


TEA SEEKING METHOD 


The assumption is made that the total external torque, T acting 
on Skylab is changing linearly with the attitude offset, and that 

the partial derivatives of the torque with respect to the offset are known. 

-ext " [Sf] &4att 3 (7 > 


3. An underlined quantity (&* is a vector with three components. 


where 



Evaluation of the on-board total angular momentum change over the 
desaturation interval, T . , results in an estimate of the total external 

torques, des 


Text - <« ' Hp>' T des 


( 8 ) 


where H is the present total angular momentum and is its past value. 

The attitude offset from the torque equilibrium attitude, assuming the 
offset is constant, is then 



( 9 ) 


Changing the attitude reference by -A 4> Q ff would, ideally, eliminate a 

further angular momentum change. However, the previously accumulated 
angular momentum away from a desired momentum state, H , has to be 
eliminated during the next desaturation interval by "~ nom 


Ho 


om 


" 3T il 1 

- - <» - » 


— nom des 


( 10 ) 


The total required attitude change is therefore 
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A± « Aimon, “ Aioff 

A! = rf?~ — 1"^~1 [(H - H ) + (H - H )] . (11) 

T des L^lJ ~ ” n ° m ~ “** 

This method for attitude change eliminates, ideally, any initial condition 
within two desaturation intervals. [The parenthetical expressions in 
equation (11) were not combined since these quantities had to be limited 
separately, ] 

In the Skylab software, all angular momentum quantities were 
normalized by the nominal angular momentum magnitude, H, of one CMG, 
and they were called e. Equation (11) then becomes 


A4> = [SLOPE! (Ae + Ae - Ae ) 
— — m — — p 


( 12 ) 


where 


[SLOPE] = 


K‘ 

des 


-1 


Ae = (H - H )/H 
— — — nom 


Ae_ = the limited value of Ae 
— m - — 


Ae^ = the past unlimited value of Ae.. 


The reorientation capability of Skylab was limited and therefore 

M had to be limited also. To avoid the possibility of a large momentum 
set overcoming the signal due to a momentum change, Ae^ is limited 

to a value which cannot command more than approximately 80 percent of 
the limit on AJ>_. 

The actual reference change is done by generating a quaternion 
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AQ al = [A^/2, / 1 - 0*25 A£ 2 ] T 

4 

and updating the reference quaternion 
< W new) * ^AL Q AL (oW) • 


STRAPDOWN DRIFT CORRECTION 


(13) 


(14) 


The TEA control mode contrins a provision for strapdown updating 
from the ground. This capability was used extensively to maintai TEA 
control. The error in the strapdown was estimated by several te uniques 
which gave either the strapdown error components in the orbit p 1 \ne or 
along the orbit normal. The only source of data for strapdown error 
about the orbit normal (or orbital Y axis) is data from the Sun sensor 
or the solar arrays. These data are available only occasionally from the 
Sun sensor when its line of sight along the vehicle Z axis comes within 
approximately 18 degrees of the Sun over a ground station. Strapdown 
error along orbit normal was estimated from the timing error of solar 
passage through the vehicle XZ plane. This event happened once per 
orbit on the daylight side. Passage occurred when the X Sun sensor 
output passed through null. Null passage was too indistinct to be 
reliable when the Sun was greater than 18 degrees away from the vehicle 
Z axis at passage. For those periods, the solar array power angle or 
the aerodynamic roll angle deviation was used to estimate the strapdown 
error. These latter techniques were less reliable and were used only 
when a consistent trend was established or when more than one indicator 
was telling the same story. 


The strapdown is the onboard vehicle attitude reference. The 
vehicle angular rates are integrated using quaternions as attitude 
parameters. Since the rate data contains errors which are integrated 
into these quaternions, they must be corrected periodically to eliminate 
offsets from the true reference. Such corrections are referred to as 
updates. These strapdown updates are applied through an update 
quaternion which is quaternion multiplied onto the strapdown quaternion. 
The vehicle strapdown quaternion is called Qyj* and the update 


quaternion is 
angle A_0. : 



This quaternion is constructed from the vector 



[0.5fVL]A0 L , 


/ 1-0.25A9 


2 jT 


(15) 


4. See Appendix E for nomenclature. 
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The quaternion update multiplication is performed as shown 


Q VI (new) « ASyj Qyjfold) * 


( 16 ) 


The onboard computer computes the quaternion and from this comes 

«AI - 5 al «J u tm 

and 

Q VA = *vi Q AI * (l8) 


The quaternion Q AL is the commanded attitude of the vehicle relative to 

the local vertical reference. The quaternion Qy A is the attitude error 

quaternion which is maintained close to identity by the onboard control 
system, A method for estimating the strapdown correction angle Ad is 
given in Appendix D. 

Originally, an automatic on-board in-the-orbital-plane strapdown 
update (with only updates perpendicular to the orbital plane from the 
ground) had been considered. It was assumed that the active TEA con- 
trol had brought the V system into alignment with the TEA at the start 
of the desaturation interval. At the end of the next desaturation 
interval any misalignment about the axis then results in a first order 

offset from the ideal TEA about the axis. Using part of the A 4> 

correction about the Z^ axis for a strapdown update should bring the 

estimate of the orbital Y^ axis into alignment with the actual orbital 

axis, since the misalignment is inertially fixed and the misalignment 

sweeps the whole orbital plane. This scheme worked as long as the 
assumption was made that the relative wind direction is always co!i»,?ar 
with the vehicle velocity vector. However, the atmosphere rotate v/ith 
the Earth and therefore the relative wind direction oscillates approximately 
two degrees about the vehicle velocity vector. The result was that the 
onboard strapdown update settled out at an offset which was furthermore 
dependent upon the TEA used. Since the ground update method (which 
did not have these shortcomings) was developed at that time, the onboard 
update was disabled by setting its gain to zero. (There was not enough 
time available to take it completely out of the ATMDC software [2]. 
Therefore the ground update method was used through the already exist- 
ing ground command capability, eliminating the need for any further 
software changes). 
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TEA OPERATION AND PERFORMANCE 


Nominal flight parameters had to be continuously generated for 
several days in advance since the conditions tended to change (for $ 
example: estimated density, estimated CP location, etc.). The data 

were given to the flight controllers in the form of tables (Table 3^ shows 
the data for DOY 171.5 through DOY 174.5 or 6/20/79 to 6/23/79 at noon 
GMT). Since the altitude decreased slowly and everybody needed time to 
get used to the TEA operation, the SLOPE was not changed during this 
time. ETLN (with components ETLN1 , ETLN2, ETLN3 in the L system) 
is the normalized nominal angular momentum. Z-AXIS BETA is the 
elevation angle of the vehicle Z axis with respect to the orbital plane. 

CD and BC are the drag and the ballistic coefficients, respectively. 

CMG MOM is the angular momentum in the CMG system and it is in per- 
cent of 3H. ETSF is the normalized and filtered total system momentum 
in vehicle components. QBLNOM is the nominal attitude reference 
quaternion (a QBL was existing already and could be used; there- 
fore, the subscrips A and B are equivalent here). 

Table 4 shows the Slope Matrix Update Log. The first three 
columns show the date when the Slope was updated (DOY 171 corresponds 
to 6/20 and DOY 192 to 7/11). The first entry is the start of TEA con- 
trol rather than the time when the initial Slope was loaded. The index 
in the fourth column is added to facilitate correlation with Table 5. 

Columns 5 through 13 show the components of the Slope as indicated by 
the heading. Table 5 shows what data wore used to generate the slopes. 
The slope generation parameters for index 2 and 3 are the same; the 
transpose of Slope 3 had been sent up (Slope 2) by mistake, whereupon 
control was promptly lost. Slope 3 was sent as soon as attitude control 
had been regained. Slope 13 was used for a higher altitude than the 
one for which it was calculated, since a 6 hr gap in the gvound coverage 
eliminated an additional slope change and it was considered better to have 
the proper slope at a lower altitude and take a reduced gain at the higher 
altitudes. 

Figure 8 shows the actual BAR angles for DOY 172:03:21 to DOY 
172:06: 12 versus time in seconds and the predicted BAR angles are 
shown as horizontal lines. It can also be noted that the BAR angle 
traces zig-zagged, indicating that the gain was too high. As a conse- 
quence the slope gain was later reduced to 0.5 from 1.0, resulting in a 
much smoother trace. The actual BAR angles are only correct to within 
the strapdown error, which, due to the availability of acquisition Sun 
Sensor (ACQ SS) information, was less than 1 degree. A very good 
correspondence between the nominal and the actual BAR angles can be 
seen, indicating that the prediction of the CM und CP locations, as well 


5. In the APL computer language, a minus sign is a superscript minus 
(to distinguish it from the subtraction operation). 
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TABLE 3. NOMINAL FLIGHT PARAMETERS 
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171.50 141.97 10.20 0.0 0.0 *1.0 0.70 6.55 154.72 

172.50 140.03 13.40 0.0 0.0 *1.0 0.65 6.55 134.79 

173.50 140.11 0.60 0.0 0.0 *1.0 0.71 6.56 134.60 

174.50 139.13 4.10 0.0 0.0 *1.0 0.01 6.56 134.51 


TABLE 3. (Concluded) 
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TABLE 4. SLOPE MATRIX UPDATE LOG 
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TABLE 5. SLOPE GENERATION PARAMETERS 
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Figure 8. Actual BAR angles (DOY 172:03:21/172:06:12). 
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as the aerodynamic coefficients, was much better than anticipated. The 
nominal BAR angles were generated assuming that the relative wind 
velocity is antiparallel to the vehicle velocity vector, but reduced by an 
amount appropriate for the Earth rate. The nominal BAR angles are 
therefore the nominal average of the actual BAR angles. Some of the 
good correspondence is due to the fact that changes in ETLN1 were used 
to influence the roll angle. Initially, a roll angle smaller than nominal 
was very detrimental since the slope matrix changed drastically such that 
a roll angle of less than 90 degrees was found to be unstable: The 

actual slope components had changed sign. An initial unfavorable rate 
gyro bias could have changed the actual roll angle by that much long 
before the ground would have had enough data to detect and correct for 
it. To eliminate this possibility, a favorable rate gyro bias compensation 
was introduced into the software such that if we had had the bad gyro 
drift, the compensation would have eliminated it. As it turned out, the 
real rate gyro bias was in a good direction and the compensation had to 
be changed to the opposite polarity (as always, Skylab seemed to refute 
Murphy's Law, since anything that could go right, would; also see the 
last paragraph of this section). The great sensitivity to the rate-gyro 
introduced drift about orbital North was due to the relatively large 
negative angular momentum bias in the Z axis (ETLN3 = -1H) which then 
was resolved into the actual axis thereby changing the equilibrium 

position drastically. This resolution also suggested the (at least 
temporary) remedy: Buck it with an appropriate amount of momentum 

bias in the X T axis. 

Li 

An overview of the BAR angles for the three weeks of TEA opera- 
tion is given in Figure 9. The bank angle is the lowest trace , the 
angle of attack is in the middle, and the roll angle is the highest. The 
horizontal scale is in days, the numbers indicate the start of a day. 

Only BAR angles from the auxiliary storage and playback recorder (ASAP) 
tapes are shown, since they were saved in our ground computer. 

TEA control was initiated on DOY 171:13:17 GMT (6/20/79) and 
Skylab impact occurred on DOY 192:16:37 GMT (7/11/79). During the 
latter half of DOY 171, it can be seen that the average roll angle steadily 
increased. This was due to the initial rate gyro bias (intentional com- 
pensation and basic bias). After corrective action was taken the angles 
settled close to their predicted values. Figures 9a and 9b show the BAR 
angles during the shift in ETLN3 from -1H to 0 between DOY 175.5 and 
DOY 180.5. With the gradual removal of the ETLN3 the sensitivity to 
rate gyro drift was removed; however, the excellent control capability 
over the roll angle (and with it, the drag) also vanished. In fact, the 
actual roll angle could later not be influenced at all (the indicated roll 
angle also showed the effect of strapdown error about orbit North and 
this fact in conjunction with the fact that the actual roll angle was 
steady, allowed the use of the indicated roll angle for strapdown update 
about orbit North). 
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Second week. 
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To illustrate the evaluation of the information on the ASAP tapes, 
DOY 187:16:56:30 to DOY 187:20:56:30 are taken as an example. Table 
6 shows the basic QBL/time (QBL = Q^ L ) information, the equivalent 

BAR angles are shown in Table 7 and graphed in Figure 10. The BAR 
angles were always evaluated for their average and a sinusoid of orbital 
frequency. The latter was fitted by a least-square fit, after the aver- 
age had been subtracted. The results for our example are shown in 
Figures 11, 12, and 13. In the lower left quadrant of each figure there 
is a table of expected angle values (right column) as a function of the 
orbital angle (left column). Since the circle charts (Fig. 14) show this 
orbital angle, the tables were used to check the actual BAR angles (live; 
from telemetry) against these reference angles. 

The sinusoids in the angle-of-attack and roll angle were due to the 
relative wind direction oscillating about the vehicle velocity direction 
(since the atmosphere rotates with the Earth). The sinusoid in the bank 
angle was due to the in the orbital-plane strapdown error and the appro- 
priate strapdown update is shown in the upper left hand corner of 
Figure 11, both with respect to orbital midnight and with respect to the 
actual time; the update is to be telemetered. For the latter a resolution 
of the error with respect to midnight had to be made. 

Application of the various methods was severely hampered by lack 
of data and large delays between the occurrence and the receipt of the 
data. If the strapdown error had drifted substantially there would have 
been no way to correct it in time to prevent loss of control. Fortunately, 
the rate gyro drift components in the orbital plane (if constant) integrate 
to exactly zero in one revolution and generally TEA control was not very 
sensitive to the drift components along the orbit normal. In-plane strap- 
down errors resulted from inaccurate navigation updates and from the 
Sun motion of about 1 degree /day. For data on strapdown errors, the 
onboard ASAP tape recorder had to be run for at least one full orbit 
(more consecutive orbits were desired for noise content reduction) during 
a quiet state; i.e. , the last slope update had to be done two orbits before 
and no other strapdown correction was permitted. Once the ASAP tape 
was dumped over a ground station, it took at least another half hour of 
data processing and parameter extraction. The results were 18 sets of 
Q al P er orbit, which then had to be fed into our Sigma 5 computer to 

be evaluated with the appropriate programs. The actual updating of the 
parameters had to wait until the next ground station. The bottom line 
was, that the delay in the actual onboard happening and the corrective 
action amounted to at least two orbits, sometimes more. 

The ground support operations were aided by the "circle charts" 
(Fig. 14) and telemetry formats (Format 8 for TEA parameters is shown 
in Fig. 15). The data on Format 8 changed every second provided the 
Skylab was over a ground station; the circle charts were applicable for 
one orbit. 


29 


TABLE 6. REPRESENTATIVE QBL SET (DOY 186:16:56/186:20:56) 
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Figure 10. Act 
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Figure 11. Btlnk angle evaluation for QBL18716 
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Figure 12. Angle-of-attack evaluation for QBL18716. 
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Figure 13. Roll angle evaluation for QBL18716 


STATION ACQUISITION PERIOD? FOR SKVLhB 

TIME OF NORAD ELEMENT SET 18:09(17 GMT 
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The circle chart has midnight at the 6 o'clock position; the outer 
numbers are minutes from midnight; the inner ones are orbital degrees 
from midnight. Orbital night time is indicated by the heavy part of the 
circle. The boxes on the outside of the circle indicate the ground 
station coverage; the ground stations are Ascension (ACN), Bermuda 
(BDA) , Goldstone (GDS), Madrid (MAD), and Santiago (AGO). 

Format 8 showed the present and the past Ae, the present attitude 
reference change AJ>_, and any impending strapdown update A 0. Also 
shown were the filtered total system momentum (in units of ffT, the CMG 
momentum (in percent of 3H>, and the inner and outer gimbal angles (to 
the nearest degree) plus then gimbal rates. 

Of further interest are the quaternions Q AL (QBL), Q VA * and 

Q VI - The first allowed evaluation of the current BAR angles, the second 

showed how well the vehicle was following the attitude reference, and the 
third indicated the attitude with respect to the SI system (useful for 
determining the closest approach of the Sun line to the center of the 
ACQ SS which then was used for strapdown update information). 

The history of the TEA control parameters as well as the solar 
elevation angle with respect to the orbital plane are shown in Figures 
16 through 20. The biasing history of the rate gyros is shown in 
Figure 18. Only one rate gyro was still operating in the Y axis (Y3), 
two each were averaged for the other two axes (XI and X3; Z1 and Z2). 
As mentioned before (cf . middle of 3rd paragraph of this section) , the 
initially introduced beneficial bias had to be taken out and replaced with 
a bias of the opposite sign. The fact that in two axes two rate gyro 
outputs were averaged allowed us to work effectively with half least 
significant bits (LSBs) by biasing only one of the averaged gyros with 
a full LSB. A trimming of the rate gyros to within about 0.25 degree/ 
orbit drift was therefore possible. 

The relatively many changes in ETLN1 (Fig. 19) reflect our effort 
in controlling the roll angle to acceptable values. ETLN2 did not 
influence the TEA since it was parallel to the orbital velocity, but it 
affected the CMG gimbal angles and was therefore changed once (in two 
steps, to minimize the transients) to improve the actual gimbal angle 
traces . 

Simulation predicted that TEA control in T121P was possible down 
to about 140 km (75 n.mi.) and in T275 down to 130 km (70 n.mi.). At 
these altitudes the GG and precessional torques were no longer large 
enough to achieve zero external torques in all axes simultaneously. Very 
fortunately, Skylab was about to reenter on the most desirable orbit, 
one with the least population density beneath it. This orbit, however, 
did go over the United States, the southern part of Canada and over 
Australia, but most of it was over water. Therefore, only a slight 
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Figure 16. Sun elevation angle beta (April 1978/August 1979). 
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Figure 18. History of rate gyro biasing 
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adjustment of the nominal impact point (which was in the Atlantic Ocean) 
was indicated, and there was no need to maneuver to the low-drag T275 
attitude (this attitude was on an unstable equilibrium in all axes, and 
the equilibria were much more precarious than the ones in T121P). The 
slight impact point adjustment could be achieved by going to a random 
tumble (with somewhat less average drag than T121P) at an altitude of 
150 km (81 n. mi.). This was done by commanding a delayed maneuver 
from the T121P attitude to the solar inertial attitude. The latter, due 
to the high aerodynamic torques, could not be held, and a random 
tumbling resulted. On the average this random tumbling had the 
desired drag, which then lengthened the lifetime sufficiently to place 
the nominal impact at the desired location. The vehicle broke apart 
later and, therefore, lower than expected. As a consequence the drag 
was reduced, shifting the reentry point farther downrange. 


CONCLUSIONS 

The TEA control was most successful in performing the mission 
that was originally laid out for it. The experience gained from the TEA 
control effort in particular and the Skylab Reactivation in general will 
be valuable in all future space activities. We found out that we could, 
in fairly short order, come up with new control methods, new procedures, 
and new tools to perform tasks and accomplish goals that would have been 
considered impossible prior to the original mission. However, this fast 
reaction capability was due to the fact that we were able to use APL on 
a Sigma V computer with interactive terminals. We can flatly state that 
without APL the development of the new Skylab control methods (EOVV 
and TEA control) would not have been possible in the extremely short 
time available (EOVV control II] was conceived, developed, simulated, 
implemented, and flown in less than three months). In many respects 
TEA control was more successful and went more smoothly than the pre- 
vious phase with EOVV. This may have been due in part to the greater 
efforts in preparation and time available for TEA. It seems also likely 
that the experience gained from EOVV carried over to TEA and sustained 
our confidence. The demonstrated use of aerodynamic torques for space- 
craft attitude control seems destined to be of value for design of future, 
low altitude spacecraft or tethered vehicles. 
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APPENDIX A 

AERODYNAMIC TORQUE MODEL 

A set of six coefficients were developed to express forces and 
torques in body coordinates as a function of the velocity of the vehicle 
relative to the air. If this velocity is V, then the forces and torques 
relative to the origin of coordinates are 

F = 0.5 6 V 2 A C p (Ala) 

T = 0.5 6 V 2 A D C T (Alb) 

3 

where 6 t Kg/m ] is the density and 

A = 79.46 m 2 (855.3 ft 2 ) reference area 

D = 10.06 m (396 in.) reference diameter 

Cp = [-CA, CY, -CN] T force coefficients 

T 

C,^ = [CL, CM, CENJ moment (torque coefficients). 


The force and moment coefficients had been calculated by Northrop 
[6] with a computer program generated by Lockheed [7] and are given 
in Table Al. These coefficients depend only on the direction of V (it is 
to be noted that the density is assumed to be small enough that free 
molecular flow assumptions apply). This direction is customarily expressed 
as aerodynamic roll angle and angle of attack, p and u respectively. 

These are defined in the Skylab vehicle as (cf. Eqs. 3 and 4): 


a = arc cos ( V x / j V ] ) 

p = urc tun (V /V ) (lour quadrant) 
y z 


(A2) 

(A3) 
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TABLE Al. FORCE AND MOMENT COEFFICIENTS 
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7 626 

2.00? 

d.237 

2.2*2 

ti 232 

2 275 

6.757 

2 .8 36 

4.18! 

1.095 

2 552 

392 

00 U 

000 

CM 

CL 

.002 

00 1. 

960 

.270 

2.366 

750 

4 279 

1 .590 

5.647 

2.0 67 

5 874 

2.262 

5.659 

2 .095 

4.746 

2.721 

3.1*73 

1 «2H 

2 .243 

357 

.000 

. O' O 0 

CM 

cr, 

. 00? 

0 00 

.609 

199 

2 419 

519 

2 415 

.991 

3 -026 

1.296 

3 12 5 

1 . 350 

3.04 0 

1 .319 

2.57R 

1 .096 

1.695 

.688 

.716 

.265 

.^00 

. 000 


CM 

CEN 

- 332 

-.298 

1 578 

-.422 

•3 . 62ti 

-.719 

■5.041 

-.522 

6 132 

-.244 

6 423 

-.001 

5.767 

.291 

4 .302 

776 

■2 .361 

.943 

- 4ti7 

.732 

304 

.436 


CM 

CFN 

- 332 

— 298 

-3 .492 

-.143 

-3 294 

.12 5 

-4.677 

879 

-5 .6 39 

1.625 

-5. 80 7 

2.041 

-S.253 

2 253 

-3.842 

2 42 8 

-2 . 885 

1 .944 

-.419 

2 .112 

304 

.436 

CM 

CUV 

- 332 

-.298 

-2 .249 

.221 

-2.490 

.92 0 

—3.778 

1.971 

-4 .170 

2 .925 

-4 . 04 5 

3.273 

-3.564 

3.345 

-2.578 

3.302 

-1.326 

2 . 537 

- 399 

1.261 

304 

436 

CM 

CKN 

- 332 

- . 298 

-.937 

.527 

-1 .632 

3 .583 

-2.263 

2 .744 

-2.280 

3.434 

-2 124 

3.606 

-1 . 858 

3.599 

-1 .353 

3.459 

-.693 

2 638 

— .0ft9 

1.356 

• 304 

.436 
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RHO * 80. ALPHA 

0 . 
28. 
48. 
68 . 
88 . 
90 . 
100 . 
120. 
140 . 
166. 
180. 

RHO * 80. ALPHA 

8 . 
20. 
40. 
60 . 
80 . 
80 . 
100 . 
126 
140. 

• 160. 
188 

RHO * 100. ALPHA 
0 . 
20. 
40. 
60 . 
80 . 
98 . 
1 0 0 
120 . 
140. 
160 . 
180 . 

RHO *128. ALPHA 
0 . 
20. 
40 . 
6U . 
80 . 
90 . 
100. 
120. 
140. 
160 . 
180 . 


TABLE Al. 


CA 

CY 

2.401 

-.002 

2.600 

- 994 

2.669 

-2.295 

2.110 

-3 711 

.812 

-4.615 

.004 

-4.715 

-.801 

-4.600 

-2.111 

-3.703 

-2.760 

-2.348 

-2.700 

-1 .011 

-2 590 

-,001 

CA 

CY 

2.401 

-.002 

2.707 

-1.048 

2.687 

-2.346 

1.998 

-3.551 

.741 

-4 . 280 

.004 

-4 348 

-.730 

-4.269 

-1 988 

-3.545 

-2.742 

-2.372 

-2.813 

-1.066 

-2.590 

- 001 

CA 

CY 

2.401 

-.002 

2.517 

- 966 

2.710 

-2.330 

2.154 

-3 778 

.832 

-4.721 

.004 

-4 857 

-.827 

-4 .734 

-2.186 

-3.825 

-2.810 

-2.392 

-2 . 630 

- 988 

-2.590 

-.001 

CA 

CY 

2.401 

-.002 

2.952 

-.985 

3 522 

-2.637 

2.854 

-4.367 

1.122 

-5.576 

.004 

-5.826 

-1.128 

-5.637 

-2.943 

-4 494 

-3.513 

-2.613 

-3.044 

-.998 

-2.590 

- 001 


(Continued) 


CN 

CL 

.002 

.000 

.227 

.126 

.468 

.256 

.733 

.475 

.903 

.640 

.922 

.664 

.900 

.639 

.734 

.474 

.491 

.265 

.248 

.141 

.000 

.000 

CN 

CL 

.002 

.000 

-.001 

.117 

- 001 

.227 

— .001 

.332 

.000 

.386 

. 000 

.391 

.000 

.385 

.000 

.333 

.000 

.223 

-.001 

.112 

.000 

.000 

CN 

CL 

.002 

.000 

-.224 

.110 

- 472 

.248 

-.751 

.433 

-.933 

.567 

-.961 

.575 

-.938 

.552 

-.768 

.395 

-.500 

.199 

-.246 

.084 

.000 

.000 

CN 

CL 

.002 

.000 

-.663 

.133 

1.655 

• 397 

2.691 

.654 

3.422 

. 722 

3.578 

.668 

3.465 

.626 

2.778 

.445 

1.643 

.204 

- 678 

.063 

.000 

.000 


CM 

CEN 

-.332 

-.298 

-.497 

.728 

-.695 

2.049 

-.874 

3.243 

-.829 

3.899 

-.730 

3.984 

-.593 

3.862 

-.329 

3.255 

-.155 

2.373 

.091 

1.225 

.304 

.436 

CM 

CEN 

-.332 

-.298 

- 317 

.706 

-.269 

2.062 

-.192 

, 3.221 

-.068 

3.865 

.000 

3 948 

.068 

3.844 

.193 

3.257 

.265 

2.378 

.305 

1.350 

.304 

.436 

CM 

CEN 

-.332 

-.298 

-.150 

.788 

.086 

2.061 

.355 

' 3.280 

.652 

4.032 

. 785 

4.181 

. 872 

4.104 

.,904 

3.503 

.735 

2.422 

.518 

1 .125 

. 304 

.436 

CM 

CEN 

-.332 

299 

-.030 

.663 

.606 

1.926 

1.525 

3.279 

2.423 

4.279 

2.755 

4.572 

2.875 

4.550 

2.722 

3.977 

1 .923 

2.634 

1.035 

1.218 

.304 

.436 
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TABLE Al. (Continued) 


RHO * 14k:. 


RHO = 1M. 


RHO = 18b. 


RHO * 


ALPHA 

CA 

C¥ 

CM 

CL 

CM 

ct;N 

3 . 

2.401 

-.00 2 

. 002 

. 001) 

- 332 

- . 293 

22 . 

3 . 340 

-.821 

-1.095 

. 095 

.0 75 

.413 

40 . 

4 .202 

-2.322 

-2.953 

. 303 

3 .117 

1 . 379 

60 . 

3.423 

-3.876 

-4.853 

.352 

2 . 80 4 

2.552 

80 . 

1.315 

-4 . 84 0 

-6 . 029 

. 363 

4 .102 

3.42 7 

90 . 

. 00 4 

-5.072 

-6 . 320 

. 264 

4 . 593 

3.735 

100 . 

-1 . 32 7 

-4 871 

—6 . 069 

.223 

4.740 

3.799 

120. 

-3.2yd 

-3 .731 

-4 661 

.129 

4 32 3 

3 388 

14 0. 

-3.981 

-2 . 194 

-2.777 

. 04 8 

3.011 

2 . 316 

J 60 . 

-3 .209 

-.783 

-1.041 

-.02 3 

1 . 376 

1.096 

]Hi> . 

-2 . 590 

* t < o 1 

. 000 

, 000 

. 304 

.436 

ALPHA 

CA 

CY 

Civ 

CL 

CM 

CEN 


2.401 

- . 0 U 2 

00 2 

. 00 0 

-.332 

-.299 

2 * . 

3.63 0 

-.473 

-1 .434 

02 1 

.169 

. 038 

4 0 

4 . 713 

-1 . 3H4 

-4 .022 

-.119 

1.598 

.449 

to 

i . 657 

-2 . 320 

-6.658 

865 

3.719 

3 .121 

00 . 

1.471 

-2 .876 

-8.215 

- 592 

5.471 

1 .839 

90 . 

. 00 5 

-2.923 

-9 . 339 

-.693 

5.990 

2.121 

1 (i 0 . 

-1 399 

-2.754 

-7 .853 

-.549 

5 . 998 

2.216 

120 . 

-3.566 

-2.148 

-6.146 

-.453 

5 . 358 

2.146 

14 0. 

-4 . 308 

-1 .264 

-3.668 

-. 324 

3. 712 

1 .662 

J 60 . 

-3 . 31 3 

- 429 

-1 .298 

-.122 

1.648 

. 846 

180 . 

—2 590 

- 001 

. 0 0 0 

.09 0 

. 304 

.436 

ALPHA 

CA 

CY 

Co 

CL 

CM 

CLN 

: .l , 

2 ■ 401 

- 0 02 

.002 

.000 

-. 332 

-.298 

20 . 

3 86 0 

- 0 o 1 

-1.618 

-.231 

. 335 

-.536 

All . 

4.873 

- oftl 

-4 .416 

-.632 

1 .794 

-.686 

60 . 

4. 01 7 

- 001 

-7 368 

-1 . 252 

4.023 

-.674 

8u . 

1.51 8 

_ i'.i.M 

• W »/ i 

-8 . 995 

-1 .573 

5-782 

-.264 

90 . 

. ii 0 5 

- 001 

-9.236 

-1.671 

6. 37 9 

. 0 00 

1 0 0 . 

-1 .406 

- 0 '0 1 

-8 . 6 Ali 

-1 . 489 

6.346 

. 248 

120 . 

-3.515 

- 001 

-6.438 

-1 . 02 8 

5 334 

.556 

140. 

-4 . 306 

- 0 O 1 

-3 . 886 

-.70 6 

3 . 821 

.767 

160 . 

-3. 344 

-00 2 

-1 . 3 9 •. ) 

-.225 

1.731 

. 519 

180 . 

-2.590 

- 001 

. t’t 0 0 

. 000 

. 304 

.436 

ALPHA 

CA 

CY 

CM 

CL 

CM 

CEN 

o • 

2 . 4 0 1 

-.00 2 

. 002 

.000 

-.332 

-.298 

20 . 

3.554 

.4 60 

-1 . 404 

-.2 37 

.116 

-.665 

40 . 

4.672 

1 369 

-3.989 

-.885 

1.569 

-1 . 590 

60 

3 . 84 0 

2.307 

—6 .62 9 

-1 .704 

3.678 

-2.284 

80 . 

1.436 

2 80 7 

-9 . 015 

-2.171 

5 . 210 

-2 .234 

9o . 

. 00 5 

2.875 

-8.203 

-2.230 

5.757 

-2.04 3 

100. 

-1 . 359 

2.674 

-7.625 

-1 .969 

5.740 

-1.687 

120. 

-3 . 321 

1.999 

-5.715 

-1 . 352 

4.938 

- 970 

140. 

-3.882 

1 136 

-3.295 

- 82 6 

3.383 

-.2 30 

160 . 

-3 246 

.418 

-1 .271 

-.280 

1.630 

.136 

1 80 . 

-2 . 590 

-001 

. 000 

. 000 

. 304 

.436 
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TABLE Al. 


^-•‘''Puv.iDUC. 

(Continued) ORIGINAL ; 


K’L'.HY OF THB 

i'Auil! 13 POOR 


RHO * 220 . 


RHO * 240. 


RHO = 260. 


RHO .* 270 . 


ALPHA 

CA 

CY 

CN 

CL 

CM 

CEN 

0. 

2.401 

-.002 

.002 

.000 

-. 332 

-.298 

20 . 

3.212 

.779 

-1 .058 

- 260 

.009 

-.824 

40. 

3.909 

2 . 137 

-2.759 

— . 832 

.832 

-1.830 

60 

3.187 

3.585 

-4 532 

-1.591 

2.383 

-3.015 

80 . 

1.224 

4.482 

-5.615 

-2.097 

3.532 

-3.302 

90 . 

.004 

4.714 

-5.902 

-2.244 

4.025 

-3.252 

100 . 

-1.225 

4 510 

-5.655 

-2.148 

4.192 

-2.918 

120 

-2.767 

3.114 

-3.928 

-1.397 

3.662 

-1.969 

140. 

-3.255 

1.778 

-2.281 

-.771 

2.585 

-.963 

160 

-3.029 

.733 

-.982 

-.302 

1. 318 

-.138 

180 . 

-2.590 

-.001 

.000 

. 000 

. 304 

.436 

*♦ 

ALPHA 

CA 

CY 

CN 

CL 

CM 

t 

CEN 

0 

2.401 

-.002 

.002 

.000 

- 332 

-.298 

20 . 

2.798 

.920 

-.632 

-.222 

-.086 

-.908 

40. 

3.063 

2.264 

-1.456 

-.580 

.346 

-1.769 

60 . 

2.473 

3.748 

-2.357 

-1.014 

1.082 

-2 788 

80 . 

.976 

4.809 

-2.993 

-1.344 

1.828 

-3.403 

90 . 

.005 

5.065 

-3.149 

-1.423 

2.166 

-3.531 

100. 

-.974 

4.832 

-3.020 

-1.373 

2.317 

-3.321 

120. 

-2.419 

3.659 

-2.312 

-1.022 

2.292 

-2.65 7 

140. 

-2.962 

2.179 

-1.406 

-.587 

1.727 

-1.575 

160 . 

-2.901 

.940 

-.649 

-.246 

1.013 

-.501 

180 

-2.590 

-.001 

.000 

.000 

.304 

.436 

ALPHA 

CA 

CY 

CN 

CL 

CM 

CEN 

0. 

2.401 

-.002 

.002 

.000 

-.332 

-.298 

20. 

2.370 

.898 

-.212 

-.126 

-.175 

-.836 

40 . 

2.480 

2.109 

-.440 

-.281 

037 

-1.836 

60 . 

1.980 

3.445 

-.702 

— . 500 

.277 

-2.874 

80 . 

.786 

4.423 

-.893 

-.667 

.594 

-3.664 

90 . 

.005 

4 .626 

-.935 

-.701 

.747 

-3.858 

100 . 

-.779 

4.431 

-.900 

-.671 

.833 

-3.698 

120. 

-1.999 

3.464 

-.717 

-.507 

.884 

-3.070 

140 . 

-2.595 

2.186 

-.468 

-.290 

.732 

-2.037 

160 . 

-2.556 

.951 

-.239 

-.343 

.518 

-.745 

180 . 

-2 . 590 

-.001 

.000 

.000 

.304* 

.436 

ALPHA 

CA 

CY 

CN 

CL 

CM ■’ 

CEN 

0 . 

2.401 

-.002 

.002 

. 000 

-. 332 

— . 298 

20 . 

2.625 

1.005 

.006 

-.121 

-.339 

-1.134 

40. 

2.628 

2.269 

.010 

-.234 

-.292 

-2.149 

60 . 

1.988 

3.516 

.009 

-.333 

-.203 

-3 . 251 

80 . 

.763 

4.360 

.004 

-.370 

-.069 

-4.011 

90 . 

.005 

4.515 

.001 

-.362 

- . 001 

-4 . 164 

100 . 

-.754 

4 359 

.003 

-.367 

.061 

-3.969 

120 . 

-2.049 

3.615 

.007 

-.317 

.177 

-3.281 

140 . 

-2.829 

2.424 

.006 

-.212 

.246 

-2.117 

160 . 

-2.858 

1.074 

.002 

-.110 

.297 

-.659 

180 . 

-2.590 

- 001 

.000 

.000 

.304 

.436 


TABLE Al. (Concluded) 


RHO * 28b . ALPHA 

0 

2b . 
40 . 
bki 
80 
y u . 
l«b 
120. 
14b . 

no . 

1 bb 

RHO - 3 ;.v). ALPHA 

o « 
2 0 . 
4ti . 
60 . 
80 . 
9 b . 
1 00 . 
120. 
14b . 
lb.:) . 

18b . 

RHO = 320. ALPHA 

V . 1 . 

2b . 
40 . 
60 . 
hb . 
90 . 
lOO . 
12b .. 
14b . 
160 - 
1 HO . 

RHO * 340. ALPHA 
Vi . 
20 . 
49. 
60 . 
80 . 
9b . 
109 . 
120 . 
1 4 9 . 
160 . 
18b . 


LA 

LY 

2 . 4b] 

-.002 

2.576 

.977 

2.706 

2 . 307 

2.164 

3.775 

. 84 7 

4.776 

. v) 0 5 

4.945 

-.839 

4.780 

-2 .204 

3.834 

-2 . 86 9 

2.421 

-2 . 760 

1.027 

-2 . 599 

— b 0 1 

LA 

CY 

2 . 4 0 1 

- . b v* 2 

2 8 a 2 

.94 3 

3.266 

2 457 

2.828 

4 . 321 

.i.125 

5.8*8 

. O 0 4 

5 . 8b 6 

] v , 

.4 *4 3 > 

5- 628 

-2 .987 

4.548 

-3 769 

2.784 

-3 . 299 

1.070 

-2 59 i 

-.001 

C.A 

L i 

2 .4 01 

-.002 

2 94 5 

72b 

3 . 699 

2 . 033 

3.18 7 

3.589 

1 .290 

4 . 74 O 

- . 0 0 1 

4 . 99 0 

-1 . 323 

4 . H4 6 

-3 .535 

3.973 

-4 533 

2.476 

-3 . 646 

. 874 

-2 . 59b 

- 001 

LA 

LY 

2 . 4 l’ 1 

- 002 

3 • i/99 

. 4 05 

4 .145 

1 . 209 

3. 407 

2.137 

1 . 371 

2 . 686 

-.003 

2 93 3 

-1 . 492 

2.911 

-3.971 

2.374 

-4 . 95 6 

1.442 

-4 . i: u 7 

. 510 

-2 .590 

- 001 


CN 

CL 

. 0O2 

. 000 

.235 

-.106 

. 493 

-.226 

780 

-.428 

.976 

-.553 

1 . 009 

-.556 

. 991 

-.542 

. HOP, 

-.378 

.526 

-.166 

.257 

-.0 70 

0 0 V) 

. 0 0 O 


CN 

Cl. 

. 002 

. 0 M 0 

. 624 

— .106 

1 .548 

-.292 

2 678 

- 536 

3.436 

-.65 9 

3.584 

-.64 0 

3.478 

-.60 6 

2 . H 3 8 

-.420 

1.785 

-.282 

.739 

— . i)4 9 

. b0 ii 

. 000 

CN 

CL 

. 002 

. 0 1 0 

.951 

— .045 

2.62 5 

-.161 

4 . 558 

-.231 

5.971 

-.295 

6.291 

I«6 

6 . 114 

-.167 

5.037 

-.073 

3.189 

. 01 3 

1 .189 

. 041 

. V) vi O' 

. 000 

O' 

CL 

.002 

. 0 W 0 

1.208 

.054 

3.566 

. 268 

5.923 

. 383 

7.762 

. 493 

8.460 

.70 9 

8.393 

. 755 

6.867 

672 

4 . 224 

.433 

1.566 

178 

. 000 

. fc)00 


CN 

CKN 

-.332 

-.298 

-.516 

-1 .0 50 

-.729 

-2 . 227 

-.935 

-3 410 

-.925 

-4 .168 

-.846 

-4 . 286 

-.70 8 

-4.078 

-.426 

-3.291 

-.221 

-2 . OHO 

.0 76 

-.690 

. 304 

.436 

CM 

CON 

-.332 

-.298 

-.951 

-1 . 1)0 ■ 

-1 . 774 

-2.271 

-2.577 

-3.671 

-2.845 

-4 .481 

-2 772 

-4.579 

-2 . 469 

—4 . 287 

-1.72] 

-3. 314 

- . 824 

-] . 846 

- 122 

-.519 

. 30 4 

.436 

CM 

CPN 1 

-.332 

-.298 

-1 .235 

-.827 

-2.778 

-1 .954 

-4.219 

-3.120 

-4 62 7 

-3.656 

—4 629 

-3.723 

-4 . 171 

-3 . 393 

-2.977 

-2.472 

-2 . 462 

-1.191 

-. 26b 

-.137 

. 304 

.436 

CM 

COM 

-.332 

-. 2 t >8 

-1 .468 

-.651 

-3.533 

-1 . 522 

-5 . 03 1 

-1 .972 

-5.83 0 

-2.132 

—6 . 0 5 1 

-2.125 

-5.535 

-1.814 

-4 . 008 

-1.034 

-1.933 

-.192 

-.439 

. 270 

. 304 

.4 36 
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KEPRODUCiniUTif OV THE 
ORIGINAL PAGE IS POOR 

The aero data provided us consisted of six coefficients C^, and C T for 

each of 231 directions defined by p = 0,20, ..., 360,90,270 and a = 0, 20, 

. . . 18° To get forces and torques for other directions we had to 

interpolate and initially we used a linear interpolation scheme. This was 

satisfactory at first but later caused problems in the calculation of TEA 

control coefficients. Ultimately a Fourier fit of the data was used. This 

model was set up such that each table value was reproduced and since 

the in-between values were now produced by a sinusoidal series, the 

variations were extremely smooth and computation problems were 

eliminated. The coefficients C. (a. , p ) can be expressed as a finite 

l j k 

series of sines and cosines in a two-step process. 

10 10 

C-(ou,p. ) = V a. (a.)cos mp. + Y b. (a.)sin mp, . (A4) 

i ] u im ] k Li im j k 

m=0 m^l 


It can be observed that the given set of angles for p is equivalent to the 
following more convenient set: p - 0,±20, . . . ,±180, ±90. We can simplify 

our task by defining two auxiliary functions G and H such that 


G(a,p) = 0.5 [C(a,p) + C(a,-p)] 

H(a,p) = 0.5 [C(a,p) - C(a,-p)] 

These functions can be expressed as 

10 

G.(a. ,p, ) = V a. (a.) cos mp. 

l j K k Li im ] K k 

m=0 

10 

H i ( a.,c k ) = l b. m (<y sin mp k . 
m=l 


( A5) 
(A6) 


(A7) 


( A8) 


The functions G.(a,p. ) and H.(a. , p, ) can be computed for each a. and 
l 3 k l ] k j 

equations (A7) and (A8) inverted to obtain a. and b. . In turn these 
^ im im 

coefficients can be expressed as a fourier series in a. Since the given 


51 


set of a's only range from 0 to 180, we can make the assumption that 

a. and b. are symmetric in a. This means we can express a. and 

lm lm J * im 

b. as cosine series: 
lm 

10 

W = <A9) 

n=0 

10 

b in ,^> = Z b im cos no ^ ♦ < A10 > 

n=l 


Using equations (A7) through (A10) the coefficients a. and b. 

im im 

can be determined. Inspection of the above will reveal that the solution 
procedure only requires inversion of the 11 x n matrix (cos ma), where 

m = 0,1,. ...10 and a. = 0,20 180,90; and the 10 x 10 matrix (sin 

idol), where m and ou are as before except the 0 is removed. The final 

form of the resulting expansion is 


10 / 

<10 

10 - \ 



/. cos na( 

Y. a. „ cos mp 

+ y. b. sin mp | 


(All) 

u \ 

imn 

u imn 1 



n=0 ' 

vm=0 

m=l / 




The subscript i ranges over the values xyz for torques and xyz for 
forces. Thus there are six components of C. and for each component 

there are 121 a. and 110 b. for a total of 231 per axis which makes 
lmn lmn r 

1386 total Fourier coefficients. 

The uncertainty in the aero torques are represented by a center 
of pressure (CP) offset. The nominal set of coefficients are computed 
for the given coordinate axes and origin. Torques are computed about 
the origin and the total force is assumed to act through it. The torque 
about the center of mass is 



T 


— o 


- R 


—cm 


X F 


-o 


( A12) 
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in the nominal case. If the center of pressure is moved by A, the 
torque becomes 


T = T + A x F 
— ae — n - -o 


(A13) 


The vector T is the torque vector required for the aero model we 
used. The nominal center of mass location, R^, , was taken from Refer- 
ence 8. Details are given in Tables A2 and A3. The units are m(in.). 


TABLE A2. NOMINAL CENTER OF MASS LOCATION 



Moment Reference Point, 
m (in . ) 

Nominal CM Location, 
m (in.) 

CM, m (in.) 

X 

90.037 

81.89210.155 

-8.14510.155 


(3544.765) 

(3224.116.1) 

(-320.66516.1) 

Y 

0 

-0.10710.053 

(-4.212.1) 

-0.10710.053 

(-4.212.1) 

Z 

0 

+0.83810.064 

(33.212.5) 

+0.83810.064 

(33.212.5) 


TABLE A3. CENTER OF PRESSURE LOCATION 



ACP . , m(in . ) 
est 

ACP acf m(in ' ) 

ACM m(in.) 

ett 

X 

0+1.524(0160) 

0.06(+2. 36) 

-8.205(-323. 03) 

Y 

0+0.762(0130) 

-0. 07(-2. 76) 

-0.037(-1.46) 

Z 

0+0.762(0130) 

-0.1K-4.33) 

+0. 948( +37 . 32) 


As can be seen in Table A3 the tolerance on the CP location is an 
order of magnitude larger than the uncertainty in the center of mass 
location. Both have the same effect on the predicted aerodynamic 
torques and as a consequence, the CM tolerance can be ignored. Before 
actual flight data were available, the large CP tolerance was used to 
assess the sensitivity of the TEA seeking method (see TEA Seeking 
Method) . After flight data were available , the actual estimate given in 
column 2 of Table A3 was used to give the effective ACM given in 
column three. 
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Zero contours for the three torque components are shown in Figure 
4. Asterisks indicate all TEA'S (the feasible ones are also labeled). 
Individual contour maps for the torque components are given in Figures 
Al, A2, and A3. The contour numbers signify actual aerodynamic 
torques in Nm for a density of 3.11E-10 kg/m*, corresponding roughly 
to an altitude of 200 km (108 n.mi.) and nominal solar activity. 7 Table 
A 4 shows contour number multiplication factors for other altitudes (con- 
tour number times multiplication factor equals torque). 

Stereo doubles (see Appendix F for an explanation) of Figures Al, 
A2, and A3 are given in Figures A4, A5, and A6. Stereo projections 
of the contours on unit spheres are shown in Figures A7 through A 12. 

TABLE A4. CONTOUR NUMBER MULTIPLICATION FACTORS 


(km) 

Altitude 

(n.mi. ) 

3 

Density (kg/m ) 

Multiplication Factor 

278 


5.310E-11 

0.17 

259 


7.855E-11 

0.25 

241 

130 

1. 197E-10 

0.38 

222 

120 

1.885E-10 

0.61 

204 

110 

3. 110E-10 

1.00 

185 

100 

5. 485E-10 

1.80 

167 

90 

1. 065E-09 

3.40 

148 

80 

2 . 505E-09 

8.10 

130 

70 

8. 920E-09 

28.70 


7. Nominal data for 6/15/79 with sun spot activity number 136.2, 

F10 = 181.66 and A = 17.96. 

P 
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Figure A4. X torque contours (stereo double). 
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APPENDIX B 

GRAVITATIONAL TORQUE MODEL 


The gravitational torques on a satellite produced by a large, 
spherical primary body are important contributors to its rotational 
dynamics. The force on a point mass m exerted by the primary K is 


F _ _ _GMm p . G = 6 . 6 72E-11 Nm 2 /kg 2 . (Bl) 

■* l**i I 3 ^ 

The satellite can be viewed as a collection of point masses. The net 
torque on this collection of masses about the origin of satellite 
coordinates is 



*2»i *h 


(B2) 


The vector r. is the position of m. relative to the origin which is at R 
relative to the primary center. Thus 




R + r. 
— o -i 


and 




x GMm. 
i 



(B3) 


(B4) 


In general |r.| << Ij^l and hence an expansion of keeping only low 
order terms becomes sufficient for most purposes. Now 






+ h.o.t. 


( B5) 
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Using the expansion of equation (B5) in (B4), 


T 

-g 


r 

-cm 


GMm 





E 


m i ^i 



+ h.o.t. 
(B6) 


We have used the definition mr = Ym.r. . Rearranging and grouping 
equation (B6) we obtain ~ cm . 1 


T 

S 


- r 
-cm 


x 


GMm 







+ h.o.t. . 


(B7) 


The term in parentheses in equation (B7) occurs in the definition of the 
moment of inertia dyadic (or tensor) 

I = ^ m i (r. 2 1 - r. r.) . (B8) 

i 


Using this definition in equation (B7) and dropping the higher order 
terms yields the gravity gradient torque expression 


T 

-gg 


- r 
-cm 


GMm „ 

I n I 3 ^ 


+ 3 


“L-J-.i 

I5»l l R o' 



(B9) 


As can be seen, if the origin is positioned at the center of mass the 
more familiar gravity gradient torque expression results: 


T 

~gg 




(BIO) 


For an orbiting body m , the orbital angular velocity magnitude is given 
by 
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Noting that Rq/IRqI is a unit vector u^, we finally write 


= 3 oj 


HR * 1 * Hr 


(BU) 


(B12) 


The torque of equation (B12) is what is commonly referred to as the 
gravity gradient torque. 
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APPENDIX C 


DETERMINATION OF TORQUE EQUILIBRIUM ATTITUDES 

The torque equation for a rotating vehicle in vehicle fixed coordi- 
nates is 


H - <*> x h 


T 

-g* 


+ T 
— «e 


(Cl) 


• 

The equilibria we seek are those which satisfy H = 0. In the case of 
Skylab , the momentum H is not only from vehicle rotation at orbital 
velocity but also includes vehicle fixed CMG momentum. Since the dis- 
turbance torques on the right hand side of equation (Cl) are smooth 
functions, we can use Newton's method to determine equilibrium positions. 
The technique starts from the assumption that 


H = F( [VL] ) 


(C2) 


i.e., the rate of change of angular momentum in body axes depends upon 
the relative attitude only. The symbol [VL) is the direction cosine 
matrix or transformation which takes a vector from the L (local vertical) 
coordinate system to the V (vehicle fixed) system: 


F » F ([VL]) + A0 • VF ; I A 6 small . (C3) 

— *~o — — o 


The basic algorithm then is 

1) Select [VLJp 

2) F = F([VL] ) 

”P P 

3) A 6 = Urn l-(VF ) -1 F ] (C4) 

— ~P “P 

4) [VLJ n = [A6 z ] 3 [A0 y ] 2 [A0 x ] 1 [VL] p 

5) If I A 0 1 < e, set [VL] = [VL] and terminate. 

1 — n 

6) Otherwise, set [VL] p = [VL] n and loop to step 2). 


PRECEDING 


PAGE Bl-ANK 


l^OT FILMED 
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The limit in step 3) is a component limit which limits the changes 
of attitude permitted in the iteration procedure above and improves con- 
vergence. With this procedure 12 equilibrium attitudes were found. 
Others may exist and an exhaustive search was not conducted. These 
equilibria could be shifted somewhat by varying the momentum stored 
in the CMG's. The total momentum H was held constant in the L frame 
to H nom * The gradient of the vector function F was computed 

numerically. 
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APPENDIX D 


STRAPDOWN ERROR ESTIMATION 


The strapdown correction angle A6^ was estimated from Skylab 

attitude history over one or more orbits in TEA. The assumption is that 
the strapdown error changes slowly over an orbit, especially the com- 
ponents of the error in the orbital plane. This assumption is justified 
by the observation that , since the vehicle orientation stays nearly fixed 
relative to the local vertical, the gyro drift resolved into the orbital 
plane is nearly sinusoidal and hence averages to near zero over an 
orbital period. The out-of-plane drift remains essentially constant and 
hence out-of-plane strapdown error grows linearly with time. The three 
angles 0 el , 0 e2 , and 0 fi3 represent the strapdown offset error. The 

error matrix is 


i°° k J 


= [6 el ] l 


te e2 J 2 


* 6 e3 ] 3 


(Dl) 


where a k subscript signifies the true quantity and O is the estimate of 
the real orbital coordinate system O^. The coordinate system O is 

defined such that O coincides with L when the vehicle crosses the ascend- 
ing node and is inertially fixed except for orbit regression. We shall 
assume orbit regression is too small to have a significant effect on strap- 
down error determination. We further assume 0 . and 0 „ are small 
angles. The position of the vehicle is then e e 


[V( V 


= [VL] 




[ % 2 ] 2 


{ W 3 


(D2) 


Using the small angle assum^t’^n , we can write 


[VO k ] = [VLJ [1 - O'] l* 0 t * 0 e2 J 2 


where 


0 ' 

-e 


Iw o tJ 2 t6 el’ °* 0 e3 ] 


<D3) 


(D4) 
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We can also write 


[VO k ] = [VL k ] [ Wo t] 2 . (D5) 

In this expression (VL^) is assumed to remain constant. We can now 
equate equations (D3) and (D5) to obtain 

[VL] = !VL k J t-e e2 i 2 11 - e e M . CD6) 


We must also assume that the gyro drift and other drift rates are small 
so @ e2 does not change significantly over an orbit. With this assumption 

[VL.][1-0'1 remains essentially constant. Since the location of the TEA 

is not known precisely, [VL k ] is unknown and of course at this point 


8 e2 is also unknown so that we may as well 




[VL k -] = [VL k ] l-e e2 J 2 


(D7) 


Thus all that is known is that (VL k ] is constant, that 0^ varies sinu- 
soidally, and that [VLl is available over one or more orbits from teleme- 
try. This matrix is used to compute BAR angles. Intuitively it would 
seem that the bank angle would contain the most direct information on 
the in plane strapdown error, so let us look first at it. 


tan 6 = 


(V W- VL k'l3> - e e'l * (VL k'll /VL k'l3 ) 9 e3 
1 + (VL k12 /VL k13 ) 6 e'l 


(D8) 


For high drag TEA, is small (the angle of attack is near 90 

degrees) . Thus the inverse tangent of equation (D8) yields 

6 * 6 k - e el = S k ' 6 el co8( “o t) + 8 e3 sln( “o ,) (D9) 

= 3 + cos(co it) + £ c sin(o> t) . (DIO) 

o C o b o 
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: t. '*ty- ■ ■.***-. ■sm** --'*. •* - 


We can extract the strapdown in-plane errors from equation (DIO) and 
construct the update vector: 

2.' * ‘ e eV °- 9 e3 ]T * [ 'V °- 6 s ,T • <»'» 


The strapdown correction A 6 L is developed from o' as 





O’ 

-e 


(D12) 


The time 1 is the time of the update relative to the time of duta. 


The out-of-plane strapdown error could only be determined from 
timing and Sun passage data. The basic procedure for thi.s was the 
same whether Sun sensor or solar array data were used. In either cose 
it was assumed the strapdown error was along the orbit normal, i.e. , 

0 i = 0 e3 = 0. The timing events used occurred when the Sun made a 

known angle from a known body- fixed direction. This can be expressed 
mathematically as 


U • U 
-s -v 


= C. 


sra 


(D13) 


This was most easily solved in the L coordinate system: 


V = (LI) tO, 0, 1] 

“D 


(D14a) 


U v - (VL] T U f 


(D 14b) 


where Up is a vehicle-fixed unit vector (see below, after equation D17c) 

The matrices [LI] and [ V L] are computed at the timing event instant. 

At this instant we can determine 0 0 : 

-e2 


'sra = < U sx* V U « J ^v 


(D15) 
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From equation (D15) we obtain 


e «2 


= arc 


sin [(C 


sra 


- U 


sy 




2 9 

1 + C<T 


]- 


a 


or 


e2 


= ir - arc 


"‘"Osra - U »y V //s ^ + C *^‘ * 


(D16) 


with 


S <p 


u 


sx 


II + U 
vx sz 


U 

vz 


(D17a) 


C <p 


V u - u u 
sx vz sz vx 


(D17b) 


4) = arc tan (S4>/C<J>) . (D17c) 

T 

The Sun sensor null crossing occurred when U_ = (0, 1, 0] and 

— r 

C = 0. The power angle was determined from sunrise on the solar 
sra fp 

arrays so that = (0, 0, lj and for various values of C gra depending 
on which array was chosen. 
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APPENDIX E 

QUATERNIONS - A BRIEF EXPOSITION 


Complex numbers of the form z = a + ib have proven to be a 
valuable concept in the study of many physical phenomena. A generaliza- 
tion of this concept which proves useful in the study of rotational motion 
is the quaternion. Recall that the imaginary unit i = /-l. Let us 
define additional units 3 and k together with the product operation o: 


i 0 i 

= i 0 j = 

k 

0 

V? 

11 

i 0 j 

= - j 0 i 

= 

k 

3 0 k 

= - k 0 j 

= 

i 

k 0 i 

11 

1 

>-> 

c. 

= 

j 


(El) 


We shall define a quaternion as any quantity of the form 

Q = Q 4 + i Q1 + j Q 2 + k Q 3 . (E 2 ) 


(See Note 1 at end of Appendix E.) By analogy to the complex number 
terminology Q 4 11 referred to as the real or scalar part of 0 . It will 
also be convenient to think of the remaining part of Q as the imaginary 
or vector part. The reason for this will become clear as we proceed. 
Let R = R 4 + iRl + jR 2 + kR 3 . The sum of quaternions Q and R is 
defined as 


s = Q + R = (Q4 + R4) + i ( Q 1 + Rl) 

+ j (Q2 + R2) + k (Q3 + R3) . (E3) 


From this definition we can see that the sum operation is commutative 
and associative. We can now give the complete definition of the jiroduct 
operation o: 
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P = Q o R = (Q4 R4 - Q1 R1 - Q2 R2 - Q3 R3) 

+ i (Q4 R1 + Q1 R4 + Q2 R3 - Q3 R2) 

+ j (Q4 R2 + Q2 R4 + Q3 R1 - Q1 R3) 

+ k (Q4 R3 + Q3 R4 + Q1 R2 - Q2 Rl) . (E4) 


With this definition we can show that o is associative and distributive but 
not commutative, i.e.,QoR#RoQ. We shall call any quaternion 
having zero imaginary part a scalar and, obviously, the algebra of 
scalars is just the algebra of real numbers. Thus, multiplication of a 
quaternion by a scalar simply results in a quaternion whose elements are 
multiplied by that scalar according to definition E4. We can now define 
the difference operation as 


D = Q - K = Q + (-1) o R 


(E5) 


For convenience, we shall always omit the o when multiplying a quater- 
nion by a scalar so that (-2) o Q = -2 Q. 

By analogy to complex algebra, let us define the conjugate quater- 
nion to Q. The conjugation operation will be denoted by ( )*. Thus 


Q* = Q4 - i Q1 - j Q2 - k Q3 . (E6) 


So far all of our definitions have been extensions of those for complex 
numbers as can be seen by assuming Q2 = Q3 = R2 = R3 = 0. Thus the 
complex number system is a subset of the quaternions. It can be easily 
seen that 


Q* o Q = Q o Q* = Q1 Ql + Q2 Q2 + Q3 Q3 + Q4 Q4 . (E7) 


Note that Q* o Q is a pure real number or scalar. With this observation 
we can define the inverse: 

Q' 1 = (l/(Q*oQ))Q* = Q*/(Q*oQ) ; Q*o Q # 0 . (E8) 


Finally then, we can define a division operation: 

Q t R = Q o R _1 . (E9) 
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We can see that Q :Q=QoQ 1 = Q o Q = 1 so that Q 1 satisfies the 
necessary properties of an inverse as long as Q £ 0 . This will be useful 
later. 

We need some additional results and definitions. First we can show 

that 


(Q o R)* = R* o Q* 


(E10) 


If the quaternion Q = Q*. then Q is necessarily a scalar. Also, if Q = 
-Q*, Q is purely imaginary or a vector quaternion. If V is a vector 
quaternion, we shall designate this by an underline as is also used to 
designate a 3-space vector, i.e., (i, j, and k will not be underlined) 


V = i VI + j V2 + k V3 . (Ell) 


For compactness of our notation we shall let 


Q = Q4 + Q , (E12) 

where 

Q = i Q1 + j Q2 + k Q3. 

Thus 

QoR = (Q4 Q4-Q - R) + Q4R+R4Q+Q *,R . (E13) 


The operations * and x are defined as for 3-space vectors so that 

Q • R = Q1 R1 + Q2 R2 + Q3 R3 (E14) 

and 

Q x R = i (Q2 R3 - 03 R2) 4 - j (Q3 R1 - Q1 R3) 

+ k (Q 1 R3 - Q3 Rl) . (E15) 
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For vectors A and B, 


AoB = -A • B + A * B 


(E16) 


In general the product of quaternions mixes scalar and vector parts 
together so that this produce is not very interesting in the study of 
rotational motion in 3-space. The triple product 


V’ = Q* o V o Q 


(E17) 


is more interesting since it does preserve scalar and vector parts of V 
without mixing them. This property is trivial for the scalar part of V 
and follows for the vector part sn.ce 


V* = (Q* o V o Q)* = ~Q* o V o Q = -V. . (E18) 


Hence as noted the triple quaternion product (E17) takes a scalar into 
a scalar and a vector into a vec tor for a ny quaternion Q. Furthermore, 
the length of the vector |v| = / V * V and 


V' • V’ = V' o V’* = Q* o V o Q o (Q* o V o Q) 


= (Q* o Q) V o V* . (E19) 


Equation (E19) indicates the triple product (El") multiplies vector length 
by the factor Q* o Q which is a real number. We note that if Q* o Q = 

1 , vector length is preserved and the vector mapping V + V ' looks like 
a rotation operator. It is a linear operator in that aA + bB^ -> a A' + bB\ 
Restricting ourselves to normalized quaternions which preserve length, 
we see that E17 is equivalent to a rotation of vector V into V’. Since 
we are looking only at normalized quaternions, we can without loss of 
generality, represent Q as 


Q = cos <t>/2 + sin i£>/2 u ; where u • u = 1 . (E20) 

The triple product equation (E17) can be combined with equation (E20) 
to give 
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(E21) 




V' = cos <p V - sin 4>u x V + (1 - cos <f>) u u * V . 


Equation (E21) is the general form of the rotation of a vector V about 
axis u through the angle -<J>. That this is true is seen by examining 
the rotation operation. Clearly any vector along the rotation axis u is 
not changed by the rotation so that if the vector V is broken into parts 
parallel to and normal to u; i.e. 


V = V + V 
— —I —1 


where V 

~ n 


V • u u and V • u = 0 . 

(E22) 


We must also have V’ = V + V ! ; where V’ • u = 0 and V = V’ . 

- -II -l -i - -I -| 

Since and V'^ are normal to u, we can express as 

V’^= xu*V^+yux(ux V A ) . '^23) 

Now V' * V = |V | 2 cosa = -y |V 2 | ■+■ y = cos a and V x v' = 

2 -i- 1 1 1 * 

[ V ( sin a u and thus x = sin u; where a is the rotation angle. Com- 
bining, we obtain 


sin a u x V - cos a u x (u x y^) (E24) 


u u • V + sin a u x v - cos a u x (u x V) . (E25) 

The equivalence between equations (E21) and (E25) for a = -<p is 
established. Thus, the mapping (E17) is equivalent to a rotation 
operator in a vector 3-space. Since the coordinate directions are also 
vectors, we can rotate the coordinate system instead of the vector. 
Rotating the vector through yields the same components V’ as 
rotating the coordinate axes through <p. Thus we can look at V' as a 
new vector formed from V by rotation and expressed in the old axis 
system or as the old vector expressed in new axes rotated relative to 
the old. We have now demonstrated that any coordinate system rotation 
can be represented by a quaternion. Note that if Q satisfies equation 
(E17), then so does -Q. Looking back at equation (E20) tells us that 
-Q corresponds to <p + 360 degrees which represents the same attitude 
4> does. 


i'i 


and 


V- = 


81 



We now look at the time variations of Q. Since Q is constrained 
to be normalized, we necessarily have 

d/dt(Q* oQ) = 0 = Q*oQ+Q*oQ . (E26) 

We see from equation (E26) that Q* o Q = -(Q* o Q)* and hence must 


be a vector. 

Let us define this vector by 



Q* o Q 

= 1/2 03 


(E27) 

■* Q = 

1/2 Q o (o (since QoQ* = Q*oQ = 1) 


(E28) 


We shall see the reason for the 1/2 factor later. When we evaluate the 
rate of change of a vector with time in two reference frames, we find 

V = Q o V’ oQ* + Q o V'oQ* + QoV' oQ* . (E29) 

Using equations (E28) in (E29) , 

V. = Q o [V' + 1/2 (to o V* - V’ o to)] o Q* (E30) 

V = Q o [V’ + w x V] o Q* . (E31) 


Now the reason for the factor 1/2 becomes clear. It is so that we can 
identify w_. Equation (E31) is exactly like the corresponding equation 
for 3-space vectors if we identify fo as the angular velocity of the primed 
reference frame expressed in primed coordinates. This identification 
follows from the fact that equation (E31) holds for an arbitrary vector V. 
Thus, w is identified as the relative angular velocity of the primed axes 
with respect to the unprimed. 

The above discussion completes the basic development of our 
quaternion tools. We now turn to the problem of developing a more con- 
venient notation. The most logical choice which comes to mind is a 
matrix representation. The quaternion Q would logically become 
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- -jr*. • 


Ql 

Q2 
Q 3 
Q4 


(E32) 


Looking back to the definition equation (E4) of the quaternion product o, 
we see that for P = Q o R we have 


" Q4 -Q3 Q2 Ql" 


’ri" 

Q3 Q4 -Ql Q2 


R2 

-Q2 Ql Q4 Q3 


R 3 

-Ql -Q2 -Q3 Q4 


R4 


(E33) 


(See Notes at end of Appendix E.) Similarly, if D = A o B o C, then 


D = AoBC = ABC = ABC 


(E34) 


This result shows that the set of matrices of the form Q have the 
properties of the quaternions and in fact comprise a matrix representa- 
tion of quaternion algebra with matrix multiplication corresponding to o. 
We can also express the quaternion product in the alternate form 


It 4 R3 -R2 

-R3 R4 R 1 

P = R2 -R 1 R4 
-R1 -R2 -R3 

(See Notes at end of Appendix E.) Here the mapping also yields a 
matrix representation except the order of the factors must be reversed. 
Thus, we have 

P = Q R = R Q 

D = ABC = CBA . (E36) 


B o C = C B 


R 1 
R2 
R 3 
R4 



Ql 


Q2 


Q3 


Q4 


= R Q 


(E35) 
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According to these definitions and results, we have 


ABC = ACB = CBA = CAB . (E37) 

1H 

Thus, an interesting and sometimes useful result is that A C = C A. 

Let us now look at our previous work and make use of these new 
definitions: 

V' = Q* o V o Q •* V' = Q* V Q = Q* Q V . (E38) 


Note that we now have a matrix formed from Q which rotates coordinate 
axes and produces V' from V. Linear vector spaces are also represented 
by matrices. The 3-space vectors, V and V’ are related as 


V' = M V (E39) 


Here M is a 3 * 3 matrix which transforms components of V to primed 
coordinates. Referring back to equation (E25) and replacing a by -<p 
we see that 

T ~ T 

M = u u - sin <p u + cos <p (1 - u u ) . (E40) 


(See Notes at end of Appendix E.) The matrix u is the so-called cross 
product matrix and just happens to be^the upper 3x3 formed by 
dropping the final row and column of u. The matrix 1 is the identity 
matrix of the appropriate size to fit the curreqj application. We have 
already shown the equivalence between M and Q* Q. Let us look more 
closely at the latter since it is 4 x 4. We can partition the double-tilde 


or double-bar matrices as 


’ Q j Q 

Q = Q4 1 + 

i 


T 1 

[-Q i Oj 


; Q - Q 4 i + 


” 2* 1 
-Q i 

i 

Q * 

T 1 

L-« ! 

0 


(E41) 


Since Q* o Q = 1, Q* = Q so that Q* Q = 1. Also, Q* = Q(transpose) 
so that 
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4 1 in *inmo4 


mmv* 


* . r Q4 2 1 

Q* Q = 


-2Q4Q+QQ+QQ 

0 


0 

1 J 


(E42) 


Hence M = Q4 2 1 - 2 Q4 Q + 2 Q Q T - Q 2 1. In expanded form 


M 


o 2 > 

yi - Q2 - 03“ + Q4“ 

20)2 01 CJ3 04) 

2(03 01 + 02 04) 


2(01 02 + 03 04) 

-01 2 + 02 2 - Q3 2 + ()4“ 
2(03 .*2 - 01 04) 


2(01 03 02 04) 

2(02 03 + 01 04) 

- 01 2 - 02 2 + 0 3 2 + 04 2 


(E43) 


From equation (E28) 





" 0 

0)3 

- 0)2 

0)1 “ 

1 ~ 

1 - ^ 

1 

- co 3 

0 

0)1 

o )2 

= 2 Q ° w 

= 2 ~ Q 

2 

•jj2 

“ 0)1 

0 

w 3 




. “ 0)1 

- oj 2 

— 0)3 

0 . 


Q1 

Q2 

Q3 

Q4 


. (E44) 


Equations (E43) and (E44) summarize the useful results from our 
discussion. 

We are now ready to consider the question of successive rotations 
applied to a coordinate reference. A coordinate frame rotation is a rigid 
displacement of all the points in the system with a fixed axis passing 
through the origin. Thus, it would seem that several successive rota- 
tions should displace every point except the origin. Let us now consider 
the coordinate frame as a rigid body and determine the most general 
displacement of it which keeps one point fixed. We must first explain 
what is meant by a rigid body displacement, A rigid boc.y displacement 
is one which preserves distances between every possible pair of points 
in the body. The displacement is mathematically represented as a vector 
function f. This function then has two basic properties: 

1) f(0) ^ 0 

2) Jf (r A ) - f(r B ) | = !£ A - r _ B i • (E45) 
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To study this in more detail, we define two additional points and r 2 
together with their images under f, f(rj) and f(i\,). Let = and 

f 2 = f(r 2 >. Let us define unit vectors 


i 




i 



-2 





— 1—1 


(E46) 


The vectors i, j, k and u ^ , u^, u^ each form orthonormal bases for 3 
dimensional space. An arbitrary vector r can be expressed as 


r = x i + y j + z k 


(E47) 


The corresponding f(r) = f u^ 
(E45) can only be satisfied if 


+ f 

y 


u + f 
—2 z 


u 


3* 


Condition 2 of equation 


f(r) = x' u_ 1 *+- y u_ 2 ± z u^ 3 


(E48) 


Thus we can define two functions f+ and f- that both satisfy equation 
(E45) and map r^ into f ^ and into f^. The function f- can be viewed as 

the reflection (x,y,z) (x,y,-z) followed by f+. We are only interested 

in continuous transitions from an initial position to a final position and 
thus reflections must be eliminated since it is not possible to go from 
(x,y,z) to (x,y.-z) continuously without violating condition 2 of equa- 
tion (E45) . Thus continuous rigid displacements can only occur in the 
form f+. This function can be written out as 


f(r) - x u 1 + y u * z u = u i + a j + a k (E49) 

i z ij x y z 


x ill 

• i + y 

u 2 • 1 + Z u 3 • 1 



X Uj 

• j + y 

u 2 • j + z u 3 . j 


(E50) 

X u^ 

• k + y 

* k + z u„ • k 

• 
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Equation (E50) can be rewritten in the matrix form 


a .= M r . (E51) 

The vectors a and r are of the same length and since this must hold for 
all pairs a and r we must have that 


M T M = 


(E52) 


Equation (E52) also implies all eigenvalues of M are of unit magnitude. 
The eigenvalues and eigenvectors of M may be complex so that if M x = 
T* T T 

X x, then x M M x = 1 = A*Ax *x. For 3-dimensional space M must 
have at least one real eigenvalue. Since M is real, its eigenvalues must 
occur in complex pairs. Therefore at least one eigenvalue of M must be 
equal to 1. The value -1 could not be acceptable since it would imply 
M u = -u which would be a reflection and already ruled out. Thus, we 
have that det M - 1. 

The matrix M is now looking very much like a rotation since the 
eigenvector u is an eigenaxis. All we must do now is to determine the 
angle of rotation. Along with the eigenvector u, let us define unit 
vectors v and w such that u, v, w is an orthonormal basis set. Also, 
v:e assume w = u x v. With these definitions we can express the matrix 
M as 


M = u u T + M 2 u v T + i*l v u T + ... + w w T . (E53) 

T 

From the fact that M u = u and that M M = 1 equation (E53) reduces to 

X t 1 x X X 2 2 

M = u u + p(v v A + w w ) + q(v w - w v ) ; p +q =1. 

(E54) 

We can now eliminate the vectors v and w from this equation by use of 
the proper function of u. Thus 

M - u u T + p(t - u u T ) - q u . (E55) 
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This completes the proof thut the matrix M is a rotation matrix. This is 
now obvious from inspection of equation (E55) by comparing it to equa- 
tion (E40) with p = cos <l> and q = sin <p. Thus the most general dis- 
placement of a rigid body (or transformation of a coordinate system) in 
which at least 1 point remains fixed is a rotation about a fixed axis 
i.e., the final orientation can be obtained from the original by a single 
rotation about the axis u through the angle <p (u and <p are determined 
from M) even though the actual motion from initial to final may have been 
more complex. 

What all the previous discussion boils down to is that the product 
of a pair of rotations is itself a rotation. Thus, if Mj and are rota- 
tions about u 1 and u 9 respectively, then Mg = is also a rotation 
through some angle <j) about some axis u„. In fancier terms the set of 

o 

rotations forms a group under matrix multiplication. 

The results of our previous discussions now suggest some new 
notation that may aid us in keeping up with the multiplicity of coordinate 
systems that must usually be dealt with in analysis of spacecraft rota- 
tional dynamics. To remain completely general, let us consider three 
coordinate frames A. B, C. We shall let the symbol [BA] represent the 
rotation matrix which transforms a vector expressed in the A frame to a 
vector expressed in the B frame. 


V (B) 


| BA | V (A) 


(E56) 


For convenience we use the notation superscript (A) or (B) etc. to 
indicate which coordinate frames the vectors are being expressed in. 

If the superscripts are not specified. »? means that the coordinate frame 
is implicit in the definition of the symbol or tb*’t it doesn't matter as 
long as all vectors are in the same frame. V.e aiv ure more interested 
in the rotations [BA] etc. There are three rotation? between pairs*. 

AB , BC, CA. From our previous work 


[CA] = ICB] [BA] . (1*57) 


Corresponding to equation (F57) is a quaternion relation of similar form. 
First, since - [ C A ] V ( ^, we have 


r(C) 


= g 


CA 


o V 


(A) 


Q 


CA 


(E58) 
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■ : ’5*' %t : '<'*jf^f • ,«►- t 


Here, Q„ A is the quaternion corresponding to [CA], Thus analogous 
to equation (E57) we hove 


V (C> = «CB* 0 Q BA* 0 ^ <A> ° «BA 0 ^CB • (E59) 

Interestingly, we see that Q^ A = o Q C jj so that the factors occur 

in reverse order from equation (E57). However, if we use the double- 
bar operator we can multiply in the same order, i.e. 

q ca = S CB Q BA • (E60) 


This result is the one which we wish to use analogous to equation (E57), 
Equations (E57) and (E60) have an easily remembered form and in fact 
behave as if multiplication cunccllcd the terms appearing on the inside. 
This makes it quite easy to construct chains of transformations to any 
desired system. In this notation we see that 


[BA] = [AB] r - [ABf 1 2 ; also 
Q BA = C *AB* " Q AB 


(E61) 


Finally, there are some useful tricks with the new notation we have 
defined. Referring to equation (E28) and adding the subscripts we have 


defined, we have U„, = 0.5 Cj t> . o n n . 

BA BA -BA 


(B) 


The vector is the 

—BA 


angular velocity of B relative to A with components in B. Consider the 
quaternion 


/\ = 1 Q 0 (C) _ 1 q / (O _ (C) 

^CB 2 ^CB : ‘-CB 2 ^CB -^CA -BA 


1 (C) 1 „ /„ * (B) . 

2 Q CB ° ^CA “ 2 Q CB ° \ Q CB ° -BA ° Q CB 


1 „ (C) 1 (B ) 

2 Q CB ° ^CA " 2 -BA ° Q CB 
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(E62) 


or in matrix form 

* 



The utility of equation { IC*J2 ) is most apparent when we use it to compute 
the attitude error of a spacecraft relative to a moving or moveable 

reference. Note that the components of and are expressed 

in different frames. Normally, . 0 ,.^ would come from rate sensors 

( H ) 

which are body fixed while ^ is a commanded maneuver rate which is 

naturally defined in the moveable reference. This equation then allows 
us to use both quantities diroeUy without cither being transformed. 


It often becomes necessary to compute the quaternion correspond- 
ing to a given rotation matrix, i.e., find Q given [BA]. We have 
developed a computer algorithm to do this. 

1) Define matrix 


S 


A 







A is the given rotation. 


2 ) 


S’ 


N »■ S 


. T 


tr A) 1 . 


© 3) 

4) 


5 ) 


I -■ max .S'jj (index o" largest element along diagonal of S’). 


o . - a'./2 / s '.. . 

j ij ii 

Q'- - •<: sgn Q, • 

J J 


1 for i 


< 0 

1 1 for i 

«4 

> 0 


Another useful and perhaps obvious tcehni<|ue is the expression of the 
quaternion resulting from a sequence of Euler rotations (rotations about 
coordinate axes): 
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„,t, w>? wr.Tir - ' 1 ' *■.' ■**'*» v*» - *&*«*■* - 



where S = sin and C = cos. The corresponding rotation is [BA] and 
is given by 


[BA] = 


[d> ]. 

lT n J i 


n 


l ¥i. 


“A 


The vectors u can be any of the three coordinate axes [1,0,0] 

If u = [1,0, 0] T 


[0,1, 0] T or [0,0, 1] T . 


then ij = 1, etc. 


We have 


added the convention that a rotation bracket with a subscript is an Euler 
rotation about the indicated axis. As an example consider the quaternion 
formed when i^ = 1, i 2 = 2, i^ = 3: 



+ 


+ 





the corresponding [BA] is 




Cf 2 C$3 
' C* 2 S* 3 

s * 2 


S4 - 2 CiJ ) 3 + S 4> 3 

S<j J 1 s 4>2 s ^3 + c< J>i ^^3 
“S4> 1 C4> 2 


j - C4>^ S<)> 2 C<j> 3 + S 4> 2 S 4> 3 
I S4> 2 S 4) 3 + S 4> 2 C4> 3 

I CtJ^ C 4) 2 


In this brief exposition, we have developed a number of useful quaternion 
results and notations. This by no means exhausts the possibilities. The 
available quaternion literature does not present the material in an easily 
applicable form and thus this short development is presented to fill that 
gap. 

NOTES : The following notes apply to the previous discussion: 

1) We use Q 4 rather than Q q for convenience. Since these 

quaternion equations will be adapted for the computer and since 0 is not 
usually allowed as a subscript it becomes necessary to use something 
else. We desire to use 1, 2, 3 for the vector components, hence is 
the real part; 

2) The symbol ~ is called double tilde and the symbol = is called 
double bar; 

3) We shall define 



which is the tilde or ci’oss product matrix for the 3-vector Q. 
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APPENDIX F 

STEREO V I SU AL 1 ZATI ON W I THOU T 0 PTI CAL A I D S 

(Cross-eyed Stereo) 


On many occasions in engineering and physical analysis it would be 
useful to be able to sketch in three dimensions. This would be especially 
helpful in the study of the rotational motion of spacecraft. To fulfill this 
wish in many cases, we have used a convenient technique which requires 
no optical devices other than one's eyes and some graph paper. All that 
is required is a stereoscopic pair of images. One additional capability is 
necessary. The observer must be able to cause the lines of sight of his 
eyes to converge; i.e., one must cross one's eyes. The stereo projec- 
tions are formed as shown in Figure FI. The images are reversed and 
viewed as in Figure F2. With a little practice, one can easily learn to 
reconstruct mentally, the 3-dimensional scene from the reversed stereo 
pairs. In Appendix G there are several stereo views of the various 
equilibrium attitudes for the Skylab vehicle. The Skylab is represented 
as an idealized "wire-frame" model which is transparent. The interested 
reader should try several viewing distances (the farther away the page 
the less crossing of the eyes is required and the easier it becomes to 
focus the images) . Squinting may also help as it increases one's depth 
of focus. When one first looks at a stereo pair, one focuses on the 
page and sees two similar but separate images. As one begins to cross 
his eyes, the two images become four. Continue crossing the eyes until 
the interior pair of images come together. Since the line connecting 
corresponding points on the images must be at the same angle about the 
line of sight as the line connecting the eyes, it may be necessary to 
rotate the page or rock the head until these two images become super- 
imposed and seem to merge into a stereo image. This technique does 
require some practice but once mastered it can be very useful for easy 
visualization in 3 dimensions. 

If a computer with plot capability is available to us , we can con- 
struct the necessary stereo projections from a set of points and lines 
that represent the object of interest. We have referred to such a 
representation as a wire frame model because of the appearance of the 
image. Let P be a representative point of the model. Each point P is 
projected into the picture plane S as shown in Figure F3. The point P 
is projected to the eyepoint E and the line PE intersects the picture 
plane S at P' . P' is the projection of P onto S . The set of all points 
P' projected from object points P together with the connecting lines form 
the desired projection. We set up a reference frame in the plane S. To 
do this, we must specify which way is up (so to speak). Let u be a 
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unit vector in this direction and u = u x l is a unit vector in S 

— r — u - 

pointing to the right. We place the origin of the S coordinate reference 

at O. Observe that r = r T , + d l\ where r and r„ are position vectors 

o -E — ~o -E 

of O and E respectively. From the geometry shown in Figure F3, we can 
see that 


P f 


+ (x\ 




dj l 


( -P ' 




From this we can compute 

x p» “ (£p> _ !lq ) ' H r 
y p . = <Ep. - £o> • Hu 


Since l ♦ u = l * u = 0, 
— -- u — -r 


v - ^ • Hr 

y P' = Vpi " -U -*u 


The set of points (Xj M , y^,) plotted conventionally forms the desired 

projection. Size can be altered by scale adjustments. These projections 
are then placed as desired. Also, the values used for d and eye separa- 
tion s are arbitrary and can be adjusted for convenience or eye comfort. 
In reul life s ^ 65 mm and d - 250 mm for comfortable reading; however, 
it may be more comfortable for d to be larger. Some initial experimenta- 
tion with this technique should establish desireable settings. 

We have developed some APE software to produce such images. 

This software was used on our Sigma V computer and is presented for 
the benefit of the interested reader with a background in APL (See 
listing at end of Appendix II . ) . 

The software consists of three functions: 1. STSETUP; 2. STROT; 

and STEREO. STSETUP allows the user to specify the various parameters 
for the observer such as d, s, 7, etc. STROT allows the user to rotate 
the line of sight and the observer position to view the object from a 
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different angle without changing the viewing distance or eye separation. 
STEREO is the function which takes a set of points connected by straight 
lines and forms the stereoscopic projections. Usually these functions are 
driven by another function which is set up to draw the object of interest. 
A function called SKYLAB does this for us in this case. 



Figure FI. Stereo projection. 



Figure F2. Stereo reconstruction by cross-eyed viewing. 
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APPENDIX G 


COMPLETE TEA SET 


(Cross-eyed Stereo Doubles) 


Figures G1 through G12 show all the available TEA's for an altitude 
of about 200 km (108 n.mi.i. The nomenclature is derived from the 
general directions of the vehicle X axis (first letter) and the vehicle Z 
axis (second letter) with respect to the L system. 

The axes can point 

N north 

S south 

U up 

D down 

F forward 

B backward 

7 

The top set of 3 x 3 numbers in each figure is the direction 
cosine matrix from the L system to the A system. Further the BAR 
angles and the x's are given. The vehicle Z axis beta is the elevation 
angle of the vehicle Z axis with respect to the orbital plane. 

The point of view is from behind the vehicle, slightly south of the 
orbital plane. The axes shown are these of the L system. The Y axis 
points to the north, the Z axis points up, and the X axis points forward 
(into the paper plane). The vehicle is moving from the lower left to the 
upper right. 


7. In APL computer language the minus sign is a superscript minus 
(to distinguish it from the subtraction operation). 
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~1H momentum in Z axis. 







0668563 0.117329 0.99004 



Figure G3. TEA - UN or T121G. 



1666563 6. 117329 6.99684 



Figure G4. TEA 










Figure G7. TEA - BN. 
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Figure G8. TEA - BS. 
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Figure G9. TEA - DN. 
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Figure G10. TEA - US. 



Figure Gil 
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Figure G12. TEA 






APPENDIX H 

TEA SIMULATION AND GROUND SUPPORT 


To assess the feasibility and refine the techniques of TEA control, 
we had to develop a simulation. Simulation tools already existed to some 
extent as a result of our work on EOVV control. These were modified 
and extended for TEA. The details of the theoretical background for 
most of the development of TEA control are presented elsewhere. We 
shall concentrate on a few additional areas. 

IBM also had Skylab simulations left over frcm original Skylab days. 
These were used to verify the control schemes we developed and also to 
verify the resulting new flight computer code. In addition, they were 
used to confirm our ground procedures and to predict TACS propellant 
usage [2]. 

The computer and the programming language are of fundamental 
importance to the success of this effort. The computer used was a 
Xerox /Honeywell Sigma V. We were most fortunate to have ready access 
to this machine together with timesharing hardware and software. The 
language used was APL, a high order, interpretive language which 
offers high computing power together with a relatively simple syntax. 

Since the language operates on the tensor or array level, it is a highly 
compressed language, allowing one or two lines of APL code to be 
equivalent to dozens of lines of Fortran. TEA control and EOVV control 
could not have been successfully accomplished in such a short time and 
low manpower effort without the efficient use of APL and the Sigma V. 

As an example, we were able to assemble virtually from scratch within 
one week a complete dynamic simulation of Skylab EOVV with 3-axis 
rotational dynamics and CMG control. 

For TEA we modified the Skylab EOVV simulation primarily by 
redefining the momentum desaturation technique. Since we hud to 
simulate many orbits of TEA operation to assess its behavior, we had to 
adopt a somewhat modified approach. To speed running time of the 
simulation functions, it was desireable to lengthen the integration step 
being used by an order of magnitude. This was done by sacrificing the 
body rotational dynamics. A slow running version of the TEA simulation 
was maintained with rotational dynamics in it to spot check this fast 
version but was not used much because it proved too slow. We assumed 
the vehicle control was perfect in maintaining the vehicle at its 
commanded attitude and thut small maneuvers were made instantaneously. 
These changes in attitude were made at the beginning of the integration 
period and afterwards only gyro drift caused attitude changes relative 
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to the local vertical reference. To improve accuracy the vehicle torque 
was computed twice during each cycle: at the beginning and at the end 

of each step. The attitude maneuver made to desaturate the angular 
momentum was assumed to be made at the beginning of each step. The 
torque was calculated for the new attitude at the beginning. Since the 
position may have changed at the end of the interval it is calculated 
again and averaged. With this technique we were able to obtain good 
results using a 300-sec integration step which was also the TEA desatura- 
tion interval. This allowed us to simulate 1 day of orbits in about 5 min, 
if we had priority demand on machine time. This was provided in the 
month prior to initiating TEA control and allowed us to develop and 
refine our operating procedures including those used to initiate TEA 
control. 

In addition to the simulation, we developed many analysis tools. 

As can be noted by examining the listings at the end of this appendix, 
there arc many functions including simulation functions, support func- 
tions, plotting functions and auxiliary functions. This set constitutes 
the library that was used in developing the TEA control schemes and in 
maintaining TEA control once it was established. Various subsets of 
these functions were used, depending on the application. 

Attitude dynamics in the slow simulation was simulated using the 
Newton-Euler rigid body rotational equations 



*Ek 


^Vk 


H 


Tk 


The total system angular momentum H was made from the sum of 

vehicle and CMG momentum. The CMG momentum was computed from the 
gimbal angles which were determined by integration of the rate commands 
coming from the CMG steering law of Reference 5. These rate commands 
were assumed to be provided by the CMG gimbal dynamics except for 
gimbal rate limits and stops on the angles. The vehicle position was 
determined by integrating the vehicle angular rate. This rate was 
obtained from the vehicle angular momentum by use of the inverse 
moment of inertia tensor or in dyadic form, 


--Vk =Vk ' ( -Tk “ ^CMGk ) 
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In the fast simulation, measured body rate was assumed controlled 
to orbital rate and the true rate was then determined by subtracting 

the gyro drift and dividing by the scale factors as in 


-Vk 


= Isk 1 


( ^o ‘ 


-Dk } 


The CMG momentum was then 


^CMGk " -Tk " =Vk * -Vk , * 


The total torque from external sources, T^, 


is made up of aerodynamic 


and gravity gradient torques as developed in previous appendices. The 
position variables used were quaternions as explained in Appendix E. 


After the TEA control was developed but before it was initiated 
in flight , the simulation was modified to convert it to a tool for the 
training of ourselves and the JSC and MSFC flight controllers who would 
be interfacing with the vehicle when the time came. We were able to 
run many practice sessions prior to actual TEA initiation. We had not 
been able to do this for EOVV and paid the price for this on several 
occasions. As a result of these training and practice sessions, we 
learned much on how to maintain TEA even in the face of large simulated 
errors between data in the flight computer models and the real vehicle. 
As a result we put ourselves in a much more favorable position to 
enhance the probability of successful initiation of TEA control. We 
learned how to prebias the rate gyros so that the strapdown would 
intentionally drift in a favorable direction. This gain in experience was 
crucial to overall success. 


On the following pages, the simulation functions are listed. The 
first four (ACR, AFMT, AGRF, AWM) are proprietary and cannot be 
listed. Only AGRF and AFMT are used by the simulation to support 
graphics and data formatting. ACR and AWM are used only for listing 
all functions in an APL workspace and are not otherwise required. 
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DATE 


WORKSPACE GLAESE TEASUP 
TIME- 10 39.7 DATE = 11 9 79 


) FNS 

NCR NFMT NGRF NWM l ADEG AECONST AETORQUE AFEA 

AFQ ALTFDOY ANYROUT AT AN ATMDC AXES BAIRSTOW BAR 

BARANGS BAhAV BARS BOX BPRBANK CHI CtiRPOL CP CPA CPASUP 
DEE DA DARROW DATE DBANK DBANKPAR DBANKSUP DEAR 

DbARA DB Ah AFAR DBARASUP DCD DDDDER DECMG DEGA DEGQ 

DELG DENS DBS AT DET DFDOY DFHEX DFOCT DPS DH 

DOYFALT DPR DPWR DPWRA DR DRAW DRAWN DRAWNS DRIFT 
DRO DTQ EAFA EDITQBL EIGVEC ENC ERASE EUL EUL 1 EUL 2 
EUl 3 FCF F1LLTADS FINDQR FIX FORMAT 8 FPWR GACCEL GBP 

GbPA GBPASUP GETSLOPE GGTORQUE HBARS HELP HESSEN 
HEXED KMAhk HOME IC IEGEN JACOBIAN LIMC LIMM 

LOUP LSTALL MAG MARK NAVVP NOM NORM NRMLIZ OCTFD 
OhBMARk OUTPUT PAR PER f OWER PRIC PUT PWR ANGLE QBLCHECK 
QbLSET QbLSETSti QDb '.DEG QDT QFA QIN QRALG QUAD QUADM 

QUANTIZE RAN 1C RANPAh RDTAPE REF REFRESH REPRIC RUN RUNON 

SCALE SEEK SE TSTATIONS SETUP SFD SKYLAB SLUPDATE SPIN 

SPLIT SSANGLE SS ANGLE 0 STAXES STEREO STPAR STRAPIS STROT 

STSETUP TACS TEAEVAL TEXPERT TOP TRANSLATE TRQ UPDATE VS 

WHATSCALE WHATW1NDOW WINDOW X XA XH XPLOT 


VLSTALLiV 1 V 

V LSTALL NLP ; 1 ; FN S ; F N ; NWM ; NFMT ; EC R ; N LR ; N LF ; N FL ; LAV 

Cl] a LISTS AUTOMATICALLY ALL UNLOCKED FUNCTIONS IN A WORKSPACE 

L2l ft ALLOWING NLP LINES OF LISTING ON EACH PAGE. 

C 3 ] SLLSTALL* \U\SNLSTALL + 2 2 p£b , 3 5 , Lb , 4 

[41 FNS* 3 NWM 6; AJ/M«-14Yb ; AJ'A17'*14 tO; AC2?*-14t5 

[5] L : NLR*-NLP 

[0] £3:2 *-FNS i 2 T 2 1 

C 7 ] -*■( 1 •-/)/ Q 

[ 8 ] FN+FNSi \ 7-1 1;NLF+0 

L 9 1 NLP*- 3 + + /(?Tl3) = l NCR FN 

[10] -*• ( /.' /• /•' > /V L h ) /M 

[11] 2 T 1 3 

[12 1 c ' V * , EN . ’ [ (.] I V ‘ 

[13] LL : h NS +■! i FNS 
[ 14 ] +L3iNLR<-NLR-NLF-l 

[15] M : * ‘ ,2 T 2 2 

[16] +2,;^ 

V 

L o T A L L bO 
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V C+A £ B 

Cl] A NORMALIZED CROSS PRODUCT 

[2] C+(Al2 3 1 3 3 1 2 ] ) -/[ 3 1 2 3 *BC 2 3 1] 

C 3 ] C-CH+/C*C)* 0.5 

V 

7 R+ADEG A;HANG;Q;S 

C 1 ] a MAKES DIRECTION COSINES OUT OP EIGEN ANGLE IN DEGREES . 

[2] W*-4p0 

C 3 ] Qlkl + '/OHANG+O .*>*RFD*MAG A 

[4] S+10RANG 
£53 $Cl 2 2)+S*A*MAG A 

C 6 ] R+AFQ Q 

7 

7 AECONST 

Cl] a COMPUTES AERODYNAMIC CONSTANTS REQUIRED FOR AERO TORQUE CO 
MPUTATIONS . 

C 2 ] KTAE+0. 5X79. 46*10. 06 

C 3 ] DB-DENS ALT+RMAG-REA 
[4] +{IBM-l)/L 

C 5 ] RHOAV+DBl 1 ] i KBULGE+DBl 2 ] *SEBKFL*C 
C 6 ] +0 

[7] LiRHOAV+RHOBASE*\+Q.B*AMP*2*B*2 

7 

7 AETORQUE\U \CMA\CNB\SNB 

Cl] a COMPUTES AERODYNAMIC TORQUE FROM RELATIVE VELOCITY AND ALT 
ITUDE . 

[2] VAE+1000*VEK+ ,*RD-OMGEK X R 

[3] U+VAEiMAG VAE 

[ 4 ] AANG+^OUll] 

[5] RANG+U121ATAN 4/ C 3 ] 

[6] CMA+2OAANG*0 , tlO 

[7] CN b+2oRANG*0 , 1 10 

[ 8 ] SNb+lOhANG*\10 

[ 9 ] AECOEFF+(CNB+ .*(CMA+ .*AN)) +SNB+ . *(CMA+. *BN ) 

CIO] -^(IBM-D/Ll 

C 11] RHO-*-( 1 - KbU LG E*2O0M0T-TH BULGE +01 ) *RHOAV 

[ 12 ] +L 2 

[13] LlibULGE+\+AMP* (1+2O0M0T-THBULGE+O1 )*3*2 

[ 14 j RHO*-( BULGE*hH0BASE* GEEKFL* 0 ) +RHOAV*SEEKFL*l 

[15] L2iTAEK+KTAE*AECCy??*RHO*VAE+ ,*VAE 

[16] TAEK<-{Z*TAEK)-{LCG+ECGNOM)X FAEK+-(*10 .06 )*3+TAEK 

7 

7 R+I AFF.A A 

[1] a ROTATION MATRIX FROM EULER ANGLES IN THE SEQUENCE I GOING 

[2] A FROM RIGHT TO LEFT FROM THE OLD TO THE NEW. 

[3] A-U[3]£'i/£ 4[3] ) + .x(lC2l£&£ Al2}' + .*IlllEUL >1[1] 

7 
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V A+AFQ Q 

[1] ft DIR COS MATRIX FROM QUATERNION 

[2] A*- 3 3tiU {A*-(S)QDT Q) + .xQDB Q) [4;4] 

V 

V ALT*- ALT F DO Y DOY \XI\YI\I 

Cl] fl INTERPOLATES TABULATED DATA TO ESTIMATE ALTITUDE FROM GIVE 

N DOY . 

[2] 

C 3 ] XI~*TAbALTNM 

[4] I + 1U+/P0Y* .>TABDOY)llO 

[5] ALT*-*\ILll + (DOY-YIlI] ) * (XI[I + l}-XI[n)*TIlI*lJ-II[I] 

V 

V /S-Ht ANYhOOT Y;I;N;P;Q 

[1] ft YIELDS POLYNOMIAL QITH SPECIFIED ROOTS WITH REAL PART X 

[2] n AND IMAGINARY PART Y. 

[ 3 ] A-\l; l*l;N*-pX;X*-,X 

L4l LI :-*L2* iQ-y[ 2 ];y^,y 
[ 5 J p-*- - 2 * A l l } ; A' **■ ( A' l i ]* 2 ) + y[ 2 l *2 

[ 6 ] A*- (A , 0,0 ) + (i'x ( 0 ,^, 0 ) )+fi*( o,C* /l) 

[7] +L 3 

[8] L2:/5^(-A'LJj x ( 0»/4) )+/l*0 

[9] LZ:1*1*\ 

[10] +2,1* i (ISiV ) 

V 

V A*-S ALAN C 

ill « 4 QUADRANT INVERSE TANGENT FUNCTION 
L 2 ] ^-(-20(7* ((5*2) +5*2) *0.5) *5* 15 

V 

V ATM DC; l KEEP; RQEKE 

[ 1 j ft SIMULATES THE SKY LAB ATM DC . 

[ 2 1 +(?<TN1 /;) / L 

L 3 1 1 KEEP* 1 Khk ; E?Y KP+-ETYK ; BET A P*-BET A ; NUZEKP-NUZEK 

[ 4 j I (JEN 

U3 EKE ■*-(%■ i k A A ) + . m A U L 3 MJZRK) + . * (EUL1-RFV*BETA) + .*(EUL1 RFDx 

h E T A /•)♦.*( A U L 3 - .*• 5 : : F. EP)+.*IKEKP 
i o i ■inu*-T*( (or i - ( (or) it/wr;r)-a7 , yA:-5fy/(P)*5W5 
I 7 j TM J .0- 2‘:V ; /, + ?0/i /J x 0 . & > ( TMID-T ) * 75P5 

( H ] A 0 t, E r, i, A /- A A 

I 9 1 t* l y A A * ( ^ - /■ 4 ’/ A A ) + .*QIN LQ E K E 

l 1 u J / : Th u*- Hi DM x A A [ ; 2 1 

[ 1 1 i ETb-'-i ET*-UCMGc\ • ti ) +-EV+1 V *■ . *0MV 
112] DEC A T 

[13] TAGS L T + ( l ii ) * IV + , *0MV ~THD 

V 
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7 Q+AXES 

Cl] fi+-5 AGBF ' * fl DRAWS X AND X AXES FOR GRAPHICS. 

7 

7 RTS+BAIRSTOW i4;AI';M;B;P;0;K;e7iS;£;A0;B;P;AP 
Cl] -►Ox\l>!V«-p4>»-l+i4U+,4 ,£+lpl£-5 ,P-+l 1 0p0«-2+M+0,ffrS+-i0 
C 2 ] -+(2 3 = pC^l 0,\-3+pB-0 1, i/V ) / 14 13 

C 3 ] BC2]«-4[2] + (-PxlpB)-0.0pBCl]+-i4Cl]-P,0pW«-3 

C4] ‘►4xiA?>*l+^+p,BCX]^C^] + (-PxBCJC-l])-Cxfl[if-2] 

C 5 ] O[2>B[2] + (-PxlpO)-0,0pO[l]-+B[l]-P 

C 6 ] •*‘6xiB>«;^+p,CC^]^BCJ]+(-PxC[J-l])-QxCCJ’-2] 

C 7 ] ■*■9 x i B * 2 + p , Q+C CB-ll-BCB-1] 

C 8 ] +1 0 , P+-P+EP+- ((B[B-l]xO[B-2])-B[B]) *P+* (£ [B - 2 ] *2 ) -£ 

C 9 ] P«-P+AP**-( (B[B-l]xC[B-2] )-BCB]x<7[B-3])*IMCCY-2]*2)-£xOCB-3] 

CIO] O«-4 + AO<-((BCB]xO[B-2])-B[B-1]x£)*0 
C 11] -►( (200<M+M+1) ,B< ( | A<2) + I AP)/16 3 

C 1 2 ] +•( (1=B) ,3<;B-B-2>/14 3,0p>l*-B,0pP+-3-0<-2+M<-0,0pBrS+-B2’S t 0£/>iD P.fl 

C 13 ] **>0 xRTS+RTS % QUAD 2 p >1 

C 14] a LIBS DELETED HERE 
C 15] +QXRTS+RTS t (— XCl]),0 
C 16 ] *SL0V OP NON CONVERGENCE * 

C 17 ] ft */* CONTAINS THE COEFFICIENTS OF THE N'TB DEGREE POLYNOMIAL 
C 1 8 ] A SUBROUTINE 'QUAD' CALCULATES ROOTS FROM QUADRATIC FACTORS 
C 19 ] a AN EXAMPLE WOULD BE SOLVING THE ROOTS FOR THE POLYNOMIAL 
C 2 0 ] A U**5) + (2xA>*4)-(4x***3)-(20x***2)-(33xjn-18 ; INPUT WOULD 
BE 

C 2 1 ] A BAIRSTOW 1 2 *4 -20 -33 -18 

C 2 2 ] A PROGRAMMED BY JOHN RIC H ARSON , 0 .C ,C . ,0.5.7, SPRING ' 7 0 
7 

7 ROUTES Ah RINiSiCiQliQ2iQ3iQH 

Cl] A COMPUTES THE BANK , ATTACK AND ROLL ANGLES FROM QUATERNION 
SET. 

C 2 ] Ql+S+RIN ; S+ipRIN-*- t RIN ) *4 
[3] Q2+S+SIRIN 
C 4 ] 03-+S + (2x5) 4BJ/V 

C 5 ] 04- (-S)tBIB 

C 6 ] ROUT+DFR* -2O(0l*2)+(-02*2)+(-03*2)+04*2 

[7] SOi/I-SOi/I, [0.1 ] (LFBx( (01x£2) -03x04 )/42MB(0 1*03) +02x04) 

[8] ROUT+(DFRx( (01*02 ) +03*04 M2MB( £ 2 x 04 ) -01*03 ) ,LllR0UT 
7 

7 BARANGS 

Cl] A LIST'S Bi4P ANGLES CALCULATED BY BAR FUNCTION IN TABULAR FOR 
M. 

[2] *14,318.2* LFMT^i. i-lipOBL) , [ 1 ]Bj4P -4 0 + OBL 

7 
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7 BARAViA 

Cl] a ALLOWS GRAPHICAL DETERMINATION OF AVERAGE BAR ANGLES AND 
[2] A AUTOMATICALLY CALLS TB XPERT . 

C 3] 'MARK ANGLES IN BAR ORDER * 

C 4] A+ DC2] 

£5] A+A % DC2] 

£6] A+-A ,0£2] 

£7] 12 1 AREA RFD*A 

£8] T EXPERT 

7 

7 BARS INTiTiSCiTSTOP 

£ 1 ] A DRAWS VERTICAL LINES ON A PLOT AT SPECIFIED INTERVAL. 

£ 2] SC+WHATSCALE 
£3] T+INTnfSCiUilNT 
£4] TST0P+INTxlSCL2]*INT 
£5] L'.MARK T 
£6] +(TSTOP*T+T+INT)/L 
V 

7 Q+BOX 

£1] fi«-9 BGRF * * a DRAWS A BOX AROUND GRAPHICS WINDOW. 

7 

7 BPRBANK 

[1] A BUILDS ARRAY PRBANK USED TO SIMULATE ASAP DATA. 

£2] PRBANK+-3 *18 +BAR FRQAL 
£3] PRATT+PRBANK121 ] 

£4] PRROLL+PRBANKlZ’^ 

£5] PRBANK+i'lS+PRT) .lO.HPRBANKlhl 
7 

7 RO+CHI RI;SiCiQl}Q2iQ3iQk 

£1] a COMPUTES CHI ANGLES FROM INPUT QUATERNION ARRAY. 

£2] Ql+S+RI ;S-*-( pRI+ ,RI ) *4 
£3] G2-*-SiS+RI 
£4] QZ<-SH2*S)IRI 
£5] 54<-<-S)+J?I 

£6] R0<-(2x(QlxQH) -Q2*QZ)ATAN ( -51 *51 ) + ( 52*$2 ) + ( -53 x$3 ) +54*54 

£7] tf5+l?5,C0.1] (2* (52*54) -51*53 )y42MJH51*51) + (-52x52) + (-$3x53)+54* 

54 

£8] RO+DFRxRO , £ 1] * 102 x (51*52) +53*54 
7 
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V P+CHRPOL AtAKiKiKiSiS 

C 13 A THIS FUNCTION USSS LEXERIER'S METHOD TO OBTAIN COEFFICIENTS 
C23 a OF THE CHARACTERISTIC EQUATION OF MATRIX A. THE COEFFICIEN 
T 

C 33 A OF THE HIGHEST POWER IS NORMALIZED TO UNITI. 

[43 A WRITTEN BI HOWARD B. WILSON , U. OF ALABAMA 

C 5 3 »‘ 0 pX*l* 0 xf«ltp 4 X 4>4 

[63 LI tP+P, (♦£«'£+l)x(*lfS)~*/P*$»l*5**S t +/ ,E*AK+AK+ ,*A 

[73 +L1*\K<N 

[83 P«-l,-P 

V 

V RX+CP EAiSBEKFLiRiRDiVLKiOMGBKiFAEKiTABKlTGGKiTOiFOlLCG 
[13 A ESTIMATES CP OFFSET FROM FLIGHT DETERMINED EQUILIBRIA. 

[ 2 3 SSEKFL+l ; R+ (RE A* ALT) *LEK [ 3 ; 3 5 RD+OMOK* (REA+ALT ) *LEKl 1 ; 3 ; VLK +1 

2 1 AFEA BA 

[33 0MGEK+0MGK*(2OLAZK)*LEKl2il 
[43 TO+TRQ VLKiLCG+Q 0 0 
[53 FQ+FAEK 

[63 FJr«'-A< 7 G«-(*F 0 +.*FQ)*F 0 X TO 

[73 'NEW SLOPE * 

[ 8 3 U+NEWSLOPE+GET SLOPE EA 
[93 * CP OFFSET *iRX 

7 

V RX+CPA ; SBBKFL iR\RDi OMGEK iFAEKiTAEKi TGGK ; TO s FO ; ACG ; VLK ; 

[13 a ESTIMATES CP OFFSET FROM SIMULATION DATA. 

[ 2 3 SFFKFL+l { FWi4G«-M>lG MFFd-MLr ) *LEKl 3 ; 3 ; RD+OMOK* (RE A* ALT )xLEK 1 1 
; 3 ; FLtf«-FL 

[ 3 ] 0 MGEK+ 0 MGK*( 2 OLAZK)*LEKl 2 O 
[43 TO+TRQ VLKiECG+Q 0 0 
[53 FQ+FAEK 

[63 RX+- LCG*-( *FO+ , *FO ) *FO I TO 

[7] SLOPE ’ 

[8] O«-FFVSL0PF«-GFrSLGPF £>i«-l 2 1 FiiFd FL* 

[93 * CP OFFSET '\RX 

7 
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7 ACP«-CPi*$i/P ; 7L ; RBOSIG j OHOK | ALT t ETLNOM 
Cl] a GROUND SUPPORT FUNCTION USED TO ESTIMATE CP OFFSETS FROM F 
LIGHT DATA. 

C2] a THIS FUNCTION CONVERTS DATA TO PROPER FORMAT FOR CPA . 

[3] 'INPUT ALTITUDE IN NAUTICAL MILES' 

C 4 3 ALT+KFM* 0 

C 5] 'INPUT ASCENDING NODE PERIOD IN SMS' 

C6] OMOK+ 02*24 60 60J.0 

£7] 'INPUT ESTIMATED SOLAR SIGMA' 

C8] RHOSIG+U 
£9] 'INPUT ETLNOM' 

CIO] ETLNOM+0 

C 11] * INPUT ESTIMATED L 2 STRAPDOWN OFFSET IN DEG' 

C 12] ETLN0M+ETLN0M+ . *EUL2 RFDxO 

[13] ' THE TRUE MOMENTUM BIAS IS '; ETLNOM 

C 14 ] • INPUT QBL ' 

C 15 3 VL+AFQti 

C 16 3 ‘THE BAR ANGLES ARE ' iDFRx l 2 1 EAFA VL 
C 17] ECP+CPA 
V 

7 DEE 

[1] 1 DRAWN' PREE ' A DRAWS EE VS TIME PLOT * 

7 

7 DA 

Cl] 1 DRAWN ' PR A ' a DRAWS STRAPDOWN OFFSET FROM SIMULATION VS TIME 
PLOTS. 

7 


7 S D ARROW A 

Cl] STEREO FP+A+.xARROWxS a DRAWS A STEREO ARROW FOR A 3-SPACE VE 
CTOR. 

7 


V DATE 

[13 ERASE* ERASES PAGE AND LABELS PAGE HEADING. 

[2] * WORKSPACE * ,(l29).5 EWM 1 

[3] * TIME = ' ; 24 60t( L (120 ) *360) *10 ; * DATE = *;100 100 100T125 
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£l] 

£2 ] 

C 3] 
£4] 

C 5 ] 

C 6 3 
C 7 ] 
£ 8 ] 
£9] 
£10] 
£ 11 ] 
£ 12 ] 

£13] 

£14] 

£15] 

£16] 

£17] 

£18] 

£19] 

£ 20 ] 

£ 21 ] 

£ 22 ] 

£233 

[24] 

£25] 

£26] 

£27] 

£28] 

£29] 

£30] 

£31] 

£32] 

£333 

£34] 

£35] 

£36] 

£37] 

£38] 

£39] 

£40] 

£41] 

£42] 

£43] 

£44] 

£45] 


DBANKiPRTiAViTSTOPiHiSiCiASP{BSPiPOiFC\FSiS2iC2\SCiCI\SI; 

PRTF 

A ESTIMATES ST RAP DOUR ERROR PROM SIMULATIOS DATA SIMULATIRG 
ASAP DUMP. 

R+qPRT+PRBARKZ 1 t ] 

PRT+DPR * 0 MO * PRT - ( PRT C 1 ] - TORB I PRT £ 1 ] - TMID ) 

TSTOP+TOP PRT 
SI+10RFD*PRTiCI+20RFD*PRT 

F0++/PRBARK £ 2 s ] { FC++/PRBARKZ 2 1 ] PS++/PRBARKI 2 ; ] *SI 
C++lCI\S++tSI\C2*-+/Cl*CI\S2++/SI*SI\SC++/SI*CI 
MAT+$(R % C t S) ,£l](C,C2,SC).C0.1]S,SC,S2 
C+MAT+.*{FQ,PC,PS) 

AV+CllliAIIP+Ci2liBtiP+Ci$') 

1 DRAWRS t PRBARKl2i')-AV' 

0 DRAW(C+(ARP*20RFD*PRTF)+BRP*l0RFDxPRTF)VS PRTF* 30x0, \L 
T STOP * 30 

'DETRL WRT MIDRIGST * ; (-ARP) t 0,BRP 
$ • 

* SUGGESTED DETRL * \{EUL2 OMOT+OMO*LT)+ ,*(-ARP) .0 ,BRP 
• • 

'TRB AVERAGE BARK ARGLB IS ' \AV 

* TRE AMPLITUDE IS * j ( (ARP*ARP)+BRP*BRP)*0 . 5 

* TRE PRASE IS * ;360 \DFR*(- ARP) AT AR BRP 

* I 

DBARKPAR 

t t 


* 15 ,F9 . 1 * LFMT 13 2+PRTF, £l . lIC+AV 
F0++/ PRATT ; FC++/PRATT*CI ; PS++/PRATT*SI ; D 
C4-MAT+.*F0,FC t PS 

A V+C £ 1 ] ; ARP+C £ 2 ] ; BR P+C £ 3 ] 

1 DRAWRS' PRATT- AV' 

0 DRAW(C+(ARPx20RFDxPRTF)+BRP*10RFD*PRTF)VS PRTF 
* TRE AVERAGE ARGLE OF ATTACK IS * ; AV 

* TRE AMPLITUDE IS * ; ( (ARP*ARP)+BBP*BRP)*0 . 5 

* TRE PRASE IS ' j 360 I DPR* (- ARP) AT AH BRP 


DBARKPAR 


*15 ,F9 . 1 ' LFMT 13 2 + PRTF . £ 1 . 11C+AV 

F0+-+ / PRROLL;FC+‘+ /PRROLL*CI ; FS+-+ / PRROLL*SI 

C+MAT+ . *FO ,FC ,FS 

A V+C £ 1 ] ; ARP+C £2 3; BR P+C £ 3 ] 

1 DRAWRS' PRROLL-AV' 

0 DRAW (C-*-(ARP* 20 RFDxPRTF)+BRP*l 0 RFDxPRTF)VS PRTF 
* TRE AVERAGE ROLL ARGLE IS ' \AV 

* TRE AMPLITUDE IS ' ; ( (ARP*ARP) +BRP*BRP) *0 . 5 

* THE PHASE IS * ; 36 0 | DFR* ( -ARP) ATAR BRP 
DBARKPAR 


*15 ,F9 . 1 ' LFMT 13 2*PRTF, [1 . UC+AV 
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V DBANKPAR 

[1] A PRINTS DATA BLOCK FOR STRAPDOVN ERROR ESTIMATE 
[23 * * 

[33 * * 

[43 *0 DESIRED UPDATE Q,I3 ,0:0, 12 ,0:0.12 ,0:0, F5 .2 ' LFMT 1 4p 

DFS T 

[53 *0 MIDNIGHT TIME 0,13, 0:0.12,0:13. 12, Q:0,F5. 2* LFMT 1 4p 

DFS TMID 

[63 *01 REV DELTA TIME 0,I2,Q:Q,I2,0:0,F5.2' LFMT 1 3p UDFS 

TORS 

[73 *0 FIRST QBL TIME 0,13 ,0 :0,I2 .0:0,12 ,0:0. F5 . 2 * LFMTt'k it 

QBL 

[8] *0 LAST QBL TIME 0,13 ,0:0,12 ,0:0, 12 ,0:0. F5 . 2 * LFMT*'* 

QBL 

[93 * * 

7 

V DBANKSUP \ N AME \ LT \0M0 lOMOT \TGMTO\GMTQiQAL' % T iPRQALiPRT iPRBANKi 
TMID ; TORB ; PRATT ; PRROLL 

[1] n GROUND SUPPORT FUNCTION USED TO ESTIMATE- STRAPDOWN ERROR F 
ROM 

[23 a ASAP DATA. THIS PROGRAM SETS UP DATA IN PROPER FORMAT FOR 
DBANK. 

[33 * HAS QBLSET FUNCTION BEEN EXECUTED WITH LATEST DATA?* 

[43 ‘ ENTER Y OR N FOR YES OR NO ' 

[53 +CYN' en)/LI t LN 

[63 LY:' INPUT GMT OF DESIRED UPDATE IN DHMS ' 

[73 T*- 365 24 60 60lD 

[83 * OBTAIN FOLLOWING FROM CIRCLE CHART ' 

[9] • INPUT MIDNIGHT GMT OP QBL SET ' 

[103 TMID+SFDU 

[113 'INPUT 1 REV DELTA TIME' 

[123 TORB+SFD 0,0 

[133 'TO GET PLOTS OF BANK , ATTACK AND ROLL ANGLES ' 

[143 ’INPUT CARRIAGE RETURNS' 

[153 PRQAL+* 4 0 iQBLift 
[163 PRT+SFD 4 0 +QBL 
[ 173 0M0*-o2*T0RB 

[183 OMOT+OMO*TORB\T-TMID 
[193 LT<- 0 

[203 PRBANK+PRT, [1] (BAR PRQAL) 

[21] PRATT+PRBANK13 ; ] 

[223 PRROLL+PRBANKlml 

[23] PRBANK +* 2 0 +PRBANK 

[24] DBANK 

[ 25 3 +0 

[263 LNi* EXECUTE FUNCTION QBLSET BEFORE USING THIS FUNCTION' 

V 
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7 DBARi INDEX ; IHAX ; ST lit PRBAR 

Cl] a DRAWS BAR ANGLES, SELECTS PROPER ORBIT TO PLOT . 

C 2] 1 DRAWN* BAR (4 % ILEG) +PRQAL’ 

C 3] SHARK 114 

C 4 ] SHARK 93 

£5] SHARK -98.5 

[63 1*0 

£73 IMAX++/ILEG*PRSTZ2 ; ] 

[8] ST<-(-IMAX)*PRSTZ2il 
£93 Ltl+I+l 

CIO] INDEX+(STZn-l) ,STZll »Sr[I]+l 
. C 11 3 0 DRAViPRBAR+BAR PRQALZ lINDBXDZlilVS PRTZINDEXl 

£123 0 DRAW PRBARZ2;lVS PRTZINDEXl 

L 13 3 0 DRAW PRBARZZi 1VS PRTZINDEXl 

. £14] +(I<IHAX-l)/L 

£153 1*1+1 

£16 3 IMSJI>*l + (Sr[I]-l ) ,5TCI].5rCl] + l 
C 173 0 DRAW(PRBAR+BAR PRQALZ ilNDEXDZlilVS PRTZINDEXl 

C 18 ] 0 DRAW PRBARZ20VS PRTZINDEXl 

£193 0 DRAW PRBARZZilVS PRTZINDEXl 

V 

V DBARA 

Cl] a DRAWS BAR ANGLES FOR ALL SIMULATION DATA, 

£23 1 PRAWNS’ BAR PRQAL' 

[3] DBARAPAR 

[4] BARANGS 

V 

V DBARAPAR 

[1] a DATA BLOCK PLACED ON BAR ANGLE PLOTS, 

[23 2T 1 3 13 

[3] *Q MIDNIGHT TIME 0, J3,a:D,I2,CI:D,I2*CI:C],F5.2*AflVr 1 4p 

DFS TMID 

[4] * Dl REV DELTA TIME 0,12 .0:0,12 ,D:0*F5 . 2 * LFMT 1 ZqUDFS 

TORB 

[5] ' flFIRST QBL TIME D, 13 *□: 0.12 ,0: 0,12 ,0: 0.F5 . 2 ' KFMTt'k If 
QBL 

[6] * HLAST QBL TIME □.13.0:0, 12,0:0, 12, □:n.F5.2*AFtfr$-4 

QBL 

£7] * ' 

7 
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V DBARASUP \NAMEQALi PRQAL j PRT ; TIVI0 j T ORB ; PSfOP 

Cl3 « GROUND SUPPORT FUNCTION FOR PLOTTING A8AP BAR ANGLE DATA . 

C 23 'm QBLSBT FUNCTION BN BN BXBCUTBD WITH LATEST DATA ?* 

C 31 * ENTER Y OR N FOR YSS OR NO* 

C4] +C v N**ti)/LY t LN 

C 5] LI i •OBTAIN FOLLOWING PROM CIRCLE CHART* 

C 6] * INPUT MIDNIGHT GMT OF QBL SET* 

ill TMID+SPDQ 

C8] * INPUT 1 REV DELTA TIME* 

1 9] TORB+SFD O.D 

CIO] *TO GET A PLOT OF THE BAR ANGLES* 

C113 * INPUT A CARRIAGE RETURN* 

£123 PRQAL+' 4 0+flflLjQ 
[13] PRT+SFD 4 0 *QBL 

C 14 3 TSTOP+TOP PRT+PRT-PRTLll-TORB\PRTlll-TMID 
CIS] DBARA 
C 16 ] ♦O 

C 17 3 LNt ‘EXECUTE FUNCTION QBLSBT BEFORE USING THIS FUNCTION* 

V 

V DCD 

Cl] 1 DRAWN* PRCD** DRAWS DRAG COEFFICIENT VS TIME PLOT. 

7 

7 X+RD DDDDER R 

Cl] ft COMPUTES 3RD DERIVATIVE OF AR FOR INTEGRATION OF ORBIT DIN 
AMICS. 

[2] X+{‘GME*R2*RMAG)*RD-OMOK*RMAG*Bl 
C 3] X+X+{(3*GMB*R2*R2*RMAG)*R)+iDELG R)+.*RD 

V 

7 DECMG 

[1] 1 DRAWN* PRECMG** DRAWS ECMG VS TIME PLOT. 

7 

7 A+DBGA RiP 

[1] P+QPA R* COMPUTES EIGEN ANGLE FROM DIRECTION COSINES. 

C 2 ] i*+2KDFJr*(-20PC4])x(3+P)*/#j4G 3+P 

7 

7 A+DEGQ R 

Cl] « COMPUTES EIGEN ANGLE FROM INPUT QUATERNION. 

C 2] A<-2*DFRx(*20RlH])*(3*R)*MAG 3t R 

7 
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V DG+DBLG R ; GNOM t OHO l ; KTOO 

Cl] A COMPUTBS NUMERICAL GRADIENT OF GRAVITATIONAL ACC BIB RAT 10 B . 
C 2 ] DG+1 000000* (OACCBL rt*lf • 8 ,00) •GNOM+OtCCBL R 
£33 DG+DO t il AnQQQOQQxiGACCBL S+O.lf-e.OJ-OJWJIf 

[4] DQ+DQ ,1000000* (GACCSL R * 0 QtlB'*)-GB0M 

; 1 v ;;^ 

V /?«■!>*»$ B\BL\BU 

Cl] R+ 2p0* INTERPOLATES DENSITY FROM TABLE OF DBBSITI VS ALTITUD 

E 

£23 

C 3] BU+HL+1 

*C4] Hll+RBOTABtBl')**((9RHOTABlBU‘)*RBOTABlBL'))*ALrTABtBtn-ALTTABl 
BL))*B-ALTTAB{BL'} 

C 5 ] Rl 2l+BULGBTABlBL1+( ( BULGBTABlBV) -BULGBTABlBLl )*ALTTABlBtn- 

ALTTABlBL1)*B-ALTTABlBLl 

V 

V DBSAT\A\A1\A2 

Cl] * API VERSION OF TEA SEEKING MET BOD USED FOR DESATU RATIOS . 

[2] EBP* RE 

C 3] EBNP+EBM 

£43 BTVBOM+AL+ . *BTLNOM 

i 5 ] BTACS+VLK+ . *ROT* . * ( WLK ) ♦ . *ETACS 

C 6] EE+BTS+STACS-BTVNOM 

[7] EBM+EELIN LIMC EE 

£83 EPHB+0 ,S*EPHLIM LIMC SLOPE* ,*EEM+BB-EBP 

C 9] ETHB+ 0 ,5*AL* .*ETBL-KETB*ETHLIMl ( ALL j 3 ] ♦ . *SL0PE * . * EBMP+EB- EBP) 

f -ETBLIM 
CIO] ETHL * 0 0 0 

C 1 1 3 EQAL+EPT t ( 1 - EPT * . * EPT+EPHB+ETBB ) *0 . 5 

[12] VL+AL+AFQ QAL+iQDB EQAL)*,*QAL 

[13] BTACS+0 0 0 

[14] EQVBK+EPBH . ( 1 - EPBH* . * EPHH ) *0 . 5 

[15] VBK+AFQ QVEK+(QDB EQVEK) * . *QVEK 

[16] VLK+VEK*,**LEK 

[17] a PRINTOUT 

V 

V Y+DET X\I\K\V 

[1] I+(pX)CK+l] 

[2] +U[*iK] = 0)/7 

[3] K*(^po).(4>(I-X)f4>JrC;^])UCAiK] 

[4] X+X-Vo .xXCKil 

[5] *(I>X*X* 1 ) / 2 

[6] +0tX++*/l 1*X 

[73 X+XlK*\I-KiK]*0 

[8] ♦(0»*/r)/10 

[9] ♦3,I[**>I[JCi]-(y/Cl]ZC^iI-IT;])Clj] 

Cio] y+Q 

[11] a COMPUTES DETERMINANT OF ARBITRARY SQUARE MATRIX . 

7 
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7 B+DFDOY A 

ill A COMPUTES DBMS FROM DOT IB DECIMAL FORM . 

£2] 5+365 24 60 60T86400M 

7 

V R+DFBBX X 

£1] J?+ 16 16 16 161*4+0 0 0 0,Xa CONVERTS BEX DATA 1BT0 DECIMAL 

7 

7 R+DFOCT X 

ill £*-8 8 8 81*4+0 0 0 0,Xa DECIMAL FROM OCTAL CONVERSION. 

7 

7 1+DPS X 

£1] 1+365 24 60 60tXa COBVBRTS 8BC0BDS IBTO DBMS. 

7 

7 DB 

[1] 1 DRAWN * OBK+ . *PRRTBK * A DRAWS BTSK IN 0 FRAME VS TIME. 

7 

7 DOY+DOYFALT ALT\XI\Y1\I 

Cl] a INTERPOLATES DOI VS ALT TABLE TO ESTIMATE DOT FROM GIVEN A 
LT. 

C 2] YI<—TABALTBM 
£3] XI+TADDOY 

£4] I+lH + /ALT*.STABALTBM)\.lQ 

£5] m+xxcjr]+((«xir)-iici])xujciii]-xici])+riciii]-iici3 

7 

7 DPB 

111 1 DRAWN ' PRLPH ' A DRAWS LPB VS TIME. 

7 

7 DPWR 

Cl] 1 DRAWN* PRBRKil\l'* DRAWS SOLAR ARRAY ENERGY VS TIME PLOT. 

7 

7 DPWRAiPRT 

til A DRAWS POWER ANGLE VS TIME FOR SIMULATED ASAP DATA. 

C2] 1 DRAWN , PRPWRAt2il * \PRT+PRPWRAil\l 

7 
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7 DR 

f 13 ft MENU OP PLOTS TO 00 THROUGH TO A8SBS8 8IMULATI0E RUH . 

[23 DTQ 
C 33 DHi Q 
[43 XBtfl 
[S3 Prf}0 
[63 XA ;□ 

[73 PPifiD 
[83 DROin 
[93 DPWR {□ 

7 

7 0+1 DRAW fi;£ 

[1] ♦((V/£«*((A/1 4*p£)v(l*ppff)AA/~l 2 3cpfi),(A/l 2*pp£) t ( l*P«fi) V 
£a~a/</ 2 2p£) , (~£<*’lAp£«-,£) »l)/4 3 2 t fl«*t0 

[23 2 LGRF L 

[33 +0,pD<-ff 

[4] 60 * ERRORDRAW \,<3 Zp'LERGTB RAHK DOHAXH * ) [fill* 3 

[5] ft GRAPHICS FOR DRAWIBG 2 OR Z DIMER SIOBAL PLOTS. 

7 

V P DRAWS NANEiLiRiIiNiNIiSiA;n 

[13 ft LABELS , SCALES , MP PfliWS MULTIPLE PLOTS OH A PAGE. 

[2] A+iNAME 

[ 33 2>*l + pd 

[43 

[53 N+(qA )*I 
[63 ei+ii*p 
[73 P<*-(P*i#I)-P-l 

[83 STRAPIS 5 
[93 DATE 
[103 2T21 21 

[113 ' FISAL VALUE OF ' ,HAME, * • ' iAU*\H) 

L 123 WISDOW 0.5 7.75 0.5 5.5 
[133 BOX 

[143 SCALE(-TSTOP*Q.l) ATOP TSTOP) ,(-11) ,11+TOPUR ,I)p4)[jS3 
[153 

[163 11:0 PPvlV >J[5375 PRTiS 3 

[173 + 

[18] +(Q*pA) / LI 

V 


7 P DR AMI 8 RANK \ Li R\ I it i HI i III At 8 

Cl] A+tHAMB* PLOTS BAKED VARIABLE(S) ¥8 TIKE OB SIBOLE PLOT IB S 
HALL WIBDOW 
C23 I1«*0P 4 

CS3 X«-Mtp4 
C4] 4«*,4 

C5] I«*(p4)»X 

C6] BI+fI*P 
C7] 5*(PKiir>-p-i 
c«] ST R API 8 5 

[9] DATE 

CIO] * PIBAL VALUE OP M4M,' > * 

Cll] ,U,l)p4CX*i#] 

C 12] • * 

CIS] WIH DOB 3.8 7.8 0.2S 4.25 
C 14 ] BOX 

Cis] 5<?4£ff<-75r0Pxo.i) t (rdP r5P0P),(-ii) t ii 
C 16] AXE8 

C 17] XI 1 0 DRAW AlSZVS PRTiSl 
C 18] 4+X+4 

C 19] ♦(0*p4)/Xl 

V 

7 D+DRIFT 

Cl] a COMPUTES ACTUAL 8TRAPD0WB ERROR FROM SIMULATIOB DATA. 

C 2 ] D+-0EK+ . « ( DEO A ( MV LX )*."VL) + .*LEK 

7 

7 DRO 

Cl] A DRAWS 4xrxr00f VS TIME WITH IBITIAL ALTITUDE AS REFEREBCE. 
C 2 ] 1 DRAWS' -REA+ALTIC-MAQ PRR * 

7 

7 070 

Cl] 1 DRAWS * PRTQGK+PRTAEK * A DRAWS TOTAL EXTERBAL TORQUE VS TIME . 
7 
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7 iA+I EAFA AtU$E 

Cl] *TBI8 PBOOBAM COMPUTES BULK* AM OLE 8 FBOH DIB C08IBB HATBICE8. 

[2] ft I BEADS FBOH LEFT TO BIDET TEE B0TATI0B8 FBOH TEE OLD TO TE 

E MEM. 

(93 ft!4(l3 18 ABOUT AXI8 I (13 ETC . 

(43 tPBOGBAMNBD El IOBB 0LAE8E 1 - 29 - 79 . 

C$3 0+9 3pl 0 0 0 

C$3 +UCl3*J(93)/L 

C 73 S+DET Ut |J(13 *1(23 ./C$33 
[93 I4«-105*4Cr(93irCl33 

C 9 3 f4*((-5«4ClCl3iJC23 3)4741 4(1(93*1(933) *»4 

C 10 3 I4*#4 , (-5»4CIC23 1 XCl3 3 )4f4f 4(lCl3ilCl33 
C 11 3 *0 

C 123 *tlC93«-6-I(l3^I(23 

C 193 f4«-204(lCi3ilCl33 
(143 S+DBT aCtlCl3,IC23»lC933 
(153 BA+BA, AtIt2Ultiy}ATAM S»4 CXC 9 3 1 XC 1 3 3 
(163 *4*<4ClCl3|l(2334r4#-£x4(.rfl3llC933).f4 
? 

7 EDITQBLiI 

(13 ft ALL0V8 OBOUBD COETBOLLEB8 TO EDIT DEL ABBA I AFTEB IT EAE E 
EBB BUTS BED, 

(23 *f#r«* 4, I OB D TO ADD* IE8EBT OB DELETE A COLUBE FBOB DEL * 

(93 +('4.n>'<O)/£4,ir,£0 

(43 LA t * IB PUT BSV QUATEBBIOB ABD OUT* 

(53 QBL+QBL, 0 

(63 *0 

(73 LI t * EBTBR BUBBBB OF COLUHBS ABE AD OF LIFE TO EE IB SEATED * 

(•3 QBL+U+QMllQBL 

(93 'BBTBB BSV QUATEBBIOB ABD OUT * 

(103 QBL+QBL, D 
(113 W+< -!♦!)♦( 23 W 
(123 +0 

(193 LDt' BBTBB BUHBEB OF COLUHB TO BE DELETED * 

(143 W*<(9,I-1)4W),(0,I*0)*W 
7 
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9 Y+B EIOYEC 4|ftI|XtI0(70tC0tV|£|ir 
C 13 * COMPUTES SIOSSnCTOE FROM GIVES MAT SIX AES El OEM VALUE. 

C 23 X*<4,I)pi,l0X*O«4*1704 

C 93 X0*0.002"*301+f#pl001 

03 *4I4£«t4(23«0 

C S3 C0#F£fXirO*O.OO2 — 5Ol+ffplOOl 

(63 C0*( (#*4(23*2)++/ ,4*4»4-4(l3*I)*0,5 

(73 #*|<X«4(23*2>+4+.«4*4-X»lf9*C0 

(•3 lf+4( 23 "If, (4,0) p£*4+ . *M 

[ o 3 loopi t r* , xmif < ( o . «xo ) -*♦ . «xo ) , (** . «xo >♦£♦ . « xo 

(103 i*x+i,opxo*4fr,opro*i*r,op»*<+/<r-xo,ro>*2)*o.s 

(113 *((X>100) ,(4*0.0001) , l)/F4X£,£00Pl,FX4X5fl 

(123 F4X£i*O,,r*(2,4>0O 

( 113 FX4X54li*0,,r*(2,4)pr 

(1*3 4t4£****4-X*O.OOOl"0O«-<+/,A"4*4-4(l3"X>*O.S 

(153 £00P2ix*x+i.opxo*F,O0iwixo-r*<«r(iftir3>«r9<*/r«y*if*.«xo>* 
o.s 

(163 ♦( (X>100) , (4*14*6 ) ,1 )/F4X£,£00F2»FX4X542 

(173 FX4X542ir*(2,4)plM0O 

7 

V 4*440 Ell 

(13 7*16* E8C0D8S DECIMAL DATA I8T0 UPLI8X OCTAL . 

(23 4*4«2*X-S 

03 *(4*0)/I 

(*3 4*4+2*X 

(S3 £i4*(Xp2>TlO. *9999+4 

(63 4*(6p9>T(16p2)i(ll0O,6pl>\4 

(73 4*50000 50000+10 10 10l%2 3p4 

9 

9 g*A44SA 

(13 4*6 A GRP 1* ERASES PAGE FOR GRAPEIC8 FU8CTI088 . 

9 

9 4*7 SUL A 

(13 *<X«1 2 3 )/£l .12 ,£3* COR STRUCT S EULER ROTATIOB MATRIX ABOUT 

AXIS I » 1,2,3. 

(23 Llt+QiR+EULl A 

03 £2t+O»4*I0L2 A 

(*3 £3:+OtW0£3 A 

9 

9 R+EUL1 A 

(13 4*3 3 P 0 a COMPUTES EULER R0TATI08 MATRIX ABOUT X AXI8. 

(23 4(lil3*l 

03 4(2 i 23*4(8 1 33*204 

(*3 4(3;23*-4(2|33*104 

9 
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7 R+BUL 2 A 

[13 R + 3 3pQa COMPUTER EULER ROT AT I OR MATRIX ABOUT I AXIS. 

[2] A[2 ;2]«*1 

[3] tf[3;3]«-J?Cljl]«-2ojl 

[4] /? C 1 1 3 ]«—/?[ 3 * i3«-*04 

7 

V R+BULZ A 

Cl] R + 3 3 pO A COMPUTES EULER ROTATIOR MATRIX ABOUT Z AXIS. 

[2] A[3;3]«-l 

[3] fl[l;l]«-J?[2«?>20A 

[4] fi[2 ; 1 ]«--]?[ 1 , 0 J+lcA 

7 

7 A+FCF Y}X t CliSliCX;SXiKtQiB 

[13 M. Sx-i + pfA COMPUTES FOURIER COEFFICIENTS FOR DATA X. 

[23 A>o(0,i2x/m#+0.5 

[3] 51«-10£ ; C1+20AT 

[4] CA+H-SA«-(l + 2xA)p0 

[5] tf«-0;yi<-2 OpO 

[6] L:A<-A.(I+.*CK) ,-r+.xS* 

[73 G«-(CtfxCl)-StfxSl 

[83 SK+(C1*SK) +S1 *CK 

[ 9 ] CK+Q 

[10] +(N*K+K+1)/L 

[ 11 ] A4-bTOUTx(+/LnA*A)*0.5 

[12] A<-((*LTOUTxl + 2*N)*0,\N) . [0.1]i« 

7 

V FILLTABSihNiRTS 

[1] A D575L0P5 DATA TO FILL STANDARDIZED TABLES OF TEA ATTITUDE. 

BIAS . ETC. 

[2] A 75 TIME. 

[3] I«-l ; TABSL+O 3 3pO; TABSLRTS+TABEAT+Z OpO 

[4] N+qTABALTSU 

[5] L: ALT+ALTIC+KFMxTABALTSUIU+I] 

[ 6 ] IC 

[7] TABSL+TABSL.lllSLOPE 

[8] RTS+QRALG SLOPE 

[ 9 ] TABSLRTS+TABSLRTS .RTSi kRTSL ; 1 ] ; 1 3 

[10] T ABE AT ^T ABE AT . .BAR QAL 

[11] +(A2l+I+l)/L 

V 


7 A+FINDQR HihNiSiLAIiMUhASiQhlDNiXKiXKO 
Cl] I+liIDN+Z 3pl 0 0 Oa PERFORMS QR TRANSFORMATION IS QR ALGORI 
TRM 

£2] N+l+pB+A* ,*A+H 

C 3] +(R*2)/0iB+B-iA*2*TRll1)-(+/TR*TR)x(S t S)f>l t S<>0 

C 4] £+<JM#*2)pO) t [l]<2,ff+2)pOsi<«-(4,(rf.2)0O) t [l](2,*+2)0O 

C 5] Li + iEPiUiAS+BlS+I+O 1 2»I1 ) /LPP 

[6] JXO«-((*/d<S*2)*0.5 )*(4SCl]*0)-4SCl]<0 

C 7] MUI+(l*AS)*XKO+ASll1 

C 8] LAI+2* 1++ / MU I*MU I 

t 9 ] QI+IDN-LAI*i\ t MUI)*.xl t MUI 

CIO] i«C;S>iflC ;$]■►. *$1 

Cll] dCS» 1+QI+ . *d[S; ] 

C 12] flCS;]^I+.xffC5s3 
C 13 ] LPPi+(I<Z) / LP 
C 14] Ail\ iI-2]«-0 
CIS] LPi+(N>I+I+l)/L 
C 16] AlNi\N-2l+0 
Cl7] A+ m 2 m 2+A 
7 

7 RX+IN FIX X 

Cl] A THIS FUNCTION SUBSTITUTES CHARACTERS IN CHARACTER VECTORS. 

C2 ] Al#Cl] IS THE OLD CHARACTER AND IN C2] IS THE NEW. 

C 3] A B.G. •— * FIX X SUBSTITUTES AN API MINUS (•) FOR THE SUBTRAC 

T 

[4] A OPERATOR MINUS (-) FOR ALL OCCURRENCES OF - IN X. 

[51 A CODED BY JOHN R GLAESE , 1/15/79 

c 6 ] L:+(~i/v[i]eff;r«-;r)/o 

C 7 ] ->£;;rUimi]]<-m2] 


7 


V FORMATS 


Cl] 

A THIS FUNCTION SIMULATES THE 

OUTPUTS ON ASCO FORMAT 8. 


C 2] 

**50 % QAM/ ATM/ ATMDC / / 

0.3*1 

' LFMT 1 SaSTATI0BiSTATI0BTABi2i 


1]{] 





[3] 

*Q 



GMT 

0.73 


,0:0.12 ,Q:0»72,Q:Q,F7.4 

* LFMT 1 

4p DFS GMTQ+T 



[4] 

*0 



TMID 

0.73 


.0:0*72 .0:0.12 .0:0. *7. 4 

* LFMT 1 

4 pDFS TGMTO+TMID • 

•TORB 


C 5 ] 

‘0 



DMID 

0.73 


, 0:0. 72 , 0:0. 72. 0:0. F7 . 4 

* LFMT 1 

4 pDFS TORB 



[6] 

*0 

TEA CONTROL PARAMETERS 

TGMT 

0.73 


.0:0.72,0:0.72 ,0:fl,F7.4 

* LFMT 1 

4p DFS TGMTQ+T 



[7] 

* * 

Y 

Z 

DELTA 

T TE 


A * jO 





[8] 

‘0 DELTA E Q.3F10.6 

* LFMT 1 

ZqLE 



[9] 

*0 DELTA BP 0.3F10.6 

* LFMT 1 

3pbEP 



[10] 

t • 





[11] 

*0 DELTA PH 0.3 FI 0.6 

* LFMT 1 

3p2 xLPBH 



[12] 

'ttDBLTA TH 0.3*10.6 

* LFMT 1 

3p2 xDFRxLTBH 



[13] 

« • 


■ 



[14] 

*0 ETSF 0.3F10.6 

* LFMT 1 

3p ETS 



[15] 

• • 





[16] 

*0 CMG MOM 0.3F10.0 

* LFMT 1 

3p( 100*3) xRTS-(*B)xIV+ .xVL+.xq 


, OMO , 0 





[17] 

• • 





[18] 

* CMG 2 •• 

• • % • • < 

• CMG 3 W%% * 



[19] 

* IGA IGR * 





[20] 

* OGA OGR' 





[21] 

•WH SPD ' 





[22] 

*7 1 0 * HDMIBfl ,75,0 MIBR 

,75,0 

JET 1234560* LFMT 

1 2p40,2633+ 


C 233 

[24] 

[25] 

[26] 

[27] 

[28] 


MIBSUM 

*722 ,OF0 (3*15 * LFMT 0 

I t 

*0 QBL 0.4F12.7 'LFMT 1 4pG*Zi 

*0 QVA □,4F12.7‘AFW2’ 1 4p0 0 0 1 

'QQVI 0, 4F12 . 7 * LFMT 1 4p«F* VL+ .xLBK+ .x^IKEK 


7 CQ+PPWR PiFOlFPiLC 

[1] C0«-P« 15 FUNCTION COMPUTES ANGLE OF INCIDENCE OF SUNRISE . 

[ 2] 7:F0«-(0. 308 x<70xC0xC0) + (-0.924*C0*(70) + (1.924x<70) -P+0.308 

[3] FP«-(0.924xCOx£70) + ( -1.848x^0) +1.924 

[4] bC+-FO*FP 

[5] CO+CO+LC 

[6] + (0.0001< U<7)/7 
7 


7 A+GACCEL R 

[1] a COMPUTES THE ACCELERATION OF GRAVITY VS R. 

[2] ER+R*RMAG+(R2+R+.xR)*0 .5 

[3] A*-(GME*J2*R2xR2 )x(Ftfx(7.5xFtf[2]*2 ) -1 . 5) “3xtftf [ 2] x 0 1 0 

[4] KTGG+3*OMOKxOMOK+(GME*R2xRMAG)*0 . 5 


7 
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v GBP B 

Cl] a C OMPUTES GIRO BIASES REQUIRED TO C ARCS h OBSERVED DRIFT. 
C 2] B+VLtt2l*RPD*B*T0RB 
C 3] B+B* O.S 1 0.5*0.1*RPD*2**10 

C 4] GBIAS+1 1 l«.*5**-0.5 1 0.5*la+0.5 
C 5] 

[6] 'TRUE BIAS ABOUT ORBITAL I ' \OMBIAS+ ,*VL*TORB*DFR 

V 

V GBPAiBthPURAiS 

Cl] A COMPUTES REQUIRED GIRO BIASES FROM DATA III SIMULATIOR 

C 2] A USIRG POWER ABODE OR SUB SEBSOR ABGLE DATA . 

C 3 ] I+PRPWRA Cl;] i*l * PRPWRA Cl;]; PWRA+* 1 1 PRPWRA C 2 ; ] 

[4] JV(I-l),/ 

C 5] -►(IC1]>0)/£1 

C 6 ] ' IBSUFFICIEBT POWER DATA TO C BARGE BIAS' 

C 7] +L 2 

C 8 ] Lis B+VL t ; 2 ] *RFD x ( OMBIAS+ . xVLxDFR) C 2 ] +8+i -/ PRPWRA C2;I] ) ♦-/ 

PRPWRAllxn 

C 9] 'DETBL2 UPDATE * * iDFRxLTULL2}+RFDxpWRA+SxT-PRPWRAlUll2'n 
CIO] B«-5*0.5 1 0.5*0. 1*/?P0*2**1O 

[11] GBIAS + 1 1 1 • . xB+0 . 5 1 0 . 5*LB+0 . 5 
C 12] 0MBIAS+ (O.ixRPDx 2*« 10 )*□♦« 

C 13] L2 : 'TRUE BIAS ABOUT ORBITAL I ' lOMBIAS* .xVlxTORBxDFR 

V 

? GBPASUP i T ; PRPWRA ; QAL ; QVA ; VL ; 

Cl] A GROUBD SUPPORT FUBCTIOB TO COMPUTE REQUIRED GIRO BIASES 

[2] A PROM POWER ABGLE OR SUB SBBSOR ABGLE FLIGHT DATA. 

[3] 'INPUT TIME OP DESIRED UPDATE' 

[4] r«-365 24 60 60lO 

[5] 'INPUT DMID IB RMS' 

[6] T0RB + 24 60 60lD 

[7] * INPUT LATEST POWER ANGLE AND GMT ' 

[8] PRPWRA +( 365 24 60 60J.1M) ,4[1 ] ; 4«-0 

[9] 'INPUT PREVIOUS POWER ANGLE AND ITS GMT' 

CIO] PJ?WJM«-((365 24 60 60ll+j1 ) , At 1 ] ) , [ 1 . 1 lPRPWRA ;4«*0 

[11] 'INPUT CURRENT AVERAGE BIASES TO BEARBST HALF LSBS ' 

[12] GBIASil 2 3 ; ]«-l 1 1*.*0 

[13] OMBIAS^(0.1xRPDx2**10)xGBIASLlil 

[14] 'INPUT QBL ' 

[15] QAL+Q 

[16] VL+AFQ QAL 

[17] GBPA 

V 


7 RBS+GBT SLOPE BAR\M\ VLK \R\RD\ SBBKPL ; OHQBK 
Cl] ft COMPUTES SLOPS MATRIX APPROPRIATE TO IB PUT BAR ABGLBS . 

C 2 ] VLK*- 1 2 1 APSA BARiSEEKFL+l 

£3} 0NGBK*-(2OLA2K)*0NGK*0 1 0 +.*LRK 
C4] IM/?fi»'ML2 , )*LffJC[3i] 

C5] RD*>(RBA+ALT)*0M0K*LBKllO 
C6] M+IB'S JACOBI AB VLK 
C 7 ] RSS*-(B*tT)*QM 

7 

7 GGTORQUBiUOV 3 

Cl] ft COMPUTES GRAVITY GRADIBBT TORQUE OH VEHICLE , 

[2] TGGK+KTGGxUOVS X IVK+ .*U0V3**VLKl it] 

7 

7 HBARS IRTiTiSCiTSTOP 

Cl] ft DRAWS HORIZONTAL BARS OR GRAPH AT DESIRED INTERVALS . 
t 2 ] SC+WBATSCALE 

[3] T+IHT*tSClZl*INT 

[4] TST0P*-INT*ISCW4IHT 
C 53 LtHMARK T 

C6] +{TST0P*T*‘T+INT) f L 

7 

7 HELP 

Cl] » THIS FUNCTION ABORTS ANY SUSPENDED APL FUNCTIONS AND ALLOW 
S THE USER 

C 2 3 ft TO REGAIN CONTROL AND START ANY INTERRUPTED OPERATION OVER 
• 

C 33 * REENTER FUNCTION AND START OVER * 

[4] e')SI CLEAR * 

7 


7 H+8ESSBN A;N il iR;RI iIDN ;LI ;S 
Cl] I+1*TRABSF0RMS INPUT MATRIX TO HESSENBURG FORM , 
C 2] N+ltpH+A 

C 33 LiLI+ixI) ,I+f | CZ+ itf-I;I] 

C 43 H*-HLLIiLI1 

C 5] R+HlSiIliHlI+liIliS+l+I+xN-I+l 

C 6 3 fl[Si]«-ff[S;]-/?«.x£f[r+i ; ] 

[73 J5fC}I*lM[;I+l]+fiCjS] + .x* 

[8] BlSiIl + 0 

[9] +IN>I+I* 1)/L 
7 


V R+HEXFD X 

Cl] fl«-16 16 16 16T^fl CONVERTS DECIMAL TO HEX DATA, 
7 


7 HHARK X 

Cl] 0 DRAW(X>X.X)VS 10 1\UVBATSCALB 

' • 7 ! ' 

7 Q+ROMB 

Cl] AGtfF 2 a J?0#$S CURSOR FOR GRAPHICS USB 

7 
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7 IC 

Cl] « INITIALIZES ALL PARAMETERS FOR START OF SIMULATION RUN . 
[2] TTAPB+MIBSUM+T+O 

C 3 ] OMBIAS+ATBL+BTACS+AB+ABP+ABM+ABMP+APBB+ATBB+O 0 0 

C 4] OMOT+OMOTIC 

C 5] ATD2+AT*2iBTLNOM+BTLNOMO 

C 6 3 OMO+OHOK+ ( GMB HREA +A LT+A LTIC ) * 3 ) * 0 . 5 

C 7 ] TBDM+OMO+ ( ( ( OHO *0M0 ) - ( LAD I* 2QLAZK ) *2 ) *0 . 5 ) -LADY* XOLAZK 

[83 TORB+(02)*OMO 
[93 KTGG<-3*0M0K*2 

[103 TMID‘-TMID-TORBxTMID>OiTMID‘—(OMOTIC-ETTM)*OMO 
[113 EN+{OBK+{EUL 3 LAZK)+.*BUL 2 LATK)i 2{3 

[12 3 IBGBN ; ENK+ENKIC ; USBK+(EUL3-PBZ ) + .*( 10 GAIK ) , 0 , 2OGAXK+GAYKI0+ 
0MSUN*T 

[13] LEK+(BUL 2 0M0T+BTIK+O* 1 ) ♦ ,*OBK 

[143 RMAG+MAG R-IREA+ALT) *LEK[3 ; 3 }RD+(OMOX+AOMOX*SBBKFL*0 ) x (REA* 
ALT)xLEKllil 
[153 a RANIC 
[163 +(SEEKFL*0)/L2 
[17] VLKO+SEEK VLKO 
[183 SLOPEO<~(H* AT)x$M 
[193 L2iVLK+AFQ(QDB QFA VLKO ) + . xQVLKOFF 
[20] BR-R*RMAG+(R2++/RxR)*0.S 
[213 EN+EN *MAG EN+R X RD 

[22] EX+EN X ER 

[23] QVEK+QFA VEK+VLK+ .*LEK 

[24] POWER 

[25] VLEST+VL+AL+AFQ QAL+(QDB QFA VLKO) + .xQALOFF 

[26] ETS<-ETSOFF+VL+.xETLNOM 

[27] TBD+TBDM*ALLi 2] 

[ 28 ] HV+IV+.xOMV<-VL+.xO t OMO t O 

[29] OMVK+{OMV*SF)-OMDK 

[30] RVK+IVK+.*OMVK 

[31] HCMGK+ti*ET*-ET S -HV *H 

[32] HTEK+{HTVK4-HVK+BCMGK) + .xVEK 

[33] PRIC 

[34] AB+ETS-AL*.*ETLNOM 

[35] GGTORQUE 

[36] AETORQUE\AECONSTiOMGEK+OMGKx 0 1 0 ; RROTAB+RBOTABQ+RBOSIGK* 
RBOTABl-RROTABQ 

[37] DERDD*-(G ACCEL R) + (.FAEK* MSKI ) + .xVEK 

[38] DERDDD+RD DDDDER R 

[39] TEKP+TBK+WEK) + .xTABK+TGGK 

[40] OUTPUT i ATHRAW+0 

[41] STATIONTAB+STATIONTABIC 

[42] STATIONTAB+B622Q \ STATIONTAB+( (STSB-GMTO) ,0) *.x25pl 

[ 43 ] STATIONTAB+STATIONTABL ; ASTATIONTABtl ; ] 3 

[44] GBIAS + 3 3 pO 
[453 SLOPE+SLOPE 0 
[46] I LEG * 1 

V 


7 IBGBB\U8liUS2\U89\SBBiCBB 


Cll 

a COMPUTES CJf] FROM GIVEM SIMULATIOB DATA . 


C 2] 

US3+USBK 


C 3] 

USl+USUMAG U81+EB X U8BK 


CO 

U82+USZ X US 1 


C 5] 

CBB+US2+ . *BB t BBTA+DFR* • 1QSBB*USB* . *BB 


C6] 

Bt/m«*<--30XCl|2]*XCl|l]>*«10-X[l|3]*SBB*CBB*((Jt[lll]*2)*mi 


2]*2)*0.5 


C 7 3 

U I1K+(U81*2QBUZBK)+US2*10BUIBK 


[8] 

UI2K+USBK X UIIK 


C 9 3 

I KEK*U I IK, il"]U I2K, i0 •l'lUSEK 


CIO] 

OBK*OBK*MAG OEK* 0 10 I If 


C 11 3 

OBK+BB ,10. 110BK 


C 12] 

0BK*(EB X 0BKl2\l),ll]0BK 


C 13] 

BTXK*(USBK+ . *0BKI 1 j ] ) ATAB USBK+ . *0BXL 3 ; ] 


C 14] 

7 

DARK*- ( ( ( (REA+ALT ) *2 ) -RBA*2 ) *0 . 5 ) * RE A* ALT 


7 

RES+A JAC0BIAB VlK\LVLK\T0\T 


Cl] 

a RU MERIC ALU COMPUTES GRADIBHT OF ABRODIEAMIC 
DIEBT 

♦ GRAY IT I GRA 

C23 

a TORQUE WITH RESPECT TO SMALL ROTATIOBS ABOUT 

BODI AXES. 

C33 

TO+TRQ VLK 


C 4 3 

LVLK+EUL1 A 


C 53 

T*TRQ LVLK+ . *VLK 


C61 

RES+(*A)*T-TQ 


C 7 3 

hVLK*EUL2 A 


[8] 

T*TRQ LVLK+.*VLK 


C9] 

RBS+RBS, 11. lH*A)*T-TQ 


CIO] 

EVLK+EULZ A 


CU] 

T*TRQ LVLK+ . *VLK 


C 12] 

7 

RES*RE$,12Z{*A)*T-TQ 


7 

R+L LIMC A 


Cl] 

7 

R*A*{/1.\A*L* PROPORT IOSAL LIMIT OB LARGEST COMPOBEBT . 

7 

R*L LIMM A 


Cl] 

R*At[/i,MAG A* La PRO PORT 10 BAL LIMIT OB MAGBITUDE.B 


7 
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1 LOOP 

ill * MAIN SIMULATION DRIVER RURCTIOS . 

£23 Lit GOTO RQUE 

til AETORQUB 
£43 POWER 

C S3 DBRDDD+RD DDDDER R 

[6] DBRDP+DSRD+RD-OMOX*RMAQ*EX\DBRDD+{GACCBL R) + {PPAXB*LBXiU ]t 
MSKI ) ♦ ( PAEi*M8KI ) + . *VBK 

[7] TiK+iTABK+TOGK)+.*fRX 
[6] OMV+Vl+.*Q t OMO'Q 

[9] OMVK+ ( OMV *8F) m OMDK+OMBIAS 

CIO] QVBK+OMVK UPDATE QVBK 
[11] OMOT+OMOT+LT*OMO 

£ 123 or+r+Af 

C i 3 3 DBRD+DBRD+ LTD 2 » < 2 « DERDD ) ♦ LT*DBRDDD 

[14] * DBR+{LT*DBRD)+iLT*LTD2*DBRDD)+hT*LT*LT*(*t)*DBRDDD 

£ IS] DBR+LTD2*DBRDP+DBBD 

C 16 ] R* ( RMAG* ( BR* 2O0M0K* LT ) ♦M* 10 OMOX* LT)+DBR 
C 17 3 RD+DBRD+0M0X*BB X R-DBR 
Cl8] REP 

1 19 3 USB K+ ( BUL3-PBZ ) ♦ . x ( 10 GAIK) , 0 , 2oOAIK<-GAIKIC+OMSUB*T 

[20] POWER { BNDKP+BBDK 

[21] BBK+2SQQIBBK+ LTD2*BBDK+BBDKP - 2 * PWRDR 

[22] AETORQUB { Afi C0PST 

[23] GGTORQUS 

[24] TBKPP+TBKP 

[25] TBKP+TBK 

[26] r*x*<rAfK+r00K)+.Km: 

[ 27] BTBK+HTBK+LTD2*TBX+TBXP-PBRC*TBKP-TBKPP 
[283 BTVX+VBX+ . *BTBX 

[29] 

[30] jyCKGfrHMTJt-SKX 

[31] ATMDC 

[32] OUTPUT 

[ 33 3 +(T<TST0Pf STATIOBTABt 1 ; 1 3 ) /LI 

[34] +(TZTST0P)/0 

[35] TTAPB+TTAPB'STATIOBTABtlill 

[36] STBAMB+STAT10BTABi2\ll 

[37] STATIONTAB+Q USTATIOBTAB t STATIOSTABiill+96220 t Q 

[38] ♦ ((♦/P/?*JM[l;]<0)*0)/0 

[39] •ATTITUDE COBTROL LOST DUE TO LOSS OP POWER 1 

[40] «► 

7 
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7 M+MAG V 

Cl] N+i+/lllV*V)*O.S*CONPUTBS MAGNITUDES OF A SET OF COLUMN VBCTO 
RS. 

7 

7 MARK N 

Cl] 0 DRAWil 0 1\ - 2 *WHATSCALB)VS Nn DRAMS VERTICAL LINE ON GRAPH 
AT N . 

7 

7 NAVUP 

Cl] A SIMULATES A NAV UPDATE FOR SIMULATIONS OVER SEVERAL DAIS . 

C 2 ] 0M0 + < < ( OMOK*OMOK ) - ( LADI *20LAZK ) *2 ) *0 . 5 ) -LADI*10LAl K 

C 3] TORB+02*OMO 
7 

7 ROM 

Cl] A RETURNS RANDOMIZED PARAMETERS TO THEIR NOMINAL VALUES. 

C 2] SP+ 1 1 1 
C 3] QVLKOFF+ 0 0 0 1 
C 43 NCG+-Q 0 0 
C 5 ] RHOSIGK+RHOSIG 

[ 6 ] STSH+ 0 

[ 7 ] BTS0FP + 0 0 0 

[ 8 3 ETLNOMO+ETLNOMNOM 
7 


142 


7 R+BORM X 

Cl] A 8BLECT8 SET OF 8AHE SB APE A8 X FROM FORMAL DI8TRIBUTI0F MI 
TR 0 HEAR 

£2] a AFD 1 ST ABOARD DEVIATXOB . 

[3] JM+/(fXplOOO),-rXplOOO)*912.87|X4’(pX),S 
7 

7 R+BRMLIl xrt«;Xtrt/txoiro 

Cl] A FORMALIZES COMPLEX EIOEBFECFOR TO OBIT MAGBITUDE. 

C2 ] /♦lttX*<X«X)+XKM4X7,OpX+»tXr.0p»*O.5>«pXr 

c s] x*(2,i)p<(xo*x>*roHi),-<XKro*ic/])-r*xo«-xc/] 

C4] R+R*(+/ ,R*R)*0.5 

7 

7 R+OCTFD X 

Cl] 1?**8 8 8 8tXa C0BVERT8 DECIMAL DATA TO OCTAL . 

7 

7 ORBMARX j T 

Cl] A MARKS TIME PLOTS AT 1 ORBIT IBTERFAL8 „• 

C 2] T+TORB-OMOTIC*OMO 
C 3] LxMARK T 
C4] +{TSTOP* T+T+TORB )/L 

7 

7 OUTPUT 

Cl] a SELECTS PARAMETERS FROM SIMULATIOB TO SAFE FOR PLOTTIBG. 

[2] PRT+PRT • T 

C3] PRTGGX+PRTGQX*TQGK 

[4] PRTAEK+PRTAEK'TAEK 
C 5 ] PRBTEK+PRBTBK , BTEK 
C6] PRBCMG+PRECMG ,BCMQK*B 
C7] PROMVX+PROMVKtOMVK 

t 8 ] PRA+PRA ,OBX+,*( HEX )+.*-DBQA{ MV LX ) ♦ . *VL 
C9] PRLPR+PRLPR % DFR*2«LPBB 
CIO] PROMOK+PROMOK , ( RD+ ,*EX) *RMAQ 
C 11 ] PRCD+PRCD ,MAG AECOEFFt 4 S 6] 

[12] PRR<-PRR>R 
C 13] PRRD+PRRD f RD 
C 14] PREBK+PREBX , XXX , EBDXP 
CIS] PRQAL+PRQAL % QAL 
C 16] PRLE+PREB , LB 
7 
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1 PAR i UK 

Cl 3 ft LIST IMPORT ART 8 MULATTO* PARAMETER OAT A. 

C23 * * 

C 3} * * 

[4] ' '\RAME 

C5J 

C 6] * • 

C 7 3 'ALT » UALTi' KM, * \ALT*MFK\* EM* 

CS] * BTLROM ■ * i BTLROM 
£93 'LSLIM « * iLBLIM 
CIO] *KtTH • * iKtTB 

C 11 3 ' EPHLIM > *i EPHLIM RAD , * \DFR*EPELM\* DEO* 

C 123 ' ETHLIM * * i ETHLIM i* RAD, * \DBR*ETRLIM\ * DEO* 

C 13] 'VLKQ « *\VLKQ 
[14] * SLOPE « * i SLOPS 

C 15 3 *QAL * *iQAL 
C 16] *QVLK « *\QFA VLK 
C 17 3 'QVLKQ ■ *iQPA VLKQ 

C 18 3 *LLK * *\LLK,MAG LLK+DBOAWD + ."VLK 
C 19 3 * QALOFP * * IQALOFF 

C 20] ’ QVLKOFF * *;$7£*0FF 

C 21 3 *ETVK » * \HTVK*R 
C 22 3 ' ETIK * *;(*?£*)♦. aflZTJr** 

[23] ’ACS a *;ACG 

C 24 3 a *iOMDKi* RAD/S * \DFR*0MDK\* DEO/S * 

[25] *SF * * ;SF 

[26] ' ET ,TST0P ,0M0TIC * * UT ,TST0P,DFR*0M0TIC 

[27] *SBBKFL,IBM * • \SBEKFL ,IBM 
7 

7 P«-PFA A 

[1] P+MAG Aft COMPUTES SPHERICAL POLAR FROM RECTAROULAR REPRESERT 
ATION. 

[2] P+P,-MOA[3]*P 

[3] P+P.R12UTAN Rl 1] 

V 

7 POWER 

[1] ft THIS FUNCTION SUPPLIED BY ROY LARIER 3-2514 

[2] ft IT COMPUTES THE NORMALIZED SOLAR ARRAY POWER FOR GIVER SUR 

[3] ft POSITION IN VEHICLE REFERENCE. 

[ 4] ENDK+USVK [ 3 ] *DARK<VLKl ; 3] ♦ . *USVX+VBK+ . *USBK 
£5] BRDK+0fBRDK-Q.S0B*(l-BRDK)*3 

7 
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m 


9 PRZC 

Cl] A ZRZTZALZSE8 PLOT STORAGE ARRAZ8. 

C 23 PROMO B+PRCD+PRT+ \ 0 

C 9 ] PRLR+PRLPE+PRA+PROMri+PRTOOK+PRTARK+PRECMQ+PBBTER+l OpO 
C 9 ) PRPURA+PR8T+PRERK+2 OpO 

Cs] mo+paa+z opo 
Co] Ff04J«*9 opo 
» 

9 Q+Q PUT R 

Cl] a 4IL0F5 Off! TO PUT LABEL 8 OR ORAPE AT P08ZTZ0R &. 

C2] g*9 AGFF OpM. 9 AOff £ 

▼ 

9 PURAEOLEiZiTTiSiZltZlPiTROlTEOEiUSiUZOiSPRiCPRiALlLEiEAREi 
8RAX 

Cl) A COMPUTES ROVER AEOLE PROM PLZORT DATA . 

C 2 3 *ZEPUT URL* 

C 9] AL + □ 

CO] 'ZRPUT OF 4* 

C 53 4 L+AFQ ( 0090 ) ♦ . «* 41 

co ) • zrput on * 

C7] If«-<04I>4.*<4F0O> 

CO) U8+LEI (9) 

C9) t/S0«-4lC3») 

CIO) •ZRPUT AEOLE OF ZECZDEECE If DEGREES * 

C 113 C5f4+2Of/0KQ 

C 123 5Ff*<mi)»09OCn>*0SC33”09OC3] 

CIS) CFf*(0fCl]*mC3])-05CS]K09OCl) 

C 19) TB0E+DPR*{-8PB ATAR CPB )4<0,Ol)+l •l**101l ( (C5F4-ffSC2]"03OC2) 
) ♦ ( ( CPf «CFf ) ♦5M«5Pf ) *0 . 5 ) T • 1 
Cis) rfof^(860*<rfof<«i90)^(*3oo«rfof>i»o)«rfof 
CIO] * UPOVBR ARGLS 0.F9.2 * LFMT IfOf C 1*4 I TB Of ) 

0 

V QBLCBECKiZ 

Cl) a TBZ8 FURCTZOR DZ8PLA18 TEE CURRERT QBL SET TO CRECK FOR ER 
HORS, 

C 2) •TBE CURRERT QBL SET ZS AS F0LL0V8t* 

* C3) * TBE COLOURS ARE DZSPLAIED AS ROUS FOR CORCZSERESS* 

C 9 ) ' * 

CS) /♦i-ltpOW 

* CO) •l9,9F13.7,I7,0*0,I2,QiD,I2,0i0.F5.2 , AFmi,Cl]0W 

V 
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. / 

/ 


9 QBLSST { PQALL t PQALI \ I 

£l) A THIS FUNCTION IS USED AS AN AID IN ENTERING ASAP QBL DATA 
USBD 

C 2] A IN ESTIMATING STRAPDOWN ERRORS . 

[3] * INPUT QBL AND GMT IN DBMS FOR EACH POINT* 

t4] * 8 ENTRIES ARE EXPECTED AT EACB POINT.* 

C5] *IF AN ERROR IS ENCOUNTERED , TIPS *L+2 AND REENTER* 

C 6] PQALL* 8 OpO 

[7] PQALI* 0 

£83 *M*\B*o.PQALI*€PQALI 
£ 9 ] £f»‘<0.01<f l- + /P$4£rCl 2 3 43*2)/# 

£103 PQALL*PQALL .PQALI 

£113 * ENTER NEXT COLUMN. D TO DELETE LINE JUST ENTERED* 

£123 *Oh T TO TERMINATE INPUT* 

£133 +( *TD*€PQALI*n)/LT,LD 

£143 +L- 1 

C 15 1 LTt'QQBL MILT WITH 0,13.0 COLUMNS. ti* LFMT -lfp PQALL 
£ 16) QBL*PQALL 

£173 *D0 IOU WISH TO SAVE THIS QBL SET?* 

£18 3 * IP IBS. ENTER I. OTHERWISE JUST RETURN* 

£19) Q>/0 

£20) QBLS*QBLS .QBL 

£21) *QBLS CONTAINS QBLSi* COLUMNS * 

£22) t*)SAVE* 

£23) *0 

£24) Mi 'INVALID DATA. REENTER* 

£25) *L+ 2 

£26) LDi 'ENTER INDEX OF LINE TO BE DELETED* 

£27) PQALL*((B.I-1)*PQALL) .(O.I*U)*PQALL 
£28) *0 LINE 0.13.0 DELETED®' EFMT I 

£ 29 J *L+ 2 

V 

7 QBLSETSH 1 T 1 Q 

£ 1 ) a THIS FUNCTION DOES SAME AS QBLSET WITH ABBREVIATED INPUT. 
£2) * ENTER INITIAL TIME IN D H M S* 

£3) T«-3OO*5P0O 

£4) QBL* 4 OpO 

£5) * ENTER FIRST THREE COMPONENTS OF QBL* 

£6) * ENTER 0 WHEN FINISHED* 

£7) LI i*(0*+/Q*Q)/L2 

£8) QBI*QBL.Q.{\-Q+.*Q)*0.% 

£ S) *L 1 

£10) L2iT*T+ZOO*\(qQBI)12] 

£ 11 j T*DFS T 
£12) QBl*M(*QBL) ,W 

V 


7 P+QDB Q 

Cl] ft FORMS Q DOUBLE BAR. 

C 2] ?«-(«,-«) CS0B3 

7 

7 Q+QDSG AiBASGiS 

Cl] C«-4pOft THIS FU ACTION COMPUTES QUATERNION FROM INPUT EIGER ARG 
LS. 

C 23 QM+2OBANG+0 ,S*RFD*MAG A 
C 3] S+IOBARG 
C 4] $[1 2 Z]+S*A*MAG A 

7 

7 P+QDT Q 

Cl] ft FORMS Q DOUBLE TILDE. 

c 2 ] p+(fl,-«)cs 0 r] 

7 

7 Q+QFA AiSiTRAiI 

Cl] ft COMPUTES QUATERNIONS FROM DIRECTION COSIJNBS 

[2] *CODED BI J GLAESE 5-12-78 

C 3] $«-(4,Cl]-4C3;2],4[l;3],4C2sl3),4C2;3],4[3;l],dCl;2], TRA++/ 1 1 

C 4 3 J>( f ( 1 l*S«-S+$S)*l-rJ?yOCll 

C 5] <MSCI;] + (-I)$0 0 0,l-riM)*2x<i4-S[l;I]-rM)*0.S 
[6] fi+exfl[4]*|«C4] 

7 

7 R+QIN Q 

Cl] /?«-«x- l *1 *1 1« TAKES INVERS OF Q. 

7 

7 RTS+QRALG A\B\EP 
Cl] A DRIVER FOR QR ALGORITBM. 

C 2] ft THIS VERSION OF THE QR ALGORITHM WAS CODED BY JOHN GLAESE 
3-20-79 

C 3] A TO USB TBIS ALGORITBM YOU ALSO NEED THE FOLLOWING FUNCTION 
S : 

C 4 ] A FINDQR , HESSEN . MAG, QUADM, SPLIT. 

C 5 ] H+HESSEN A;EP+1B* 12 
[6] RTS+SPLIT B 
7 

7 Z+QUAD A;BiW 

Cl] -OxiO<V<-i4[2]-i4[l]x<-B+-jiCl]*2)*2,Z«-iO 
[2] +0 t Z+(B+W) t 0,(B-W+( \W)*0.5),0 

C 3] Z«-fl,(-l7) I V) *0 . 5 

[4] A COMPUTES ROOTS FROM QUADRATIC FACTORS. 


7 


7 R+QUADM M\D 

[1] ft COMPUTES EIGENVALUES FOR 2*2 MATRICES. 

[2] J?Wl 1*M 

X 33 R+R t (MCI; 1]*MC2; 2] )-M[l* 2]*Jf[2{ 1] 

[4] ZMFCl]*2)-4xl?[2] 

C 5 3 +(DiO)/L 

[6] R+O .SxRLU t D t Rll1 % -D+{-D) *0 .S 

C 7 3 +0 

[8] LiR*-0.S*(Ril1-D) ,0 , 0?[1 ]+/?«*( (F[l]i0)-J?[l]<0)x ( | P)*0. 5) .0 
7 

7 R+LSB QUANTIZE X 

Cl] ft QUANTIZES INPUT INTO MULTIPLES OF LSB. 

C 23 R+LSB * L 0 . S+X*LSB 
7 

7 RANICiR 

Cl] ft COMPUTES RANDOM INITIAL CONDITIONS ON SELECTED PARAMETERS 

[23 ft TO TEST OUR ABILITY TO MAINTAIN TEA CONTROL IN FACE OF 

C 33 ft UNCERTAIN DATA. 

[4] SF«- 1 1 l+0.01xl-0.02x?3pl00 

[5] QVLKOFP+RFD* 0 . 5 *(LEK+ . *USEK*NORM 1 ) + 0 . 5x (LBK+. *IKBKIU l*NORM 
1 ) + LEK+ . xIKEKl2 ; ] *NORM 1 

[6] fiyLKOFF^AFOFF, ( l- + /fiFM0FF*2 ) *0 . 5 
C 7 3 ACG«-2 1 1x0.0254x1 S*NORM 111 

C 8 3 RHOSIGMORM 1 

[9] ETSQFF-*-0 .1x1-0. 002x?1000 1000 1000 
V 


V F4A/P/1F 


Cl] 

ft 

DISPLAYS VALUES 

C 2 ] 

’□ 

SF 

C 33 

*□ 

CP 

C 4 ] 

*□ 

OFF 

C5] 

*□ 

ETS OFFSET 

C 6 ] 

*0 

SUNSPOT SIGMA 


OF RANDOMIZED PARAMETERS. 

= 0, 3F8 . 2 * LFMT 1 3p(100xSF-l) 

= □ » 3 F8 . 2 ' LFMT 1 3p-ACG 

= CI.3F8.2' AFM!T 1 3p DEGQ QVLKOFF 

= 0, 3F8 . 2 ’ AFMT 1 ZqETSOFF 

= □, F8 . 2 ’ AFMT RR0SIGK 


V 
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V RDTAPBiliTTiSiIliIlPiTaOiTBOEiUSiUZOiSPaiCPBiALiLBiMAMSiSMAX 
Cl] a SIMULATES READIES ASAP TAPE DATA AMD COMPUTIBQ SUE SEE 80R 
OR 

C2] A POWER ABGLE DATA AS WELL AS QBL DATA . 

C 3] I++/TTAPB*T-l*OOiEAMB+l Ep'POVER' iI*P*0 
C 4] ♦(I>0)/L1 

C 5 ] • EO TAPS DATA AVAILABLE ' 

£63 • * 

£73 *Lk 

£83 LUTT+TTAPBin 
C9] I++/PRT*TT-TORB 
CIO] ■►( J>0) / LZ 

C 11 3 * BO TAPE DATA AVAILABLE 1ST' 

£123 +Lk 

£133 LZiI+I % ILBG-+fPRTiTT 

C 14 3 $+PP*0msI[l]*O.ilC2]-I[l]]jPflP/JrMmiCl]*O.tlC2]-I[l]],[ 
0 .1] (BAR PP^LCilCl]+0.iJC2]-lCl]])Cl5] 

£ 15 ] IHtO{M-l 
£ 16 3 Il*-<Il-l) + <((Il-l)$S)i*0)il 
C 17 ] +((Il<pS)xIl>l)/£2 

CIS] 'BO SUE RISE DETECTED * 

£193 +LH 

C 20 3 L2 s TB0+0 ; CSRA+FPWR S[I1] 

£213 IS^PPPC ;ICl]+Il-l] 

£223 LB+LB*MAG LB 

£233 LB+UB X PRRDi ;I[lWl-l] ) .[O.lHP 
C 24] LBlUl+LEtUltMAG LEili} 

£25 3 Lff«-(L£[ls]I LEi2i ] ) .CllLff 

£26 3 03O«-/L[3;];/L«-(S0L1 /£[3] )♦. * (EUL2 /L[2] ) + .*EUL1(AL+RFD* % BAR 

PP$/L[{I[l]+Il-l])[l] 

C 27 ] US+LB+ . *(BUL3-PBZ ) ♦ . * ( 10 GAIK ) , 0 , 200/ IK+GAXKIC+OMSUN*PRTtIl 1 ] + 
11-1] 

£ 28 3 SP0*-<0S[l]x03O[l])+0S[3]x03O[3] 

£29] 0P0*— (-0S[l]x03O[3] )+05[3]x03O[l] 

[30] TBOE+DFRxi-SPB ATAB CPB )+{ OO 1) + 1 •1x*1O(0S/?/-05[2]x03O£2])*( 
(CPB*2 ) +5P0*2 ) *0 . 5 

[31] THOE<-TBOElUi\TBOBl 

[32] * 5/1 ,□ ABGLEn t PB.2 t n READ AT GMT □.I3,□s□,I2,0s^,I2,□:□,P6.2 , 
EFMT(NAME \ 1 SpTHOE t DPS GMTO+TLDT+PRTlItll+Il-1]) 

[ 33 ] PRPWRA+PRPWRA ,PRTllZ 1 ]+Jl] t TBQB-TBO 

[34] L4j-*-(l/Sm3]<0.94)/0 

[35] CSRA+USVKl 2] 

[ 36 ] LE+LEK 

[37] 03 0«-/L[2} ];/£«-( EULl /L[3 ] ) + .*(BUL2 /I[2]) + .xBULl(AL+RFD* ,BAR 
PRQALl i • 1+pPPr] ) [ 1 ] 

[ 38 ] US+LE+.xUSEK 

[39] SP0«-(0S[l]x03O[l])+0S[3]x03O[3] 

[40] <7P0«-(0S[l]x03O[3])-0S[3]x03O[l] 

[41] TBOE+DFR* ( -SPB ATAB CPB) + ( O. 01 )*l • lx» loll ( (.CSRA-USZ2 ]x03O[2] 
) ♦ ( ( CPB *CPB ) +SPB *SPB ) *0 . 5 ) [ • 1 

[42] * LIVE SUN SENSOR DATA ’ 

[ 43 ] rpo^(360xrisroF<«i80)+(»360xrffos>i80)+rpoff 

[44] *0500 SENSOR ANGLE a,P8.2'APWr TBOBllH jTflOfi] 

7 
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V REF . 

[13 a GENERATES ALL REQUIRED ROTATIONAL TRANSFORMATION^ FOR SIMU 
LATION . v , 

[23 ER+R*RMAG+(R2++/R*R)*0,S 

[33 EN+EN*MAG EN+R X RD 

[43 BX+BN X SR \ 

[53 LEK+EX ,[l3E0 ,[0. HER \ 

[63 VBK+AFQ QVEK 

[73 VLK+VBK+.**LBK 

V ■ 

1 REFRESH 

[13 a ERASES THE DISPLAY AND PUTS LATEST DATA EVALUATION BLOCK U 
P. 

[23 ERASE 
[33 TEAEVAL 
[43 RDTAPE 

V 

7 REPRICiN 

[13 a THROWS AWAY OLD DATA FROM PLOT ARRAYS SO THAT ENOUGH MEMOR 
Y IS 

[23 a IS AVAILABLE TO CONTINUE . 

[33 N++/FRT <T-2*T0RB 

[43 PRT+N+PRT 
[53 PRCD<-N\PRCD 
[63 PROMOK+N \PROMOK 
[73 PRHTEK+N*PRHTEKlN+O t N 
[83 PRECMG+N \PRECMG 
[93 PRTAEK+N+PRTAEK 
[103 PRTGGK+N \PRTGGK 
[113 PROMVK+BtPROMVK 
[123 PRA+N*PRA 
[13 3 PREPH+-N +PRLPH 
[14 3 PRENK+N +PRENK 
[153 PRR+N+PRR 

[ 16 ] PRRD+-N+PRRD 

[ 17 ] PRQAL+N*PRQAL 

[ 18 ] PR&E+N+PRLE 

[19] N++ / PRSTll ; "\<T-2*TORB \ PRST[2\ ]«*F/?5T[2; 3 — A? C 2 3 
[203 PRST+{Q ,N) iPRST 

7 
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7 RUN I BANK 

Cl] ft TO SIMPLIFY RUNNING TBS SIMULATION FOR TRAINING PURPOSES 
[2] ft THIS FUNCTION SERVES AS AN EXECUTIVE TO DRIVE ALL TBE VARI 
OUS 

C 33 ft SUBFUNCTIONS REQUIRED SO TBAT TBE USER DOES NOT BAVE TO BE 
[4] A AN APL EXPERT. 

C 5] IC 
[63 LOOP 
[7] +(T*TSTOP) /Q 

[63 RDTAPE\TBAEVAL 

[9] V* 

[10] * DO YOU VI SB TO CBANGE AN I PARAMETERS NOV OR TERMINATE RUN? ' 

[11] •IF IBS , BNTBR T t I OR F TO TERMINATE .CHANGE OS EVAL. A FUNCT 
ION ' 

[123 'IF YOU TERMINATE , YOU MAI RBSUME BY TYPING RUNON AND RETURN 

t 

[13 3 '•IF NO CHANGE IS DESIRED , BIT CARRIAGE RBTURN TO CONTINUE .' 
[14 3 * IF AN ERROR IS ENCOUNTERED * TYPE +6 AND RBTURN 

[153 * * 

[16 3 LO: + Ul<QNAMB)v((NAME‘+ t n)€ t W)*O)/03 
[173 +CTYF*£NAME)/0,L1 % L2 
[183 +2 

[193 LI: •ENTER NAME OF VARIABLE TO BE CHANGED' 

[20 3 NAMB+€\3 

[213 + ( 0* ( 1 *NAME ) € ' EABCDEFGHIJKLMNOPQRSTUVVXYZ ' ) / ERROR 

[223 * * 

[233 * ENTER NEV VALUE(S)' 

[243 * * 

[25] €NAMB.'+U' 

[26] 0#BIilS[l]«-(O.lxf?FD*2*-lO)*O.5x<?BIi4S[l;l3+GBIi4S[3;l] 

[27 3 0AfBI4S[2]+(O.lxBF0x2**lO)xGBr>«S[3{2 3 

[28] 0M£I>15[3]+(O.lxFF0x2**lO)xO.5x0BljlS[l;3]+GBIv!S[2;33 

[293 * ANY OTHER CHANGES ?' 

[30] +7 

[313 ERROR :* IMPROPER VARIABLE NAMB ' 

[32] *00 YOU WANT TO STOP OR GO ON ?' 

[33] * * 

[34] -*-7 

[353 02 : 'ENTER FUNCTION WITH INPUT OR OTHER EXPRESSION .' 

[36] * AFTER FUNCTION COMPLETION , HIT CARRIAGE RETURN TO RESUME SIM 
ULATION . ' 

[37] *IF AN ERROR IS MADE HERE , TYPE +02 AND RETURN 

[38] * YOU MAY HAVE TO EXECUTE RUNON ' 

[39] «□ 

[403 * ANY OTHER CHANGES?' i 0 

[41] +7 

[42] 03 : 'ENTRIES RESTRICTED TO T % Y OR F. ' 

[433 +00 
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9 RUBOB i NAME 

Cl] A THIS FUBCTIOB DOSS SAMS AS RUB EXCEPT THE RABDOMIZBD PAR AM 
ETERS ARB 

[23 A ROT CHARGED FROM THEIR CURRENT VALUES. 

[33 A IC 
[43 LOOP 

[53 +(T*TSTOP)/0 ~ 

[63 * RDTAPBiTEAEVAL ‘ - 

[73 * * ' 

[83 • DO YOU WISH TO CHARGE ARI PARAMETERS BOW OR TBRMIBATB RUNT' 

[93 * IF IBS t ENTER T, Y OR F TO TERMINATE , CHARGE OR BVAL. A FUBCT 

IOB ' - * 

[103 * IF YOU TBRMIBATB , YOU MAY RESUMB BY TYPIBG RUBOB AHD RETURB 

[113 'IF BO CHARGE IS DESIRED , HIT CARRIAGE RETURB TO COBTIBUB.' 

[12] 'IF AB ERROR IS EBC OUST BRED , TYPE -*-6 ABD RETURB 

[13] ‘ ‘ 

[14] LQ:+((l<pNAME)vi(RAME+ t n)£'TYF')=0)/L3 

[15] +('TYF'eRAMB)/O t Ll % L2 

[16] -*-2 

[17] LU' EBTER NAME OF VARIABLE TO BE CHARGED ' 

[18] RAMB+c Q 

[19] -►( 0 = ( 1 +NAME ) e ' EABCDEFGHIJKLMNOPQRSTUVWXYZ ' ) /ERROR 

[ 20 ] * * 

[21] 'ENTER NEW VALUERS)' 

[ 22 ] * ’ 

[23] eNAME, '+Q' 

[24] OMBIASi 1 ]«-(0 . l*RFD*2* m 10 ) *0 . S*GBIASL 1 ; 1]+GBJ4S[3 { 1 ] 

[25] OMBIASt2]*-(.Q . l*RFD*2* m 10 ) xGBIASl 3 ; 2 ] 

[26] CA/BI>l5[3]^(0.1xi?FDx2*-10)x0.5xGBI i 4S[l;3]+Crfli’i45[2i3] 

[27] 'ANY OTHER CHANGES?' 

[28] ->7 

[29] ERROR: ' IMPROPER VARIABLE NAME' 

[30] 'DO YOU WANT TO STOP OR GO ON?' 

[31] * * 

[32] ■+■7 

[33] L2: 'ENTER FUNCTION WITH INPUT OR OTHER EXPRESSION .' 

[34] 'AFTER FUNCTION COMPLETION , HIT CARRIAGE RETURB TO RESUMB SIM 
ULATION. ' 

[35] ’IF AN ERROR IS MADE HERE , TYPE +L2 AND RETURN.' 

[36] 'YOU MAY HAVE TO EXECUTE RUNON * 

[37] eD 

[38] 'ANY OTHER CHANGES ? ' ;□ 

[39] +1 

[40] LZ: 'ENTRIES RESTRICTED TO T t Y OR F.' 

[41] +L 0 
7 


V &+SCAL8 & 

Cl] *2 +v/4«*(a/ 0 4 6*pff) ,~A/</3 2p(ff*,fl) ,1,2 
C 2] *0*fl*3 LGRF & 

C 33 4*0006+' BRRORSCALB *.,(2 6p * LEHGTBDOMAIE * ) [fill i ] 

[4] 4 ALLOWS USSR TO SET SCALE FOR GRAPHICS PROGRAMS, 

V 

V RE SUSSEX VLKiTO 

Cl] 4 TRIS FUECTIOES DRIVES TBE ITERATIVE PROCESS WHICH SEEKS A 
TEA, 

C 2] 0MGBK+{2oLAZK)*0MGK*Q 1 0 +.*LEK 

C 3] M*lB-6 JACOB I AS VLK 

C 4 3 70*7*4 VLK 

CS] □*A7ff*-70lMf 

[6 ] RBS+VLK+(BUL1 A7*Cl] )+.*(BUL2 A7BC2] )+.*(BUL3 A7*[3] )*. *VLK 
C7] *<1B*8SIMG A7ff)/1 

V 

V SET ST AT I0SS ; SAME 

Cl] 4 SETS UP A SERIES OF SIMULATED GROUSD STATIOBS OVER WHICH 

C 2] 4 DATA IS AVAILABLE FROM THE SIMULATIOS SO DATA IS AVAILAB 

LE 

C 3] 4 BLSEWBBRB . THIS IS FOR GROUSD COSTROLLER TRAIBISG . 

[4] * ISPUT STATIOS AOS TIME IS HRS MISS SECS * 

[5] BAMB+U 

C 6] NAME +( 365 24 60 SOlDOI .SAMBll 2] , 0 ) ,**#*[3] 

C 7] S7A7IO*7ABI<>S7ii7I0*7ABI<7,*ilAf* 

C 8 ] *1 

V 
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V SETUP 

C 13 A SIMPLIFIES SETT IMG UP SIMUL AXIOM TO MATCH FLIGHT COMDITIOM 
S, 

[2] •INPUT LAMBDA Y* 

C 3] LAYK+RFD*U 

[4] • IMPUT GAMMA I' 

[ 5 ] GAYK+GAYKIC+RFD*n 

[6] * IMPUT GMT TGMT AMD TMID CORRESPOMDIMG TO ABOVE DATA ' 

C 73 * IMPUT GMT * 

[8] GMTOO+SFDU 

C 9 3 * IMPUT CURRBMT ALTITUDE IM MM * 

CIO] ALTIC+KFMxQ 

C ii3 * imput rcwro* 

C 123 TGMTO+SFDU 

[13] * IMPUT TMID * 

[14] TMID+SFDU 

[15] TMID+TMID+GMTOQ-TGMTO 

[16] * INPUT REGRESSION RATE DEG /DAI* 

[17] LADY+RFD*U* 86400 

[18] * IMPUT DESIRED GMT * 

[19] GMTO+SFDU 

[20] LAYK+LAYK+LADYxGMTO-GMTQO 

[21] GAYKIC<-GAYK*-GAYKIC+OMSUN*GMTO-GMTOO 

V 

V R+SFD DBMS 

[1] A COMPUTES TIME IN SECONDS FROM IMPUT DBMS . 

[2] £<-365 24 60 6Ol0£AfS 
7 
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V SKILAB ARAMEihA 

Cl] a DRAM 8 I DR All It D SKILAB VESICLE IE STEREO PROJECT 10 RS . 

C 2.3 STRAPJS 5 

£31 "SCALE $C 

C 43 DATBiA+MtAMAME 

£ S3 STEREO XAXZ8 

£6] STEREO LETTERX+3 6p(6p6),12p0 

£7] STEREO TAXIS 

£83 STEREO LETTERI+B 6p(6p0) , (6p6) ,6p0 

£ 93 8TERE0 ZAXZS 

£10] STEREO LETTERl+3 6p(12p0),6p6 

£11] STEREO A+,*PXERD 

£12] STEREO A+.*MXBRD 

C 13] 1+1 

£14] Lx8TERE0 A+.*(BULX l*0*k)+ ,*SIDE 1 

C 153 ♦(nw*D/i 

£16] STEREO A*.* ATM 
£17] STEREO A+.*BASE 
£18] STEREO A+.*OWS 
£193 f * 

£ 20 ] 'EQUILIBRIUM ATTITUDE ', ARAMS i* A 
£ 21 ] 1 ’ 

£22] 'BARK Alt GIB , ARGLE OP ATTACK, ROLL AROLE • '{DPR* 1 2 1 
A 

V 

7 X5L SLUPDATE IiSiMiRiSL 

£1] 4 SLOP! MATRIX UPDATE DATA IR PROPER PORMAT POR UP 

LIRE. 

£2] DATE 

£3] 'DWW1,2,3 * 0.316.2*0 DOT « 0.16.2,0 ALT • O.F5.0.0 BETA 

> 0 , f 6 . 2 ' LPMT 1 SpTABBIASlhl ,TA8D0IM ,TABALTRMtI) ,TABBBTAt 
I] 

£4] ’ # 

£5] * KSL « 'iKSL 

£6] * SLOPE MATRIX UPDATE DCS COMMARD 

S' 

£7] ' * 

£8] * 52002 UPDATED SLOPS' 

£93 * * 

£10] • ’ 

£113 SL+TABSLIIh)*KSL 
£12] $W10.5«-1000k5L)*1000 

£ 13] M+l 
£14] LMxR+X 
£15] LRxS+BRC SLCMtRJ 

£16] *129,15,0 SLOPSUJ, 11,0, 0,11,0) > 0,P6.3' 6PMT(S[lhM \R i 

SLlMiRD 

£ 17] *129,15 'LPMT 5C2] 

£18] * * 

£19] ♦ (3*f«-f*l)/IF 

£20] +(3*M+M+X)/LM 

7 
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7 SPZH KiXiI 

Cl] 1*1 A TBIS FURCTIOB ADJUSTS TBS SSBD FOB TBS RABDOH NUMBER OS 

BBRATOB 

C 23 X*T 1000000 

C33 *(RU*I* 1)/2 

V 

7 R+SPLIT B\9 \LD\I\TR\ TRT ; TRP ; I ID 

Cl] A SPLITS OFF LINSAR OR QUADRATIC FACTORS FOR TBS QR ALOORITB 

M. 

C2] R* 0 2p0 

C 3] T0Pt+( 0 1 2«f )/0,L3.M;I^0;ir^(pff)Cl];r^r/?r^0 0 

C4] LQiB*FIBDQR Bi 1*1*1 

C 5] ♦<Z>30)/L4 

C 6] Z0«-'2*l 1*1 0+fl 

C 7 ] *(BP> \ LD) / LI ,L2 

C8] ♦(</ \LD) fLLl 

C 9] TRT*TR*BlBiR),QiTRP*TRT 

CIO 3 *►(0.25 >{MA0 TR-TRP)*MAG TR)/LQ 

C ii 3 tr*o o 
C 12] +L0 

C 133 Lll\TR*2 2pQUADN *2 *2*BiTRP*TRT 
C 14 ] TRT+TR+'TRlUklTRhn-BlBtinO 
C 153 *(0 ,25>(AM<? TR-TRP)*MAG TR)/L0 

C 16 3 TR* 0 0 
C 17 ] *L0 

C 16 3 LU*TOP\B**l •UBiR*R t ll')Bl!li!n t Q 
C 193 Llt*TOP\B*'2 -2+£f;ff«-tf,Cl32 2pQUADM *2 *2iB 
C 203 Lkt 'SLOB OR HORCONVERGBHCB * 

C 21 3 'OF TBB MATRIX U' \B 
C 22 3 R*(B, 2)p0 

7 
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7 S8AM0LE\8\A 

Cl] A COMPUTES SUM SEMSOR AM OLE FROM FLIGHT DATA . 

C 23 A THIS FUMCTIOM ALLOWS IMTSRPOLATIOM BETH EMM MORE TRAM OME T 
IMS 

[3] A POIMT FOR IMCREASED ACCURACY . 

C4] 'I*P&r QUATSRMIOM DATA AS REQUIRED ' 

C 5] 'IMPUT X 8UM 8EM80R READIES' 

c 6 3 s+n 

C 73 A+SSAMGLSQ 

C 8 ] *IF IOU WISH TO IMTERPOLATE 8UM SEMSOR MMOLES,* 

C 9 ] 'BIT CARRIAGE RETURM | OTHERWISE EMTER T TO TERMIMATE * 

CIO] -►(-•fc.fl)/!/ 

C 11 3 * SUM SEMSOR AMGLB * tA 

C 123 +0 

Cl 3] LIx'IMPUT 2MD SET OF DATA ' 

C 14] 'IMPUT X SUM SEMSOR READOUT FROM FORMAT 7* 

C 153 5+5,0 
C 163 A+A,SS AMGLB 0 

C 17 ] 4+4Cl3 + (ilC2]”ACl3 )*-5Cl ]*5C2]-5Cl 1 

C 183 * IMTERPOLATBD SUM SEMSOR AMGLE ' \A 

7 

7 TH0B+SSAHGLB0iIiTTiStIliIlPiTB0iUSiU30iSPRtCPBtALiLBtMAMBi 
SMAX 

Cl] A THIS FUMCTIOM COMPUTES THE SUM SEMSOR AMGLB FROM FLIGHT DA 
TA . 

C 23 'IMPUT QBL' 

C 3] AL+ 0 

C 4] • IMPUT QVA ' 

C 53 AL+AFQ(QDBn) + .*AL 
C 6 3 'IMPUT QVI* 

C 7] Lf+<*/lZ;)*.K(jlFG0) 

C 8 ] 55+L5[;3] 

C 9 ] £/30+4L[2 ; ] 

CIO] CSRA+Q 

C 113 5P5+<05tl]*';3OCl])«>55[3]*53OC3] 

C 123 CFff+(f/5Cl]*tf30C3])-l/5C3]*530Cl] 

C 133 TBQB+DFR*(-SPH ATAM CPB)*(Q ,Ol)+l •l*«101l ( (CSRA-USl2)*U30l2l 
) *((CPH*CPB)+SPB*SPB)* 0 . 5) r-1 
C 143 TBOB+(nQ*TBOB<'lBO) + (*nO*TBOB>160)+TBOB 
C 153 TBQB+THQBtUl\THOB'} 

7 

7 STATES 

Cl] A TBIS FUMCTIOM DRAWS TBS STEREO AXES AMD LABELS THEM WITH S 
TBRBO LETTERS . 

C 2] STEREO XAXIS 

C 3] STEREO LETTER!* 3 6p(6p6),12p0 

C 4] STEREO TAXIS 

C 5 ] STEREO LETTBRI+Z 6p(6p0) ,(6p6> ,6p0 

C 6 3 STEREO ZAXIS 

C 7 ] STEREO LBTTBRZ+Z 6p(12p0),6p6 


7 
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7 STEREO PiBiXRiXLihM 

Cl] A DRAWS STBRBO PROJECTIONS OP POIBT 8BT P USING PREVIOUS LI 8 
PBC1PIBD 

C 2 ] A STS LOCATION t SBPARATIOB , LIBB OP SIGHT % BTC . 

[3] p^(fi£.ci]ae.co.i]Mfi)4..K? 

C4] ^(PC3;]-IUC33)*C 

c$] Xfr-<PCl;3*a-IX4tl3>*#»XJ?«-<P[lt3-g«IXlCl3>** 

[6] r^(PC2s3-mt23)*ir 

[7] WINDOW 3.5 8 0.5 5 

[8] 0 DRAW T VS XL 

[9] WINDOW 0 4.5 0.5 5 

[103 0 DRAW I VS XR 

V 

V STPAR 

[1] * PRINTS OUT CURRENT STBRBO PARAMBTBR8 , 


[23 * LINS OP SIOBT * ifcgfi 

[33 'OBSERVERS RIGHT ' iRl 

[4] * OBSERVERS UP * l!LZ 

[53 'AVERAGE BIB POSITION ' iElH*LG& 

[63 ' BIS SEPARATION 

[7] * PAGE DISTANCE ' jfi 

7 


7 g+STRAPIS fl 

Cl] a SE7S UP PROPER STRAPPING POR GRAPHICS FUNCTIONS AT AN INST 
ALLATION . 

[23 a >17 M5FC 17 SHOULD BE SET TO 5 AND MUST BB BXBCUTBD ONCE AT 
LEAST 

[3] A FO/? EACH TERMINAL SESSION. 

[4] ♦m*i(7<fl)v$>fl+l*fltg«-i0 

[53 ♦0;fi+ll tGRP B 

[63 ERR: 'NO SUCH STRAPPING ' 

7 

7 STROTiTRiTRliTH 

[1] a MOVES LOS AND BIS ON SURPACB OP SPHERE POR BETTER VIEWING. 

[2] 'INPUT DESIRED AZIMUTH AND ELEVATION OP LOS ' 

[3] 'OBJECT TO OBSERVER , ANGLES IN DEGREES' 

[4] TR+RFD**\ 1*$0 

[53 70«-l*FFF-£{2$ 

[63 TR1+IBUL2 THlll)+.*BUL3 7ff[23 

[7] TR+M(EUL2 TRlil)+.*BULi 7F[23 

[8] 7F«'7J?'* . *7F1 

[ 9 ] ItiS+TR+.*lC§, 

[ 10 ] m+TR+.*{lZ 

Cii] &Z+TR+ •*&! 

7 
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V 8T8ET08 

ft 8BT8 08888888 8A8AM8T888 808 888880 8 188 180. 
U8+Q\Q+'0B8B88i8 LOCATION 
Mi+0\Q+'0B8B888B 1188 08 8I08T I8t* 

gf+’D 1 1 > ' 08888888 88BTICAL DI88CTZ08 ( 8880 80T 88 EXACT)* 
M&+WHAO Mi 

ILt+1Lt-m*LU+-*1Ll 

{Lt-!Ll**AQ ii 

fi«-0.5«0»Ck*ir* 888A8ATI08 ZB 0888 0BZT8 * 

I MZIC3 3 *818-808888 DZ8TA8C8 
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V TACS BTS 

Cl] * TBIS FUNCTION SIMULATES TBS SKXLAB TAC8 FIRING LOGIC 
C2] -**(1.92 >MAG STS) /O 

C 3] FIRB+t 3pO 

C4) ♦( ( ( |5?5[2] ) > ( I fr5[3] )♦! 3 F) a( | ff5C2] )> t ffSCl] )/Ll 

C 5] +(<(<!*r5[3])+ff3fi)*!ffr5C2])v(!*mi])ft|ffW2])/£3 

C6] ♦O 

C7] Lli*<mC2]*0)/L2 

C 8 ] FIJ?B[4t]«-2 1 1 

[9] A *JBT 4* 

CIO] +1 7 

Cu] £2tm£Cit]^i i i 

[12] a VM 1* 

[13] +L7 

[14] I3^(m[l]<ffBCl])/ii5 

[15] ♦(STS[3]iO)/L4 

[16] FIRBl 6;]*1 1 1 

[17] A WET 6* 

[18] +L7 

[19] I4*FIFF[2*>1 1 1 

[20] A WET 2* 

[21] +L 7 

[22] I5^(SrS[3]*0)/£6 

[23] FIFF[3;]«-1 1 1 

[24] a * JET 3’ 

[25] +L 7 

[26] L6*FIFF[5;]M 1 1 

[27] A WET 5* 

[28] L7:HTEK+(*VEK)+.*BTVK+HTVK+LBTVK++/NIB*LBVS X MIBS**FIRB 

[29] a T * , ;2V LHTVK * * j LBTVK \ ' W0K ■ ' iMAG BT0K*B 

[ 30 ] BTACS+ETACS-LHTVK*H 

[31] BTS+BTS+ tHTVK*H 

[32] Q+MIBSUM+MIBSUN+l 

[33] a 'MIBS* * \HIBSUM 

[34] +1 
7 

7 

[1] A THIS FUNCTION PRINTS A STATION BLOCK OF DATA TO BBLP BVALU 
ATS 

[2] a TEA CONTROL PERFORMANCE . 

C3] ‘ * 

[4] 52MriGB[S27MJIfS|],,<'X7,0 GMT • 0,13,0:0, 12, 0*0, 12. 0*0. F6. 2* 
LFMT 1 4p DFS GMT0*T) 

[5] *UTIMB TIL NEXT STATION 0, 12. 0*0, 12, 0*0. F6. 2,0 PRESENT A 
irO,F6.0 *A FMT 1 4p ( 24 60 BOjSTATIONTABi 1 * 1 ]-D t MFK*ALT 

[6] *0 ELAPSED SIM TIME * 0.13 .0*0.12 ,0:0,12 .0:0, F6 . 2 * EFMT 1 4p 
DPS T 

[7? 4 0 BANK, ATTACK . ROLL ANGLE • 0.3F10. 2 * EFMTMBAR 4 •1+FJT04& 

[8] PRST+PRST , T , PFT i T 

[9] a •RQRBITAL DRIFT ANGLE D,F10.2*Am • UPRCBI 
7 
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7 TEX PERT iTitCGiSLPiAiOHGEKiRiRDi SEEKPL 
Cl] A TRIS function ESTIMATED THE CP OFFSET AMD recommended a mo 
MBSTUM BIAS 

C 2 3 a TO HELP US ACRISVE OUR DESIRED ATTITUDE . 

C 3] SEE IF 1+ 1 

C 43 OMGBK+< 20LAZK) *0HGK*LBKZ2 * ltR+(RBA+ALT)*LBKiZi ] iRD+OMOK * ( REA ♦ 

[5] LCG+-CPASUP 

[6] a * in PUT CURRENT SLOPE MATRIX * 

1 7 ] A SLP^O 

C 8 3 ETLNOM+ 0 0 0 

£93 T+TRQ AL+ 1 2 1 4F** RFD*BATNOM 
CIO] * SUGGESTED MOMENTUM BIAS * 

C ii j i o jt 

7 

V k+TOP AiLiP 

Cl] A DETERMINES RANGE OP VARIATION AND ROUNDS OFF TO CONVENIENT 
NUMBER. 

C 2] /?<»-PU<-10*ll0«P<*l7 I .A 

C 3] P«-1.5 2 2.5 3 4 5 6 7 8 9 10 

C 4] p<-pciwp*p]*i; 

7 
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[13 

C2] 

[33 

[43 

[S3 

[63 

[73 

[83 

[93 

[103 

[113 

[123 

[133 

[143 

[153 

[16] 

[173 

[183 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 

[25] 

[26] 
[27] 


TRANSLATE 

« PROVIDES CONVERSION TABLE TO CONVERT SKYLAB TERHINOLOGI TO 
a SIMULATION TERMINOLOGY . 

*ETLN I * \ETLN 1 
*BTLN 2 = \ETLN2 
* ETLN3 = * t ETLN3 
'KDETBX- \KDETRX 
* RDETNZ = * tKDETHZ 
* DBTBL1 - ' % DETBL1 
'DETBL2- \DBTUL2 
' DETBL3- * t DETBL3 
*DEPBL1- ' t DEPBLl 
* DEPBL2~ * t DEPBL2 
'DEPBL 3= ' % DEPBL 3 
•0SL1 * ' t DELl 

'DEL2 ~ ' % DEL2 
'DEL 3 = *,D5L3 

* TDESAT- * % TDESAT 
' DEPRC3L- ' t DEPBC3L 
•PHD1DX= ' t PHDlDX 
'PHD1DY= ' t PHDlDY 
' PHDlDZ- ' .PHDlDZ 
'PBD2DX- ' .PBD2DX 
* PBD2DY- * .PHD2DY 
' PHD2DZ- * % PHD2DZ 
'PBD3DX- ' .PBD3DX 
' PHD3DY- ' .PHD3DY 
'PHD3DZ= * .PHD3DZ 


7 RES+TRQ V LK\OMV K\HVK ; VEK\RHOT AB 

[1] A COMPUTES EXTERNAL TORQUES. USED IN SEEK AND JACOBIAN. 

[ 3 3 AETORQUE ; GGTORQUE ; . *LEKiAECONST 

[4] 0WA>yLtf+.xO,0M0*,O 

[5] 

[6] Jf H*VLK+ .xETLNOM 
V 


V QNEW+OMA UPDATE QOLD iQOM; LPBI \MAGOM 

[1] A QUATERNION INTEGRATION FUNCTION . 

[2] UOM+OMAiMAGOM+MAG OMA 

[3] hPHI+MAGOM*LTD2 

L ] QOM+(UOM*10LPHI) ,2obPHI 
[5] QNEW+{QDB QOM)+.*QOLD 

V 
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V £*£ VS fitCtfitff 

Cl3 "►2*v/£4.(A/£*i$£«.i,(j«.p£«.,£) ,£*l*pfl«-(*2 + l l»pfl)pfl) »(l<PP£)v2<p 

pfi 

C 2 3 *Oxpfl«»(((i*pfl),£)pff),tl]<£*ir£)pfc 

C 3 3 .(2 ZIq'AH argument of vs is of improper rare arguments of vs 

ARE NOT CONFORMABLE ')[(*£) lit 1 

C4] a TAKES 2 VECTORS OF SAME LENGTH AND CONSTRUCTS 2*lf MATRIX F 
OR PLOTTING . 

7 

7 fr-WHATSCALE 

Cl] £*11 LGRF 1A DISPLAIS CURRENT SCALE SETTING. 

7 

7 g+WRATWINDOW 

Cl] £*11 A GRP 2 a DISPLAIS CURRENT WIN DOM SETTING. 

7 

7 fl * WINDOW & 

[1] *2+v/fl«-(A/0 4*pfl),*/</2 2p(ff**ff),i4 

C 2 ] *pfl*2 LGRF fi 

C 3 3 £*0pO*6$’ ERRORWINDOW *,,(2 6p f LENGTH DOMAIN ' ) [fill; ] 

C4] A ALLOWS USER TO SET NEW WINDOW SPECIFICATION . 

7 

7 C+-A X B 

Cl] A CROSS PRODUCT 

C 2] 0((1<K13/1)*(-14>[13B)>-<14>C13B)X-I*[i3i« 

7 

7 XA 

Cl] A CONSTRUCTS CROSSPLOTS OF THE IN PLANE STRAPDOWN ERROR. 

C 23 1 XPLOT'PRA C3 It]* 

7 

7 XH 

Cl] A CROSSPLOTS IN PLANE COMPONENTS OF HTEK. 

C2] 1 XPLOT' (OEK*.*PRHTEK)lZ 1;]* 

7 
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V P XPLOT NAMB\S\N\I\NI\I1 

[I] a MAKES CROSS PLOT OP DESIRED DATA ON SQUARE WINDOW . 
C 2] ST RAP IS 5 

C 33 DATB 
C 43 I+(f>€NAMB)l2l 
C 53 NI+U*P 

[6] S+-(P* \NI) -P- 1 

[7] WINDOW 2 7 0.25 5.25 

[ 8 ] BOX 

C 9 ] scalec -jri ) ,n , ( -n ) ,n^?r/c 1 3 1 €Namb 

[10] AXES 

[II] 0 DRAW(eNAME)itSl 

[ 12 ] * 

[13] * FINAL VALUE * , NAME , * * * ; , ( ciMWP) [ ; • lfpeiMWff] 

V 
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